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PREFACE 


This  paper  contains  the  proceedings  of  the  Low-Speed  Boundary- 
Layer  Transition  Workshop:  II,  held  at  the  offices  of  The  Rand  Cor- 
poration, Santa  Monica,  California,  September  13-15,  1976.  The  Work- 
shop was  jointly  sponsored  by  the  Defense  Advanced  Research  Projects 
Agency  and  the  Office  of  Naval  Research,  and  organized  by  The  Rand 
Corporation.  The  proceedings,  which  have  been  given  a limited  circu- 
lation to  ARPA  and  ONR,  are  here  reproduced  in  Rand's  Paper  series 
for  worldwide  distribution. 

The  primary  goal  of  the  Workshop  was  to  provide  a forum  to  exchange 
current  research  results,  to  assess  recent  progress  in  laminar-flow 
technology,  and  to  establish  new,  viable  research  goals.  Consistent  with 
this  principal  theme  was  the  desire  to  provide  an  opportunity  for 
researchers  and  engineering  designers  to  interact  and  to  develop  a 
mutual  appreciation  of  their  respective  efforts  in  improving  the  knowl- 
edge of  transition.  This  paper  was  prepared  from  the  camera-ready 
abstracts  provided  by  the  authors,  and  should  be  useful  to  hydrodynamicists , 
designers  of  submersibles , and  others  engaged  in  fluid  mechanics  research. 
Other  related  Rand  publications  include: 

R-1752-ARPA/ONR,  Low-Speed  Boundary -Layer  Transition  Workshop , 

June  1975. 

R-1789-ARPA,  Controlling  the  Separation  of  Laminar  Boundary 
Layers  in  Water:  Heating  and  Suction , September  1975. 

R-1863-ARPA,  The  Effects  of  Wall  Temperature  and  Suction  on 
Laminar  Boundary -Layer  Stability,  April  1976. 

p 

R-1898-ARPA,  "e  Stability  Theory  and  Boundary -Layer  Tran- 
sition, February  1977. 

R-1966-ARPA,  The  Buoyancy  and  Variable  Viscosity  Effects  on  a 
Water  Laminar  Boundary  Layer  Along  a Heated  Longitudinal 
Horizontal  Cylinder,  February  1977. 

R-2111-ARPA,  Entry  Flow  in  a Heated  Tube,  June  1977. 

R-2165-ARPA,  Approximate  Methods  for  Calculating  the  Properties 
of  Laminar  Boundary  Layers  in  Water  (to  be  published). 
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INTRODUCTION 

The  Low-Speed  Boundary-Layer  Transition  Workshop  was  held  on  Sep- 
tember 13  to  September  15,  1976  at  the  Rand  Corporation,  Santa  Monica. 
The  purpose  of  the  Workshop  was  to  provide  a forum  to  exchange  current 
research  results,  to  assess  recent  progress  in  laminar -flow  technology 
and  to  establish  viable  new  research  goals.  In  support  of  these  goals, 
one  specific  aim  was  to  provide  an  opportunity  for  researchers  and  en- 
gineering designers  to  interact  and  to  develop  a mutual  appreciation 
of  their  respective  efforts  in  improving  our  knowledge  of  transition. 

Research  into  the  fundamental  aspects  of  boundary-layer  transition 
has  advanced  relatively  rapidly  during  recent  years.  Improved  experi- 
mental protocol  is  providing  a wealth  of  data  that  both  substantiates 
and  even  surpasses  available  theory.  The  catalyst  for  the  rapid  devel- 
opment in  theoretical  research  has  been  improved  analytical  techniques 
and  faster,  larger  computers.  The  large  number  of  research  papers  pre- 
sented at  the  Workshop  attest  to  the  vigor  of  the  current  revival  in 
transition  research. 

Concurrent  with  advances  in  research  have  been  advances  in  the  en- 
gineering art  of  applied  laminar-flow  technology.  A principal  reason 
for  this  renewed  interest  is  the  reduced  availability  of  energy.  Con- 
sequently the  desire  to  improve  performance  and  reduce  energy  consump- 
tion for  both  civilian  and  military  applications  has  inspired  design  and 
development  studies  of  long-range  aircraft,  high  speed  ground  transpor- 
tation, and  high  performance  underwater  vehicles.  Designs  are  available 
that  would  exploit  laminar  flow  technology,  and  performance  estimates 
have  been  made.  However,  this  still  appears  to  be  an  area  of  highly 
specialized  expertise  and  intuition.  What  is  needed  is  a routine  pro- 
cedure for  designing  laminar-flow  vehicles  and  estimating  their  perfor- 
mance. Several  such  procedures  were  introduced  at  the  Workshop. 

The  Workshop  was  the  second  low-speed  boundary-layer  transition 
workshop  with  joint  sponsorship  of  ARPA  and  ONR  and  organized  by  Rand. 
Attendance  was  by  invitation,  and  there  were  110  participants  and  32 
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invited  papers  were  presented  in  the  six  sessions.  The  papers  were 
divided  according  to  subject  matter  and  this  resulted  in  the  follow- 
ing sessions: 

I.  Numerical  Investigations  of  Navier-Stokes  Equation 

II.  Transition  and  Stability  Theory 

III.  Predictions  of  Transition  by  Model  Equations 

IV.  The  Effects  of  Free  Stream  Disturbances  and 
Surface  Roughness  on  Transition 

V.  Experimental  Investigations  of  Transition 

VI.  The  Design  and  Parametric  Studies  of  Laminar  Flow 
Bodies 

This  proceedings  volume  includes  the  meeting  program,  list  of  attendees, 
and  abstracts  of  the  papers. 
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Program 

4 

LOW-SPEED  BOUNDARY-LAYER  TRANSITION  WORKSHOP:  II  « 


Monday,  September  13,  1976 


8:30  a.m.  WELCOME:  Dr.  E.  C.  Gritton,  Head,  Physical  Sciences  Department, 

The  Rand  Corporation 

Laminar- Flow  Research 

Dr.  G.  L.  Donohue,  Defense  Advanced  Research  Projects  Agency 
9:00  a.m.  KEYNOTE  SPEAKER:  Dr.  E.  Reshotko,  Case  Western  Reserve  University 

Morning  Section:  I.  NUMERICAL  INVESTIGATIONS  OF  NAVIER-STOKES  EQUATIONS 
Chairman:  W.  S.  King,  The  Rand  Corporation 

1-1  Numerical  Investigation  of  Nonlinear  Wave  Interaction  in  a Two-Dimen- 
sional Boundary  Layer , J.  W.  Murdock  and  T.  D.  Taylor,  The  Aero- 
space Corporation 

1-2  Numerical  Simulation  of  Transition  Phenomena  in  Incompressible  Flows, 

H.  Fasel,  Institut  A fur  Mechanik,  Universitat  Stuttgart 

1-3  Fundamental  Dynamical  Aspects  of  Indeterminacy  and  Incompleteness  of 
the  Navier-Stoke8  Equations  Bearing  on  Transition,  P.  Lieber, 

K.  S.  Wen.  and  J.  S.  Marks,  University  of  California,  Berkeley 

Afternoon  Section:  II.  TRANSITION  AND  STABILITY  THEORY 
Chairman:  L.  M.  Mack,  Jet  Propulsion  Laboratory 

II-l  Numerical  Simulation  of  Transition,  S.  A.  Orszag,  Massachusetts 
Institute  of  Technology 

II-2  The  Effects  of  Acceleration  on  the  Growth  Rate  of  Small  Disturbances 
in  a Laminar  Boundary  Layer,  C.  von  Kerczek,  Naval  Ship  Research 
and  Development  Center 

II-3  On  Boundary -Layer  Transition  Studies,  S.  I.  Cheng,  Princeton  University 

II-4  Transition , Pressure  Gradient,  Suction,  Separation,  and  Stability 
Theory,  J.  L.  van  Ingen,  Delft  University  of  Technology, 

The  Netherlands 

I 1-5  Nonparallel  Stability  of  Boundary  Layers,  W.  S.  Saric  and  A.  H.  Nayfeh, 
Virginia  Polytechnic  Institute  and  State  University 

II-6  Wave  Mechanical  Aspects  of  Transition , M.  T.  Landahl,  Massachusetts 
Institute  of  Technology 
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Tueeday,  September  14,  1976 


Morning  Section:  III.  PREDICTIONS  OF  TRANSITION  BY  MODEL  EQUATIONS 
Chairman:  M.  V.  Morkovin,  Illinois  Institute  of  Technology 

III-l  Theoretical  and  Experimental  Transition  Studies  at  ONERA/CERT, 

R.  Michel,  ONERA-CERT-DERAT 

III-2  Turbulence-Model  Predicted  Effects  of  Surface  Heat  Transfer  on 

Boundary -Layer  Transition , D.  C.  Wilcox  and  T.  L.  Chambers,  DCW 
Industries,  Inc. 

II 1-3  On  the  Application  of  Second-Order  Closure  Models  to  Boundary -Layer 
Transition , M.  L.  Finson,  Physical  Sciences,  Inc. 


IV.  THE  EFFECTS  OF  FREESTREAM  DISTURBANCES  AND  SURFACE  ROUGHNESS  ON  TRANSITION 
Chairman:  M.  V.  Morkovin,  Illinois  Institute  of  Technology 

IV- 1 An  Analytical  Study  of  the  Effect  of  Surface  Roughness  on  the  Sta- 
bility of  a Heated  Water  Boundary  Layer , D.  R.  S.  Ko  and  M.  Kosecoff, 
Physical  Dynamics,  Inc. 

IV-2  On  the  Effect  of  Freestream  Turbulence  on  Boundary -Layer  Transition, 

L.  M.  Mack,  Jet  Propulsion  Laboratory 
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IV- 3 The  Effects  of  Finite-Amplitude  Disturbances  on  the  Mean  Flow  Pro- 
files in  a Laminar  Boundary  Layer,  C.  L.  Merkle,  D.  R.  S.  Ko,  and 
T.  Kubota,  Flow  Research,  Inc. 

IV-4  Fluctuations  in  a Boundary  Layer  Introduced  by  Traveling -Wave  Irro- 
tational  Freestream  Disturbances , H.  L.  Rogler,  Case  Western  Re- 
serve University 


Afternoon  Section:  V.  EXPERIMENTAL  INVESTIGATIONS  OF  TRANSITION 
Chairman:  A.  M.  0.  Smith,  University  of  California,  Los  Angeles 

V-l  Stability  of  a Heated  Water  Boundary  Layer  for  Non-Uniform  Wall 
Temperature  Distributions , A.  Strazisar  and  E.  Reshotko,  Case 
Western  Reserve  University 

V-2  Transition  in  an  Axisyrmetric  Boundary  Layer  with  Wall  Heating, 

S.  J.  Barker,  Poseidon  Research  and  University  of  California, 

Los  Angeles 

V-3  Pressure  Fluctuations  in  the  Transition  Regions  of  Forebodies  of  Revo- 
lution, T.  T.  Huang,  Naval  Ship  Research  and  Development  Center 

V-4  A Comparison  between  Measured  and  Computed  Locations  of  Transition  on 
Nine  Forebodies  of  Revolution,  J.  L.  Power,  Naval  Ship  Research 
and  Development  Center 
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Wednesday,  September  15,  1976 

VI.  THE  DESIGN  AND  PARAMETRIC  STUDIES  OF  LAMINAR  FLOW  BODIES 
Chairman:  J.  Aroesty,  The  Rand  Corporation 

VI-1  The  Combined  Effects  of  Pressure  Gradient  and  Heating  on  Boundary- 
Layer  Stability  and  Transition,  A.  R.  Wazzan  and  2.  Gazley,  Jr., 

The  Rand  Corporation 

VI-2  Laminarization  of  Water  Boundary  Layers  over  Bodies  of  Revolution, 

A.  M.  0.  Smith  and  A.  R.  Wazzan,  University  of  California,  Los  Angeles 

VI-3  A Study  of  Heat  Addition  Effects  on  Laminar , Water  Boundary  Layers 

of  Axisymmetric  Bodies,  J.  J.  Eisenhuth,  Pennsylvania  State  University 

VI-4  Design  of  a Laminar , Low-Drag  Research  Vehicle,  R.  S.  Scotti, 

B.  H.  Carmichael,  and  D.  A.  King,  Rockwell  International 

VI-5  Axisymmetric  Vehicle  System  Design  for  Total  Propulsion  Power  Optimi- 
zation, F.  R.  Goldschmied,  Westinghouse  Research  Laboratories 


INFORMAL  PRESENTATIONS 

"Simple  Relations  for  the  Stability  of  Heated  Laminar  Boundary  Layers  in 
Water:  Modified  Dunn-Lin  Method,"  J.  Aroesty,  et  al..  The  Rand  Corporation. 

"Primary  and  Secondary  Instabilities  in  Unstably  Stratified  Shear  Flows," 
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NUMERICAL  INVESTIGATION  OF  NONLINEAR  WAVE  INTERACTION 
IN  A TWO-DIMENSIONAL  BOUNDARY  LAYER 
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1-1  Numerical  Investigation  of  Nonlinear  Wave  Interaction 
in  a Two-Dimensional  Boundary  Layer 

John  W.  Murdock  and  Thomas  D.  Taylor 
The  Aerospace  Corporation 
P.  O.  Box  92957 
Los  Angeles,  California  90009 

The  present  work  is  the  initial  phase  of  a study  whose  purpose  is  to 
develop  a computer  code  to  solve  the  Navier-Stokes  equations  (or  an  appropriate 
simplified  version  of  same)  in  the  Reynolds  number  range  of  transition  and  to 
study  the  various  stages  of  boundary  layer  transition.  The  unsteady,  two-dimensional, 
Navier-Stokes  equations  as  well  as  parabolized  vorticity  equations  have  been 
solved  for  flow  over  a flat  plate.  (The  parabolized  vorticity  equations  were  obtained 
from  the  Navier-Stokes  equations  by  neglecting  the  fourth  order  streamwise 
derivative  of  the  stream  function.)  The  physical  problem  considered  is  equivalent 
to  that  frequently  used  in  wind  tunnel  investigations  of  linear  stability  and  transition 
in  which  a two-dimensional,  time -periodic , disturbance  is  introduced  with  a 
vibrating  ribbon  in  the  boundary  layer.  This  situation  is  achieved  numerically  by 
imposing  an  upstream  boundary  condition,  at  a fixed  location  on  a plate,  which 
consists  of  the  sum  of  the  Blasius  solution  and  a time-periodic  solution  to  the  Orr- 
Sommerfeld  equation. 

The  dependent  variables  in  the  study  were  the  stream  function  and  vorticity; 
the  independent  variables  were  orthogonal  parabolic  coordinates.  The  equations 
w ere  solved  with  a so-called  spectral  scheme  which  was  selected  because  it  can 
be  more  accurate  for  comparable  machine  storage  than  a finite  difference  method. 

The  spectral  scheme  was  implemented  by  expanding  the  dependent  variables  in 
CJicbyshev  polynomials  in  both  space  dimensions.  The  solution  was  developed  by 
substituting  the  expansion  into  the  flow  equations  and  equating  the  coefficients  of 
equal  order  polynomials.  Since  the  coefficients  are  a function  of  time  the  result 
is  a set  of  coupled  ordinary  differential  equations.  These  equations  have  been 
solved  by  both  implicit  and  explicit  integration  schemes  for  the  complete  equations 
and  the  parabolized  equations. 

There  are  four  main  results  or  conclusions  which  have  been  achieved 
to  date.  The  first  is  that  it  is  possible  to  solve  the  two-dimensional  Navier- 
Stokes  equations  with  non-trivial  boundary  conditions  using  an  orthogonal  function 
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expansion  in  both  space  dimensions.  (A  c nnpletc  description  of  the  numerics 
will  be  contained  in  a forthcoming  paper.  ) 

The  second  conclusion  reached  in  this  study  is  that  for  the  type  of 
flows  considered  the  parabolized  vorticity  equations  are  an  adequate  model  of 
the  flow.  This  conclusion  is  substantiated  by  the  fact  that  these  equations  were 
used  to  derive  a modified  Orr-Sommerfeld  equation,  the  solution  of  which  agreed 
with  the  solution  of  the  usual  Orr-Sommerfeld  equation  to  4 significant  figures  for 
a typical  test  case.  This  is  not  surprising  if  the  order  of  magnitude  of  the  various 
terms  in  the  Orr-Sommerfeld  equation  is  considered.  (See  for  example  Gaster  (1)  .) 
Figure  1 shows  a comparison  of  the  Navier-Stokes  calculations  of  Fasel  (2)  with  the 
present  parabolized  vorticity  computations.  The  agreement  is  excellent;  similar 
agreement  has  been  found  between  the  present  Navier-Stokes  and  parabolized 
vorticity  solutions.  This  result  is  important  since  it  has  been  possible  to  generate 
solutions  to  the  parabolized  vorticity  equations  in  as  much  as  one -twentieth  the  time 
required  to  solve  the  same  physical  problem  with  the  Navier-Stokes  equations. 

The  saving  arises  because  the  additional  boundary  condition  which  must  be  imposed  on 
the  Navier-Stokes  equations  in  general  creates  a tliin  boundary  region  at  the  down- 
stream boundary  of  the  computation.  For  numerical  stability  this  boundary  region 
must  be  resolved;  this  requires  a fine  spatial  resolution  and  an  associated  decrease 
in  the  time  step  which,  in  general,  substantially  increases  running  time. 

The  third  result  of  this  work  is  that  the  nonlinear  effects  are  initially 
destabilizing  relative  to  the  linear  effects.  That  is  to  say  that  the  boundary 
layer  is  more  unstable  to  large  disturbances  than  to  small  ones.  This  is  an 
important  result  since  previous  analytic  nonlinear  stability  analyses  apply  only 
to  parallel  flow  and,  in  addition,  are  series  expansions  about  the  critical  point. 

The  present  results  are  for  a non-parallel  boundary  layer  and  may  be  obtained 
for  arbitrary  Reynolds  number  and  disturbance  wave  form. 

The  nonlinear  staDility  effects  have  been  investigated  by  inputing  Orr- 
Sommerfeld  disturbances  of  an  amplxt  ~.e  of  about  five  percent  of  the  mean  flow 
at  = 10^  and  solving  for  the  steady  state  disturbance  downstream.  In  Figure 
2 the  computed  linear  (A  = 0.  001)  and  non  linear  (A  = 0.  08)  wave  forms  arc  compared. 
Although  the  terms  "linear"  and  "non  linear"  are  used  to  distinguish  between  the 
two  computations  described  herein,  both  results  are  solutions  of  the  non  linear 
parabolized  vorticity  equations.  Figure  2 shows  significant  distortion  from  a 
modulated  sinusoidial  wave  form  by  nonlinear  effects  as  the  wave  propagates  from 
R = 105  to  2.  2 x 105. 

X 
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Inspection  of  Figure  2 demonstntes  that  it  is  difficult  to  physically  inter- 
pret the  nonlinear  effect  by  simply  inspecting  the  variation  of  velocity  as  a 

function  of  R . (Note,  however,  that  the  broadened  positive  peaks  and  large 

x 

narrow  negative  peaks  are  qualitatively  similar  to  the  hot  wire  measurements 
of  Klebanoff  and  Tidstrom  (3)  taken  in  a boundary  layer  excited  with  a large 
amplitude  Tollmien-Schlicting  wave.  ) Because  the  boundary  layer  is  being 
excited  at  a fixed  dimensionless  frequency,  u)  = 13.19,  it  may  be  Fourier  transformed 
in  time  once  all  transient  effects  are  washed  downstream.  The  Fourier  amplitude 
of  the  perturbation  velocity  is  illustrated  in  Figure  3 at  the  same  station  in  the 
boundary  layer  as  Figure  2.  Figures  2 and  3 are  at  a station  well  down  in  the 
shear  layer  ( T]  = 0.  2);  the  Fourier  amplitude  at  a station  near  the  maximum  of  the 
input  Orr-Sommerfeld  solution  is  shown  in  Figure  4.  The  nonlinear  curves 
(A  = 0.  08)  in  Figures  3 and  4 show  that  up  to  = 1.  3 x 10^  the  primary  mode 
is  changed  only  slightly  from  the  linear  behavior  while  the  second  mode  grows 
dramatically.  Thus,  in  this  rarg  e the  nonlinear  effects  are  destabilizing.  Beyond 
this  point  the  relative  stability  becomes  to  some  extent  a matter  of  definition. 

The  magnitude  of  the  mean  perturbation  from  the  Blasius  is  also  shown 
at  selected  points  in  Figures  3 and  4.  It  is  clearly  less  important  than  the 
second  harmonic.  Higher  harmonics  are  smaller  than  the  mean  and  are  therefore 
not  shown. 

Figures  3 and  4 show  the  Fourier  amplitude  of  the  perturbation  as  a 

function  of  R ; it  is  also  of  interest  to  consider  the  variation  across  the  boundary 
x 1 

layer  Asa  point  of  reference  consider  first  the  linear  behavior.  Figure  5 

" 5 

shows  the  disturbance  amplitude  as  input  at  R = 10  ; the  corresponding  curve 
5 x 

at  R^  = 2.  2 x 10  shows  the  peak  is  shifted  in  T)-space  as  the  wave  propagates 
downstream  and  is  amplified.  The  Fourier  amplitudes  of  the  various  nonlinear 

5 

modes  at  R^  = 1.  3 x 10  are  shown  in  Figure  6.  The  primary  mode  shape  is 

changed  only  slightly  from  the  linear  one.  The  general  shape  of  the  secondary  is 

similar  to  the  primary  but  the  peak  and  phase  reversal  point  are  closer  to  the 

wall.  The  amplitude  of  the  mean  is  positive  near  the  wall  with  a smaller  negative 

region  adjacent.  (Only  absolute  values  arc  shown  in  > e figure.)  Figure  7 shows 

5 

a similar  result  at  R^  = 2.  2 x 10  . In  this  case  the  normalized  amplitude  of  the 
nonlinear  primary  is  changed  somewhat  from  the  linear  curve.  The  mean 
contribution  again  has  a positive  region  near  the  wall  followed  by  a now 
significant  negative  region. 

Up  to  this  point  the  results  have  emphasized  the  Fourier  amplitude;  it 
is  also  of  interest  to  look  at  the  phase  relationship  between  the  primary  and 

. 1 
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secondary  disturbance.  The  usual  analytic  studies  of  nonlinear  stability  theory 

assume  that  the  secaidary  is  a harmonic  in  space  as  well  as  in  time.  Figure  8 

shows  a typical  phase  relationship  between  the  two  disturbances  (where  for  convenience 

both  amplitudes  are  normalized  to  unity).  The  results  show  that  the  secondary  is 

very  nearly  a spatial  harmonic  of  the  primary.  However,  a particularly  interesting 

feature  of  this  curve  is  the  rather  abrupt  change  in  the  relative  phase  at  about 
5 

R = 1.  7 x 10  which  is  about  the  location  of  the  change  in  shape  of  the  secondary 
x 

disturbance.  Although  this  point  needs  further  study,  a tentative  finding  of  this 
work  is  that  small  changes  in  the  relative  phase  appear  to  be  associated  with 
significant  changes  in  the  growth/decay  behavior  of  the  secondary  wave. 
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Fig.  5.  Change  in  Linear  Profile  with  R 
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1-2  NUMERICAL  SIMULATION  OP  TRANSITION  PHENOMENA 
IN  INCOMPRESSIBLE  FLOWS 

H.  Fasel 

Institut  A fUr  Mechanik 
UniversitSt  Stuttgart 
Stuttgart,  Germany 


For  the  investigation  of  stability  and  phenomena  of  initial 
transition  in  incompressible,  two-dimensional  flows  an  im-  ' 

plicit  finite  difference  method  for  the  solution  of  the 
complete  Navier-Stokes  equations  was  developed  (1).  The 

applicability  of  the  numerical  method  for  such  investiga-  i 

tions  was  demonstrated  by  simulations  of  Tollmien-Schlichting 
waves  in  boundary  layer  flows  on  a flat  plate.  In  the  numeri- 
cal model  these  waves  were  produced  by  introducing  periodic 
disturbances  at  the  upstream  boundary  A-D  of  a rectangular 
integration  domain  (Figure  1).  The  reaction  of  the  flow  to 
such  disturbances  was  then  directly  determined  by  solving  the 
Navier-Stokes  equations  within  the  rectangular  domain  using 
appropriate  boundary  conditions. 

For  perturbations  of  very  small  amplitudes  (maximum  disturb- 
ance amplitude  of  stream-wise  velocity  component  0.05  per- 
cent of  free  stream  velocity)  extensive  comparative  calcula- 
tions showed  good  agreement  with  linear  stability  theory  and 
experimental  measurements  (2,3).  Amplification  of  the  dis- 
turbances occurred  approximately  in  the  same  region  of 
Reynolds  number  and  disturbance  frequency  as  in  the  experi- 
ments and  in  linear  stability  theory,  particularly  when  non- 
parallel  flow  effects  were  taken  into  account  in  the  linear 
stability  theory  analysis  (4).  A comparison  of  the  numeri- 
cally determined  amplification  rates  and  amplitude  distribu- 
tions of  the  disturbances  with  linear  stability  theory  and 
experimental  measurements  was  equally  satisfactory. 
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The  usefulness  of  the  numerical  model  for  stability  and 
transition  studies  has  been  further  tested  by  applying 
it  to  another  type  of  flow,  i.e.  the  plane  Poiseuille 
flow  between  two  parallel  plates.  It  was  possible  to 
demonstrate  that  for  small  sinusoidal  disturbances,  which 
were  introduced  at  a fixed  downstream  location,  space-wise 
(in  downstream  direction)  amplification  of  the  disturbances 
can  occur  within  a certain  Reynolds  number  - frequency 
region.  In  Figure  2,  for  example,  amplification  curves 
for  two  typical  test  cases  are  shown,  one  for  a stable 
and  the  other  for  an  unstable  disturbance.  This  evidence 
of  spatially  growing  disturbances  in  plane  Poiseuille  flow 
was  also  verified  experimentally  by  Nishioka  et  al.  (5)* 

In  the  development  of  the  numerical  method  it  was  possible 
to  avoid  any  restrictions  with  regard  to  amplitude  or  form 
of  the  perturbations.  Thus,  the  numerical  model  lends  it- 
self for  the  investigation  of  a number  of  other  aspects 
of  stability  and  transition  for  which  treatment  with  other 
theoretical  models  becomes  extremely  difficult.  One  of  the 
main  points  of  interest  in  this  field  is  focused  on  the 
investigation  of  later  stages  of  the  transition  process 
which  follow  the  amplification  of  unstable  Tollmien-Schlich- 
ting  waves. 

For  an  investigation  of  the  two-dimensional  non-linear 
development  in  the  transition  process,  numerical  experiments 
are  performed  in  which  the  boundary  layer  flow  is  disturbed 
at  the  upstream  boundary  A-D  of  the  rectangular  integra- 
tion domain  (Figure  1)  with  periodic  perturbations  of  larger 
amplitudes.  In  these  experiments,  disturbance  amplitudes  of 
the  stream-wise  velocity  component  are  increased  up  to  5 per- 
cent of  the  free  stream  velocity.  In  Figure  3 typical  results 
are  shown  for  calculations  with  a Reynolds  number  (based  on 
displacement  thickness)  of  635  and  a frequency  parameter 

(F  = 6v/U*)  of  1.32* 10*  . In  these  calculations  the  loca- 
00 


-31- 


tion  of  the  upstream  boundary  corresponds  to  a point  on  the 
lower  branch  of  the  neutral  curve  of  linear  stability  theory. 

The  amplification  curves  for  different  disturbance  amplitudes 
which  are  plotted  in  Figure  3 clearly  indicate  a dependency 
of  the  amplification  rate  on  the  amplitude  of  the  disturbances. 

As  expected,  the  influence  of  the  magnitude  of  the  disturbance 
amplitude  for  amplitudes  below  one  percent  is  practically 
negligible.  However,  for  larger  amplitudes  an  increase  in 
disturbance  amplitude  has  considerable  influence  on  the  ampli- 
fication rate  of  the  disturbances. 

In  these  investigations  strong  emphasis  is  also  placed  on  a 
thorough  study  of  the  role  of  other  harmonic  disturbance  wave  ' 

components  in  the  transition  process.  The  numerical  results 
indicate  the  existence  of  other  harmonic  wave  components  super- 
imposed on  the  primary  disturbance  wave.  The  influence  of  these 
wave  components  seems  to  become  more  and  more  dominant  when 

< 

disturbance  amplitudes  are  increased. 

Another  effort  involves  a detailed  study  of  the  influence 
of  a single  roughness  element  on  the  transition  process  in 
a flat  plate  boundary  layer  flow,  corresponding  to  the  exper- 
imental work  of  Klebanoff  and  Tidstrom  (6).  Because  of  the 
simpler  geometry,  however,  a backward  facing  step  was  chosen 
as  roughness  element,  as  shown  in  Figure  4.  In  this  study  the 
step  heights  are  ranging  from  one  half  to  two  displacement 
thicknesses.  By  means  of  comparative  calculations  with  iden- 
tical flow  conditions  for  a flat  plate  with  the  step  and  with- 
out the  step  the  effects  of  the  roughness  element  on  the  tran- 
sition processes  can  be  clearly  observed.  For  example,  when 
random  disturbances  are  introduced  at  a fixed  downstream 
station,  a certain  selection  process  takes  place  for  both 
cases  where  disturbances, that  are,  according  to  linear  stabil- 
ity theory,  most  unstable  for  the  particular  Reynolds  numbers, 
are  singled  out.  However,  the  selection  process  is  much  more 
pronounced  when  the  roughness  element  is  present. 
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Figure  2.  Amplification  curves  for  spatially  growing 
disturbances  in  plane  Poiseuille  flow 
(Disturbance  amplitude  for  maximum  of  stream- 
wise  velocity  component  0.05%  of  center- 
line  velocity  Uc  ). 


Figure  3*  Amplification  curves  for  different  disturbance 
amplitudes  (maximum  of  stream-wise  velocity 
component  in  relation  to  free  stream  velocity). 


Figure  4.  Integration  domain  for  flow  over  a backward 
facing  step  (single  roughness  element). 
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Figure  1.  Integration  domain  for  flat  plate 
calculations 
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1-3  FUNDAMENTAL  DYNAMICAL  ASPECTS  OF  INDETERMINACY  AND  INCOMPLETENESS 
OF  THE  NAVIER-STOKES  EQUATIONS  GEARING  ON  TRANSITION 

by 

Paul  Lieber,  K.  S.  Wen  and  J.  S.  Marks 

The  application  of  the  Navier-Stokes  equations  to  the  calculation  of  actual 
flows  and  in  particular  transition,  by  analytical  and  numerical  procedures,  is 
predicated  on  the  supposition  that  they  are  in  fact  complete  in  specifying  the 
dynamical  constraints  necessary  for  so  obtaining  this  information.  We  demonstrate 
on  the  basis  of  the  first  principles  of  mechanics  which  circumscribe  and  thus 
limit  the  dynamical  information  content  of  the  Navier-Stokes  equations,  that  this 
basic  presupposition  must  be  fundamentally  questioned.  This  is  done  by  demonstrat- 
ing the  primacy  and  autonomy  of  force  relative  to  the  principles  of  mechanics  and 
that  this  important  fact  carries  over  into  the  Navier-Stokes  theory  as  well  and 
accounts  for  the  most  fundamental  dynamical  aspect  of  indeterminacy  and  incom- 
pleteness of  this  theory.  This  fundamental  aspect  of  indeterminacy  which  is 
predicated  on  the  primacy  of  the  forces,  allows  and  calls  for  the  specifications 
of  additional  dynamical  contraints  which  augment  and  complete  the  Navier-Stokes 
theory  for  the  purpose  of  calculating  actual  realizable  flows.  We  present  in 
outline  analytical  and  numerical  results  obtained  by  the  application  of  hydro- 
namical  variational  principles  that  were  conceived  and  designed  for  this  purpose. 

Certain  experimental  results  which  bear  both  directly  and  indirectly  on 
transition  and  which  experimentally  display  other  significant  aspects  of  incom- 
pleteness and  indeterminacy  of  the  Navier-Stokes  equations  are  described  as  well 
as  the  experimental  arrangements  that  produced  them.  The  theoretical  consider- 
ations that  led  to  the  conception  of  these  unconventional  hydrodynamical  experi- 
ments and  to  the  expectations  that  led  to  the  identification  of  the  results 
obtained  are  outlined. 
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II-l  Numerical  Simulation  of  Transition 

Steven  A.  Orszag 

Department  of  Mathematics 

M.I.T. 

Cambridge,  Mass.  02139 


In  this  paper  we  shall  summarize  and  compare  results 
obtained  from  two  distinct  theoretical  approaches  to  the 
transition  problem.  In  the  first,  the  three-dimensional 
time-dependent  Navier-Stokes  equations  are  solved  in  a plane 
channel  with  appropriate  boundary  conditions.  The  evolution 
of  three-dimensional  disturbances  is  studied  with  regard  to 
cross-stream  wavelength  selection,  spanwise  and  boundary-layer- 
wise  flow  variations,  etc.  The  effects  of  surface  curvature 
will  also  be  discussed.  In  the  second  study,  a highly 
simplified  set  of  Galerkin  equations  was  used  to  predict 
transition.  Comparisons  with  the  direct  numerical  simulations 
will  be  reported. 
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THE  EFFECTS  OF  ACCELERATION  ON  THE  GROWTH 
DISTURBANCES  IN  A LAMINAR  BOUNDARY 
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II-2  THE  EFFECTS  OF  ACCELERATION  ON  THE  GROWTH  RATE  OF  SMALL 
DISTURBANCES  IN  A LAMINAR  BOUNDARY  LAYER 

by 

C.  von  Kerczek 

David  W.  Taylor  Naval  Ship  Research  and  Development  Center 
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The  effect  of  small  acceleration  on  the  amplification  of  disturbances  in  a 

laminar  boundary  layer,  as  might  occur  in  the  latter  stages  of  the  flight  of  a 

buoyancy  propelled  body,  is  considered. 

Consider,  for  example,  a circular  cylinder  whose  radius  is  15  feet  and 

terminal  velocity  is  40  knots  (i.e.,  Reynolds  number  = 67  x 10^).  We  suppose 

-0  1 

that  it  accelerates  according  to  the  following  law  UjC'C')  =(l  - e ) where 

is  the  dimensionless  velocity  of  the  cylinder  at  an  instant  of  dimension- 
less time  U . The  velocity  scale  used  is  the  terminal  velocity  40  knots  and  the 
time  scale  is  the  radius  of  the  body  divided  by  the  terminal  velocity.  The  above 
acceleration  is  a guess  of  what  might  be  required  in  order  to  be  near  the  terminal 
velocity  after  a distance  of  about  1000  feet. 

The  instantaneous  inviscid  velocity  distribution  on  the  surface  of  the 

-0  IT 

cylinder  is  then  U(x,T)  ■ 2(1  - e ’ ) sin  x where  x is  the  arclength  along 

the  circumference  of  the  cylinder  measured  from  the  front  stagnation  point.  If 
one  considers  the  temporal  growth  of  small  disturbances  in  the  boundary  layer, 
the  analysis  of  [1]  justifies  the  use  of  quasi-steady  stability  theory  to  deter- 
mine growth  rates.  We  assume  that  validity  of  the  quasi-steady  assumption  to 
compute  temporal  growth  rates  allows  the  same  approximation  for  computation  of 
spatial  growth  rates.  Thus,  we  assume  that  the  amplitude  A of  a disturbance 


in  this  accelerating  boundary  layer  is  given  by 

X.  = 

where  -o<^  is  the  instantaneous  growth  rate,  R^^  is  the  instantaneous  length 
Reynolds  number,  RjV  is  t*le  instantaneous  displacement  thickness  Reynolds 
number  and  x^  is  the  instantaneous  point  of  neutral  stability.  The  amplification 
(1)  must  be  calculated  following  the  disturbance  with  its  group  velocity,  see  [2]. 


I 


However,  this  second  requirement  may  be  relaxed  in  this  case  because  the  group 
velocity  scales  on  U (e.g.,  calculations  typically  show  group  velocity  near 
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0.2  U)  and  hence,  the  disturbances  traverse  the  girth  of  the  cylinder  before  the 
boundary  layer  changes  appreciably. 

As  a sample  calculation,  Figure  1 shows  the  results  of  the  unsteady  boundary 
layer  calculation  (by  the  time-dependent  Pohlhausen  momentum  integral  method) 
on  a sector  ahead  of  but  near  the  top  (x  = 0.5)  of  the  cylinder.  The 
is  the  Pohlhausen  parameter  defined  in  unsteady  flow  as 
where  £ is  a boundary  layer  thickness  and  V is  th«  kinematic  viscosity.  The 
various  curves  give  A-  at  different  instants  of  time.  The  amplification  (1)  was 
computed  at  the  instant  X,  = 12  which  is  when  the  body  has  reached  about  70  per- 
cent of  its  terminal  velocity.  At  this  instant  the  disturbance  with  nondimen- 
sional  frequency  0.7  x 10  6 traverses  the  arc  (moving  with  its  group  velocity) 
from  x = 0.4  to  0.5  in  Ar<o.s.  In  this  increment  of  time,  the  boundary 
layer  has  little  perceptible  change  as  can  be  inferred  from  Figure  1.  This  gives 
an  a posteriori  justification  of  the  quasi-steady  assumption.  Figure  2 shows 
JLv\  A/A  ,as  a function  of  x at  tJ  = 12.  The  curve  marked  (a)  is  what  is  obtained 
in  a steady  flow  having  the  terminal  velocity.  The  curve  marked  (b)  is  what  is 
obtained  when  the  flow  is  steady  at  70  percent  of  terminal  velocity  (i.e.,  effect 
of  acceleration  of  the  flow  on  A.  is  neglected) . The  curve  marked  (c)  is  what  is 
obtained  when  boundary  layer  acceleration  is  taken  into  account.  The  conclus- 
ion is  that  the  acceleration  mainly  affects  boundary  layer  profile  shape  (i.e., A-  )• 
Even  though  this  effect  on  A.  is  fairly  small,  it  can  have  a large  effect  on  the 
amplification  in  the  region  of  the  body  where  A would  be  near  zero  in  a steady 
flow  (i.e.,  near  x ■ 0.5). 
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II-3  ON  BOUNDARY  LAYER  TRANSITION  STUDIES 
Sin-I  Cheng 
Princeton  University 

The  development  of  turbulent  spots  in  the  late  stages  of  transition  is  as 
important  as  spot  generation  in  transition  studies.  Laminar  instability  is 
important  but  not  all  encompassing  even  in  the  early  stages  of  burst  formation. 

A turbulent  boundary  layer  is  visualized  as  an  aggregate  of  distinct  and 
overlapping  turbulent  spots  convecting  downstream  in  their  asymptotic  states. 

The  turbulent  front  of  each  asymptotic  spot  develops  in  its  local  laminar  environ- 
ment independently  of  its  neighbors  until  overlaps  occur,  where  upon  nonadditive 
union  of  spots  and  annihilation  of  fronts  result.  Incipient  transition  is  reached 
when  sufficient  number  of  bursts  and  turbulent  spots  are  present  so  that  the  sum 
of  the  transit  times  of  all  these  spots  under  asymptotic  development  will  exceed 
unity  (intermittency  factor  = 1)  somewhere  downstream.  The  resulting  outgoing 
turbulent  front  (mean)  will  eventually  redch  its  equilibrium  position  about  half 
way  toward  the  outer  edge  of  the  turbulent  boundary  layer  under  the  large  con- 
vective effects  there.  The  ingoing  front  propagates  into  the  thin  wall  layer 
where  it  is  promptly  dissipated.  A fully  turbulent  boundary  layer  then  results. 

As  such,  the  mean  turbulent  properties  within  a developing  spot  are  those  in  the 
fully  turbulent  boundary  layer  further  downstream. 

Kinematics,  dynamic  and  energy  considerations  of  the  asymptotic  development 
of  a spot  and  its  interaction  with  the  "mean  flow"  permit  an  estimate  of  the 
asymptotic  states  and  of  the  transit  time  of  a spot.  Such  estimates  are  in  sub- 
stantial agreement  with  observations.  A criterion  of  incipient  transition  is 
then  formulated  in  terms  of  the  transit  times  of  all  the  asymptotic  spots  resulting 
from  the  population  of  bursts  and  spots  at  a given  station.  A burst  is  presumed 
to  be  generated  when  the  disturbance  amplitude  exceeds  what  is  needed  to  extract 
sufficient  energy  from  the  mean  flow  for  establishing  a self-sustaining  turbulent 
front  of  an  asymptotic  spot.  This  cut  off  amplitude  increases  with  the  scale  and 
frequency  of  the  disturbance  to  enhance  such  selectivity  in  burst  formation. 

If  external  disturbances  are  sufficiently  large  and  numerous,  the  proposed 
transition  criterion  provides  specific  functional  relations  between  transition 
Reynolds  number  and  the  external  disturbances  that  describe  adequately  the  obser- 
vational data.  If  insufficient,  transition  will  require  additional  spots  from 
growing  laminar  disturbances  with  sufficiently  large  amplification  rate.  The 
many  faces  of  transition  can  then  be  put  in  proper  perspective  to  further  better 
understanding  and  the  construction  of  some  practically  useful  transition  criterion. 
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II-4 

TRANSITION,  PRESSURE  GRADIENT,  SUCTION,  SEPARATION  AND 
STABILITY  THEORY. 

J.L.  van  Ingen 

Department  of  Aerospace  Engineering,  Delft  University  of  Technology, 
Kluyverweg  1,  Delft,  The  Netherlands. 


SUMMARY. 

A semi-empirical  method  is  presented  for  the  prediction  of  transition 
in  two-dimensional  incompressible  flows  with  pressure  gradient  and 
suction.  Included  is  the  case  of  the  laminar  separation  bubble,  where 
transition  is  preceded  by  laminar  separation. 

The  method  employs  linear  stability  theory  to  calculate  the  amplification 
factor  a for  unstable  disturbances  in  the  laminar  boundary  layer. 

(o  is  defined  as  the  natural  logarithm  of  the  ratio  between  the 
amplitude  ot  a disturbance  at  a given  instant  or  position  to  the 
amplitude  at  neutral  stability).  It  is  found  that  at  the  experimentally 
determined  transition  position  the  calculated  amplification  factor 
for  the  critical  disturbances  attains  nearly  the  same  value  (about 
10)  in  many  different  cases  for  flows  with  low  free  stream  turbulence 
levels.  An  attempt  is  made  to  include  the  effects  of  higher  free 
stream  turbulence  levels  by  allowing  the  critical  amplification  factor 
to  decrease  with  increasing  free  stream  turbulence. 


NOTATION. 


The  symbols  used  are  the  conventional  ones 
stability  theory.  To  avoid  confusion  a few 
specifically  below, 
c reference  length 


m 


0^  dU 
V dx 


for  boundary  layer  and 
of  them  are  mentioned 


U velocity  at  edge  of  boundary  layer 
reference  velocity 
U - U/U 

00 


I 


> 
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x or  s distance  along  contour  of  body 

x distance  along  chord 

c 

X 

X — 

c 

s s/c 

subscript  sep  refers  to  conditions  at  separation. 


r 
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LINEAR  STABILITY  THEORY. 


In  linear  stability  theory  a given  two-dimensional  laminar  main 
flow  is  subjected  to  sinusoidal  disturbances  with  a disturbance 
stream  function: 


'J' 


4>(y)  ei(arttt) 


(1) 


For  the  spatial  mode  U)  is  real  and  a is  complex  a = a + i a. . 

“Cl  • x 

This  leads  to  a factor  e i in  the  disturbance  amplitude  and  a 

follows  from:  tt 

x U^c  x 

o = / -a.dx  = 10  / T.U  dx 

Xq  1 V rr 

xo 


(2) 


where  xq  is  the  streamwise  position  where  the  disturbance  with 

frequency  u>  is  neutrally  stable. 

T is  defined  as:  -a  9 

T = — 106  (3) 

R0 

In  the  temporal  mode  the  same  expression  (2)  for  o is  found  with 
a different  definition  for  T. 

It  is  clear  that  0 is  a function  of  x and  to  for  a given  boundary 
layer;  a can  be  calculated  as  soon  as  stability  diagrams  are 
available  for  the  velocity  profiles  for  successive  streamwise 
positions  x. 

For  a long  time  Pretsch' stability  diagrams  for  the  temporal  stability 
of  the  Hartree  similar  velocity  profiles,  have  been  the  only  source 
of  detailed  stability  data  for  flows  with  non-zero  pressure  gradient [6] . 
Results  for  the  spatial  stability  of  the  Hartree  flows  have  been 
given  by  Wazzan,  Okamura  and  Smith  [^7j  and  Kummerer  [8]  , stability 
diagrams  for  the  reversed  flow  solutions  of  the  Falkner-Skan  equation 
fiave  been  obtained  by  Taghavi  and  WazzanfllJ. 
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STABILITY  AND  TRANSITION  OF  THE  FLAT  PLATE  BOUNDARY  LAYER. 

Fig.  la  shows  a for  the  flat  plate  according  to  Pretsch  for 
different  non-dimensional  frequencies  ^ . The  envelope  of  these 
curves  gives  the  maximum  value  of  a for  each  streamwise  position. 

In  what  follows  we  will  in  general  mean  this  maximum  value  when 
we  mention  a.  The  curve  labelled  3 in  fig.  lb  is  the  envelope 
according  to  [7J  and  {^8j ; the  curve  labelled  2 will  be  discussed 
later.  A well  known  result  for  the  experimentally  determined 
transition  region  is  due  to  Schubauer  and  Skramstad  [12J . They 
find  for  low  free  stream  turbulence  levels  Reynoldsnumbers  at 
beginning  and  end  of  transition  equal  to  2.8  x 10^  and  3.9  x 10^ 
respectively.  To  these  Reynolds  numbers  correspond  certain 
values  Oj  and  a ^ for  cr  which  are  indicated  in  table  1. 

FIRST  VERSION  OF  THE  PREDICTION  METHOD  (1956). 

The  present  author  used  Pretsch  charts  in  £l]  to  calculate 
amplification  factors  for  an  airfoil  section  (EC  1440)  at  different 
values  of  angle  of  attack  and  Reynolds  number. 

It  was  shown  that  Oj=7.6  and  CT2“9.7  gave  a reasonably  accurate 
prediction  of  the  transition  region.  Smith  and  Gamberoni  [3], 
defining  a transition  point  rather  than  a transition  region  found 
that  a»9  would  correlate  different  transition  experiments  reasonably 
well. 

Although  it  is  clear  that  a transition  criterion  should  be  based  on 
the  actual  amplitude  of  the  disturbance,  rather  than  on  an  amplification 
ratio,  the  method  has  been  used  extensively.  Its  success  may  have  been 
due  to  the  fact  that  the  initial  disturbances  - due  to  free  stream 
turbulence  for  instance  - have  been  about  the  same  for  the  cases 
investigated. 

Another  way  to  explain  the  success  of  the  method  may  be  that  a is 
a suitable  factor  in  which  different  factors,  known  to  influence 
transition,  may  be  correlated. 

SECOND  VERSION,  ALSO  APPLICABLE  TO  FLOWS  WITH  SUCTION. 

In  1965  the  present  author  extended  the  method  to  the  case  of  two- 
dimensional  incompressible  boundary  layers  with  suction  [2j. 

Since  at  that  time  the  Pretsch  charts  were  still  the  only  source 
of  detailed  information  on  amplification  rates,  some  drastic 
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simplifying  assumptions  had  to  be  made.  First  it  was  assumed  that 
all  possible  stability  diagrams,  including  those  for  suction 
boundary  layers,  formed  a one-parameter  family  with  the  critical 
Reynolds  number  as  parameter.  Furthermore,  it  was  assumed  that  the 
critical  Reynolds  number  could  be  determined  from  an  approximation 
formula  due  to  Lin.  The  suction  boundary  layer  was  calculated  using 
a two-parameter  method  of  integral  relations.  This  necessitated  a 
new  "calibration"  of  the  transition  prediction  method  against  the 
flat  plate  without  suction,  leading  to  curve  2 in  fig.  lb  with 
Oj=9.2  and  (^“l  1.2. 

To  facilitate  the  amplification  calculations  using  a computer 
Pretsch'  charts  have  been  brought  in  a tabular  form.  Fig.  2 shows 
an  application  to  the  EC  1440  airfoil;  some  results  for  an  airfoil 
with  suction  through  a porous  surface  are  shown  in  figs  3 and  4. 

In  view  of  the  many  simplifying  assumptions  which  had  to  be  made  the 
correspondence  between  theory  and  experiment  may  be  considered  to  be 
good. 

Since  1965  this  version  of  the  method  has  been  included  in  a computer 
program  for  the  analysis  and  design  of  airfoil  sections  [[13J.  The 
streamwise  position  for  the  end  of  the  transition  region  (determined 
by  c^)  has  been  used  as  the  starting  point  for  the  turbulent  boundary 
layer  calculation. 

It  has  been  found  that  an  improved  transition  prediction  could  be 
made  by  allowing  the  value  for  <?2  to  vary  from  11.2  for  favourable 
and  zero  pressure  gradient  to  about  20  for  boundary  layers  near 
separation.  (In  the  last  version  is  again  more  nearly  constant). 

In  general  the  position  of  transition  was  predicted  within  a few 
percent  of  the  chord.  An  example  of  application  of  this  airfoil 
analysis  program  taken  from  £l4^  is  shown  in  fig.  5 . The  airfoil 
investigated  is  that  of  the  horizontal  tailplane  of  the  Italian  sail- 
plane M300  "Aliante".  The  airfoil  was  designed  by  cambering  the 
NACA  63^-018  section.  The  tailplane  is  produced  through  an  extrusion 
process  which  caused  appreciable  surface  waviness.  An  actual  specimen 
of  this  tailplane  was  tested  as  a two-dimensional  model  in  the  low 
speed  wind  tunnel  of  the  Department  of  Aerospace  Engineering  at  Delft. 
It  was  found  that  the  surface  waves  caused  early  transition  in  a 
certain  angle  of  attack  range;  this  could  be  remedied  by  smoothing 
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the  forward  part  of  the  surface.  The  calculation,  starting  from 
the  airfoil  coordinates  for  both  conditions,  predicted  this  change 
quite  well. 

It  should  be  stressed  again  that  the  present  method  may  be  considered 
as  a method  to  correlate  different  transition  experiments.  The 
calculated  amplification  factors  need  not  have  a precise  physical 
meaning.  It  is  however  a definite  advantage  of  the  method  that  linear 
stability  theory  is  used  which  has  proved  to  be  a valuable  tool  to 
describe  the  early  phases  of  the  transition  process.  It  should  also 
be  observed  that  inaccuracies  in  one  of  the  elements  of  the  method 
(viz.  boundary  layer  calculation;  calculation  of  the  critical 
Reynolds  number  using  Lin's  formula;  the  stability  diagrams  used) 
may  have  been  neutralized  by  inaccuracies  in  another  element.  Hence 
if  any  element  is  changed,  a new  calibration  is  necessary. 

An  important  imperfection  in  the  second  version  of  the  method  was 
that  the  stability  characteristics  in  laminar  separation  bubbles 
were  obtained  by  extrapolation  from  the  attached  flow.  This  may 
have  been  the  cause  of  the  high  values  of  required  to  predict 
accurately  the  end  of  the  transition  region  in  boundary  layers  near 
to  or  after  separation. 

A SHORT  CUT  METHOD  TO  PREDICT  TRANSITION  IN  SEPARATION  BUBBLES. 

In  [5]  the  present  author  published  a short-cut  method  to  predict 
transition  in  separated  flow.  The  method  is  based  on  the  stability 
diagrams  for  reversed  flows  due  to  Taghavi  and  Wazzan  [llj  and 
some  additional  calculations  by  the  present  author  for  the  limiting 
stability  characteristics  when  Kg"***  using  the  inviscid  stability 
equation  (Rayleigh  equation).  The  following  assumptions  are  made 

1)  U,  6 and  R^  in  the  separation  bubble  are  independent  of  x and 
equal  to  their  values  at  separation.  Then  a constant  value  of 

(i)0 

0)  also  means  a constant  value  of  — jj. 

2)  The  separation  streamline  is  straight,  and  leaves  the  wall  at 
an  angle  y determined  by: 

■«<»  ■ 37-  <*> 

^ v'sep 

where  B is  a constant  equal  to  17.5. 
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3)  The  Reynolds  number  is  so  high  with  respect  to  the  (very  low) 
critical  Reynolds  number  that  the  stability  characteristics 
are  given  with  sufficient  accuracy  by  the  limiting  values 
determined  from  the  inviscid  stability  equation. 

(l)0 

Then  -0^9  only  depends  on  the  value  of  and  the  velocity 
profile  shape  parameter. 

Finally  we  introduce  the  shape  parameter  z«g  x msep*  where  g is  the 

height  of  the  separation  streamline  above  the  wall  devided  by  0 and 
2 

0 

mse  ' ~ v~  dx  3t  seParat^on*  Then  the  integration  w.r.t.  x in  (2) 
can  be  replaced  by  an  integration  w.r.t.  to  z leading  to: 


a 


(V..,  s 


B.m 


sep 


(-a^9)dz 


(5) 


(a  similar  result  may  be  obtained  for  a small  region  upstream  of 

To0  , 

separation  when  integration  w.r.t.  i,  - — ^ is  used). 

The  inviscid  instability  for  different  values  of  the  Hartree  parameter 

4 

8 is  shown  in  figs  6 and  7.  Values  of  10  / (-ou9)dz  are  shown  in 

fig.  8 for  different  values  of  ~ together  with  the  envelope  giving 
the  maximum  value  I of  the  integral  as  a function  of  Z.  (See  also 
table  2).  Hence  in  the  separation  bubble  we  have: 

I 


10"*  (Rq) 


a ■ 


B.m 


8 sep 


(6) 


sep 


Using  this  short-cut  method  it  was  found  in  [5]  that  a 2 for 
separation  bubbles  on  an  airfoil  in  a small  "noisy"  tunnel  was 
about  12.5  (fig.  9).  For  separation  bubbles  on  a circular  cylinder 
with  a tapered  tail  in  the  large  low  turbulence  wind  tunnel,  values 
of  a 2 between  13.2  and  15.7  were  found,  depending  on  the  wind  speed. 
Using  the  same  Short-cut  method  Van  der  Meulen  [15]  obtained 
for  a body  of  revolution  in  a small  high  speed  water  tunnel. 


PRESENT  STATUS  OF  THE  TRANSITION  PREDICTION  METHOD. 

All  stability  data  obtained  from  [7,8, 1 l]  and  the  inviscid  stability 
calculations  mentioned  in  the  preceding  section,  have  been  reduced 
to  a table  containing  about  300  numbers. 

Using  this  table,  the  amplification  rate  T can  easily  be  obtained  for 
any  velocity  profile,  as  soon  as  the  critical  Reynolds  number  is  known. 
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The  present  author  employs  a boundary  layer  calculation  method  ^5  J 
which  for  attached  flow  is  similar  to  Thwaites'  method.  It  contains 
an  extra  parameter  however,  which  makes  the  prediction  of  the 
separation  position  as  accurate  as  for  Stratford's  two-layer  method. 
In  separated  flows  an  integral  method  is  used  in  which  the  shape  of 
the  separation  streamline  is  prescrioed.  Both  for  attached  and 
separated  flow  the  primary  profile  shape  parameter  is  m/mse^.  The 
critical  Reynolds  number  is  a function  of  m/m^^;  this  function  is 
assumed  to  be  equal  to  that  obtained  for  the  Falkner-Skan  solutions. 
From  calculations  with  the  full  method  it  has  been  found  that  the 
short-cut  method,  described  in  the  preceding  section,  gives  a very 
good  approximation  in  separation  bubbles.  Furthermore  it  has  been 
found  that  the  values  of  Oj  and  c^,  when  transition  occurs  near 
separation  are  much  nearer  to  the  flat  plate  values  than  for  the 
second  version.  It  can  now  be  expected  that  Oj  and  will  be  more 
or  less  constant  for  flows  with  the  same  initial  disturbances. 
However,  and  a ^ may  have  to  vary  with  the  level  of  initial 
disturbances  due  to  free  stream  turbulence  and  noise. 

From  curve  3 in  fig.  1 and  table  1 it  follows  that  Oj  ■ 8.3  and 

0 2 “ 10.4  if  Schubauer  and  Skramstad's  transition  results  for  the 
flat  plate  are  used.  From  Spangler  and  Wells'  experiment  on  a flat 
plate  in  a tunnel  with  reduced  background  noise  [l6,17j  and  from 
the  authors  own  experiments  somewhat  larger  values  for  Oj  ( 12  ) 

and  (14.5  ) would  be  obtained,  jaffe,  Okamura  and  Smith  [9] 
applied  their  solution  technique  for  the  Orr-Sommerfeld  equation 
to  velocity  profiles  that  had  been  obtained  numerically  for  two- 
dimensional  and  axi-symnetric  flows.  They  find  Cj-8.3  for  the 
Schubauer  and  Skramstad  results  and  0j“11.8  for  Well's  results; 
for  a large  number  of  flows  with  pressure  gradient  0j  values  ranging 
from  6.8  to  12.1  were  obtained.  A good  overall  correlation  of 
transition  position  was  obtained  using  o^lO. 


RELATION  BETWEEN  a] , 02  AND  FREE  STREAM  TURBULENCE. 

Although  it  is  clear  that  the  initial  disturbances  cannot  be  sufficiently 
characterised  by  the  r.m.s.  value  of  free  stream  turbulence  alone,  it 
will  be  attempted  in  the  present  section  to  find  a relation  between 
Oj,  02  end  the  r.m.s.  free  stream  turbulence  Tu.(in  X) . 
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In  many  different  papers  relations  between  Tu,  Rg  or  Rx  at  transition 
have  been  given  for  the  flat  plate.  The  measured  transition  positions 
may  be  converted  to  a-values  using  curve  3 from  fig.  lb.  Then  a will 
decrease  when  Tu  increases;  fig.  10  shows  a collection  of  these  data; 
for  Tu  > 0.1%  the  relation  used  by  Mack  in  fig.  3 of  [l8j  can  be 
approximated  by: 

Oj  -2.13 -6. 1810  log  Tu  (7) 

while  for  a reasonable  approximation  is: 

Cj  ° 5 -6.1810  log  Tu  (8) 

For  values  of  Tu  < 0.1%  there  is  much  more  scatter  because  in  this 
region  sound  disturbances  may  become  the  factor  controlling  transition 
rather  than  turbulence.  We  may  also  use  the  relations  (7)  and  (8) 
for  Tu  < 0.1%;  but  then  we  should  define  an  "effective"  value  for  Tu. 

Of  course  this  does  not  solve  the  problem  because  we  can  only  define 
an  "effective  Tu"  for  a wind  tunnel  after  transition  experiments 
have  been  made  in  that  same  tunnel. 

At  the  time  of  writing  this  abstract  some  additional  measurements  in 
the  low  speed  low  turbulence  wind  tunnel  of  the  Department  of  Aerospace 
Engineering  are  being  evaluated.  Some  preliminary  results  show  that 
the  "effective  Tu"  even  may  increase  at  the  lower  windspeeds  where  Tu 
decreases.  It  is  thought  that  this  is  due  to  the  fact  that  the  critical 
frequencies  in  the  boundary  layer  may  better  be  matched  to  the  wind 
tunnel  noise  spectrum  at  lower  speeds. 

For  the  time  being  it  is  suggested  to  use  (7)  and  (8),  assuming  an 
effective  Tu  equal  to  0.1%  for  modern  wind  tunnels,  resulting  in 
Oj  ■ 8.3  and  * 11.2. 
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curve  no 
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°2 

stability  diagram  used 
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7.6 

D 

Fretsch,  flat  plate  (6=0);  version  1. 

2 

9.2 

11.2 

Pretsch,  stability  diagram  for  *^log(— ) . “2.345 

which  according  to  Lin’s  ormula  wouldv 
apply  to  the  flat  plate  velocity  profile  in 
version  2. 

3 

8.3 

10.4 

from  £ 7^J  and  [ 8j 

Table  1 ; Critical  values  for  o at  beginning  (Oj)  and  end  of  the 
transition  region  (o?)  on  a flat  plate  according  to  different 
stability  calculations.  Transition  Reynolds  numbers  2.8  and  3.9  x 10 
according  to  £l2j. 
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.360 
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.682 

483 
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4.249 

1331 

Table  2:  z and  I as  a function  of  the  Hartree  shape  parameter  6 
for  reversed  flows. 
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Fig.  10:  Relation  between  Oj,  0 ^ and  Tu  for  the  flat  plate. 
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n-5  Nonparallel  Stability  of  Boundary  Layers 

William  S.  Saric  and  Ali  Hasan  Nayfeh 
Department  of  Engineering  Science  and  Mechanics 
Virginia  Polytechnic  Institute  and  State  University 
Blacksburg,  Virginia  24061  U.S.A 

This  paper  investigates  the  stability  of  flows  that  are  not  parallel. 
That  is  to  say,  the  mean-flow  velocity  is  a function  of  two  dimensions 
and  thus,  it  has  a streamwise  as  well  as  a transverse  component.  In  this 
work,  we  confine  our  attention  to  the  following  class  of  flows:  boundary 
layers  with  variable  pressure  gradients  and  blowing  or  suction.  A common 
feature  of  these  flows  is  that  they  are  weakly  nonparallel  and  thus  many 
previous  stability  analyses  have  assumed  the  primary  flow  to  be  parallel 
that  is,  the  flow  is  one-dimensional  and  Identical  at  every  streamwise 
locati on. 

Recently,  considerable  attention  has  been  given  to  linear  stability 
theories  of  shear  flows  that  account  for  nonparallel  flow  effects  (Bouthier 
1973,  Nayfeh,  Saric  A Mook  1974,  Gaster  1975,  Saric  & Nayfeh  1975)  because 
classical  linear  stability  theories  which  treat  the  primary  flows  as  quasi- 
parallel flows  have  not  produced  satisfactory  results.  For  a flow  past  a 
flat  plate  (Blaslus  flow),  the  critical  Reynolds  number  predicted  by  parallel 
* stability  theories  is  about  30*  above  the  experimental  results  of  Schubauer  & 

Skramstad  (1947),  Ross,  Barnes,  Bums,  A Ross  (1970),  Wortman  (1955),  and 
Strazlsar,  Prahl , A Reshotko  (1975),  as  shown  in  figure  1.  A survey  of 
stability  analyses  and  their  application  to  transition  is  given  by  Reshotko 
(1976). 

For  a parallel  mean  flow,  the  spatially  varying  disturbance  stream- 
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function  can  be  expressed  as 

Y = c(y)exp{i  (k0x  - cut) } (1) 

where  c is  the  mode  shape,  u>  is  the  disturbance  frequency,  and  k0  is  a 
complex  constant.  Since  the  mode  shape  t is  assumed  to  be  independent  of 
x,  the  streamwise  behavior  of  any  disturbance  quantity,  such  as  the  velo- 
city, the  pressure,  and  the  kinetic  energy,  is  governed  by  the  exponent 
k0.  In  particular,  the  amplification  and  attenuation  rates  of  any  dis- 
turbance are  given  uniquely  by  the  imaginary  part  of  k0.  On  the  other 
hand,  for  a near  parallel  flow,  the  disturbance  stream  function  has  the 
form 

v = Aodyixjexpli  j (k0  + ek^dx  - iu>t}  (2) 

Thus,  the  effects  of  nonparallelism  are  to  make  k0  be  a function  of  the 
long  scale  x*  = ex,  where  e is  a small  parameter  characterizing  the  non- 
parallelism of  the  mean  flow,  to  produce  a correction  ek1(xJ)  to  k0,  and 
to  make  the  mode  shape  c vary  in  the  streamwise  direction.  Hence,  the 
streamwise  variation  of  each  flow  quantity  depends  on  its  distance  from 
the  wall.  Moreover,  at  each  distance  from  the  wall,  the  different  flow 
quantities  vary  differently  in  the  streamwise  direction.  Thus,  unless 
these  factors  are  taken  into  account,  meaningful  comparisons  cannot  be 
made  between  theoretical  and  experimental  results.  Figure  1 shows  good 
agreement  between  the  nonparallel  neutral  stability  curve  calculated  by 
Sarlc  & Nayfeh  (1975)  and  the  results  of  various  experiments. 

Recent  experiments  by  Strazlsar,  Prahl  & Reshotko  (1975)  in  a water 
tunnel  and  by  Kachanov,  Kozlov  & Levchenko  (1975)  In  a wind  tunnel  have 
provided  additional  data  which  can  be  used -for  comparison.  Figure  2 shows 
a comparison  of  the  maximum  amplification  rates  calculated  by  using  parallel 
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and  nonparallel  theories  and  those  determined  experimentally  by  Strazisar, 
Prahl  & Reshotko.  We  note  that  the  maximum  amplification  rates  calculated 
by  using  the  nonparallel  theory  is  approximately  l*j  times  those  calculated 
by  using  the  parallel  theory  for  low  Reynolds  numbers.  Figures  3 and  4 
show  good  agreement  between  the  amplitude  histories  calculated  by  using  the 
nonparallel  theory  and  those  determined  experimentally  by  Kachanov,  Koslov 
& Levchenko.  Thus  the  theory  is  in  good  agreement  with  available  experi- 
ments for  neutral  stability,  amplification  rate,  and  amplitude  history. 
Comparison  between  other  nonparallel  theories  are  made  in  the  paper. 

Our  nonparallel  stability  model  has  been  extended  to  include  the 
effects  of  pressure  gradients,  suction  or  blowing,  and  higher-order  boun- 
dary-layer terms. 


j 

i 
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Reynolds  Number 

Figure  1.  Comparison  of  experimental  results  and  nonparallel  theory 

for  Blasius  boundary  layer.  From  Strazisar,  Prahl,  & Reshotko 
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Section  II 

WAVE  MECHANICAL  ASPECTS  OF  TRANSITION 
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II-6  Wave  Mechanical  Aspects  of  Transition 

M.  T.  Landahl 
Mass.  Inst.  Technology 


Synopsis 

The  final  stage  of  transition  of  wall  bounded  shear  flows  is  usu- 
ally characterized  by  a sudden,  almost  explosive  onset  of  a short-wave 
secondary  instability  with  a characteristic  scale  much  smaller  than 
that  typical  of  the  unstable  Tollmien-Schlichting  waves.  Such  onset 
is  found,  for  example,  in  turbulent  spots  (see  e.g.  picture  by  Elder, 
J.F.M.  9^,  35),  and  in  the  formation  of  low-speed  "spikes,"  (Klebanoff, 
et  al.,  J.  Fluid  M.  12^,  1).  The  propagation  of  such  disturbances,  con- 
sidered as  wave  packets,  can  be  analyzed  by  aid  of  kinematical  wave 
theory  applied  to  propagation  in  a slowly  varying  nonhomogeneous  medium. 
This  theory  may  be  considered  as  an  asymptotic  one  to  first  order  in  a 
small  parameter  measuring  the  ratio  of  secondary  wave  length  to  the  char- 
acteristic scale  of  the  primary  inhomogeneity.  For  the  case  when  the 
large-scale  in  homogeneity  is  a travelling  wave  of  permanent  shape,  a 
critical  condition  for  strong  secondary  instability  was  found  to  arise 
(Landahl,  J.  Fluid  M.  5(5,  775)  when  the  group  velocity  of  the  secondary, 
locally  unstable  wave  is  equal  to  the  phase  velocity  of  the  large-scale 
(primary)  wave.  When  this  condition  is  met,  space-time  focusing  of  the 
small-scale  disturbance  on  the  large-scale  wave  will  take  place.  The 
simple  theory  was  found  to  give  results  in  good  agreement  with  measure- 
ments for  the  prediction  of  the  location  of  the  first  appearing  spike, 
as  well  as  for  the  frequency  of  the  secondary  instability  (Landahl,  loc. 
cit.).  Some  new  results  are  presented  from  breakdown  calculations  for 
a simple  stratified  shear  flow  model  (Fig.  1).  Breakdown  can  occur  for 
both  stable  and  unstable  primary  waves.  The  maximum  wave  slope  attain- 
able for  the  most  unstable  pure  shear  wave  (zero  Richardson  number)  be- 
fore breakdown  is  found  to  be  about  30°.  The  simple  breakdown  condition 
can  be  generalized  to  hold  for  arbitrary  time-dependent  three-dimensional 
primary  flows. 
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Its  implication  for  dispersive  and  nondispersive  systems  is  discus- 
sed and  a comparison  made  with  ordinary  ray  (space-)  focusing  (Fig.  2). 

New  developments  in  kinematical  wave  theory  for  general  classes 
nearly  conservative  fluid-flow  systems  (W.  Chin,  M.I.T.  Ph.D.  Thesis, 
1975)  show  that  shock-like  transition  giving  discontinuities  on  the 
scale  of  large-scale  flow  are  possible  for  systems  that  are  supercriti- 
cal in  the  sense  of  the  breakdown  condition.  The  discontinuities  arise 
both  in  the  large-scale  flow  and  in  the  amplitude  of  the  wave  train 
(Fig.  3).  Such  a transition  could  be  regarded  as  a generalization  of 
the  classical  hydraulic  pump  or  gas  dynamic  shock  wave,  and  is  proposed 
as  a unified  mechanism  for  fluid-flow  transition  phenomena. 


.1 


.2 


.3 


.4 


Largest  attainable  wave  slope  for  most 
a mpl  if  i eo'  prim  a ry  wave:  ( J0-  0 , k0--  0.4  ) 

k<A=0-62  ('30°  ) 

I’ifi.  1. 
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Section  III 

THEORETICAL  AND  EXPERIMENTAL  TRANSITION 
STUDIES  AT  ONERA  / CERT 
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IH-1  THEORETICAL  AND  EXPERIMENTAL  TRANSITION 
STUDIES  AT  ONERA  / CERT 

R.  MICHEL 

Departement  d'ASrothermodyanmique,  ONERA  / CERT,  Toulouse,  FRANCE 


In  the  beginning,  the  study  has  been  undertaken  essentially  with  the 
aim  of  defining  means  of  calculation  less  empirical  and  more  coherent  than  the 
available  criterions,  for  predicting  the  onset  and  development  of  transition, 
under  the  complex  influence  of  the  various  factors  which  can  appear  in  prac- 
tical problems. 

For  controlling  the  hypotheses  and  results  of  such  a calculation  method, 
it  was  felt  necessary  to  carry  out  an  experimental  investigation  giving  suffi- 
ciently detailed  informations  about  the  development  of  mean  and  fluctuating 
characteristics  of  the  boundary  layer  from  the  laminar  to  the  turbulent  regime  ; 
this  study  led  in  fact  to  observations  obviously  concerned  with  the  fundamentals 
of  the  transition  process,  and  involved  as  a logical  consequence,  special  and 
we  hope,  interesting  developments. 

I - PREDICTION  OF  TRANSITION  BY  MODEL  EQUATIONS. 

On  the  theoretical  side,  the  technique  makes  use  of  transport  equations 
for  different  turbulent  quantities,  like  turbulence  intensity  K,  turbulent 
shear-stress  u'v',  and  the  dissipation  of  turbulence  e. 

A two  equations  model  (K,  e)  was  more  intensively  developed;  the  mode- 
lization  is  based  in  a first  stage  on  the  structure  of  established  turbulent 
shear  flows  ; empirical  functions  using  essentially  the  properties  of  the  vis- 
cous wall  region  of  a turbulent  boundary  layer,  are  introduced  after  that  for 
taking  account  of  the  effects  of  viscosity  upon  turbulence. 

Very  encouraging  results  were  obtained  by  showing  at  first  that  the  metho< 
is  actually  able  to  predict  a damping  for  a given  perturbation  if  it  is  intro- 
duced at  a sufficiently  low  Reynolds  number,  and  after  that  amplification  at  a 
Reynolds  number  corresponding  to  the  experimental  transition  Re-number  (fig.  1). 


Different  applications  have  been  made,  which  led  to  a good  agreement 
with  experiments  for  the  onset  of  transition  as  well  as  for  its  development  up 
to  the  turbulent  regime  for  problems  like  : 
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- transition  as  function  of  external  turbulence  on  a flat  plate  in  incom- 
pressible flow  (fig.  2), 

- combined  effects  of  pressure  gradient  and  external  turbulence, 

- combined  effects  of  pressure  gradient,  turbulence  and  wall  heat  transfer 

(calculation  of  transitional  boundary  Layers  on  the  upper  and  lower  sur- 
faces of  turbine  blades,  fig.  3). 

2 - FUNDAMENTAL  EXPERIMENTAL  STUDY  OF  TRANSITION. 

On  the  experimental  side,  the  transition  on  a cylinder  aligned  along  the 
axis  of  an  axisymmetric  channel  has  been  studied  (Velocity  : 33  m/s  ; zero  pres- 
sure gradient  ; external  turbulence  : 0,25  Z) , involving  at  first  a detailed 
determination  of  the  mean  velocity  and  fluctuations  profiles  in  the  transition 
region. 

A good  agreement  with  a calculation  by  the  two  equations  model  was  obtai- 
ned for  the  velocity  profiles  and  mean  flow  characteristics  of  the  boundary 
.layer  (fig.  4).  A general  coherence  was  observed  also  as  concerns  the  develop- 
ment of  turbulence,  but  the  experiments  revealed  complex  phenomena,  especially 
related  with  a marked  intermittency  that  the  present  calculation  model  is  ob- 
viously unable  to  represent. 

A careful  investigation  could  show  in  evidence  Tollmien-Schlichting  waves 
and  their  amplification  ; a spectral  analyzis  of  the  fluctuation  led  to  interes- 
ting observations,  by  showing  notably  the  coherence  between  the  properties  of 
the  waves  observed  and  the  results  of  instability  theory  (fig.  5). 

A specially  significant  aspect  of  the  study  is  concerned  with  the  forma- 
tion of  turbulent  spots  and  with  intermittency.  A systematic  conditional  sampling 
of  the  signal  delivered  by  the  hot  wire  (fig.  6)  was  performed,  and  led  to  de- 
tailed informations  for  : 

- the  distribution  of  the  intermittency  factor  (fig.  7), 

- the  passage  velocities  of  successive  laminar  and  turbulent  flows, 

- the  classification  of  turbulent  spots  and  the  establishment  of  well 
defined  characteristics  for  the  formed  spots, 

- the  determination  of  a statistical  type-picture  for  the  sequence  of 
laminar  and  turbulent  flows  which  could  lead  to  a modelization  of  the 
intermittency  phenomenon  for  future  analytical  developments. 
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Damping  and  amplification  of  turbulent  kinetic  energy  in  a 
laminar  boundary  layer  (Tu  = 0,02  %)  Calculation  by  means 
of  the  K,  e model  method.  e 
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4 J } E*parim«nts 
- Calculation  (k.S  model) 


FIG.  2 * Displacement!  thickness  Reynolds  Number  at  beginning  of 
transition  (flat  plate  in  Incompressible  flow). 
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Section  III 


TURBULENCE-MODEL  PREDICTED  EFFECTS  OF  SURFACE  HEAT  TRANSFER 
ON  BOUNDARY-LAYER  TRANSITION 
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III-2  TURBULENCE-MODEL  PREDICTED  EFFECTS 
OF  SURFACE  HEAT  TRANSFER  ON 
BOUNDARY-LAYER  TRANSITION 

by 

D.  C.  Wilcox  and  T.  L.  Chambers 
DCW  Industries,  Inc. 

Sherman  Oaks , California 

Effects  of  surface  heat  transfer  on  boundary-layer  transition  have 
been  analyzed  using  a second-order  closure  turbulence  model.  The 
primary  objective  of  the  program  has  been  to  use  the  model  to 
analyze  effects  of  surface  heating  on  hydrodynamic  boundary-layer 
transition.  As  a secondary  objective,  the  model  has  been  used  to 
study  transition  sensitivity  to  suction,  surface  roughness  and 
freestream  turbulence. 

With  no  transition-specific  modifications,  the  turbulence  model 
predicts  salient  features  of  incompressible  boundary-layer  transi- 
tion including  sensitivity  to  freestream  turbulence  and  surface 
roughness.  As  shown  in  Figures  1 and  2,  in  addition  to  accurately 
predicting  effects  of  freestream  turbulence,  the  model  accurately 
simulates  transition  width,  skin  friction  overshoot  and  transi- 
tional velocity  profiles. 

With  transition  modifications  based  on  linear  stability  theory, 
the  model  accurately  predicts  transition  sensitivity  to  surface 
heat  transfer,  pressure  gradient  and  suction.  Figure  3 compares 
predicted  and  measured  transition  Reynolds  numbers  for  surface 
roughness  effects  and  for  the  latter  three  effects  on  aerodynamic 
boundary  layers.  Of  particular  interest  to  the  main  project 
objective,  consistent  with  measurements,  the  model  predicts  that 
heating  destabilizes  an  aerodynamic  boundary  layer. 


Figure  4 presents  results  of  computations  for  four  hydrodynamic 
bodies  with  and  without  surface  heating.  Again  consistent  with 
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measurements,  the  model  Indicates  that  small  heating  rates  have  a 
strong  stabilizing  effect  on  hydrodynamic  boundary  layers. 

In  summary,  the  model  accurately  simulates  the  various  effects 
pertinent  to  hydrodynamic  boundary-layer  transition.  The  primary 
unresolved  question  is  the  model's  sensitivity  to  freestream  turbu- 
lence and  the  model's  Inability  to  dJ  :tly  account  for  Tollmien- 
Schlichting  waves.  These  problems  a:  undoubtedly  tied  to  the 
masking  of  spectral  information  atter.  i ir.V,  long-time  (Reynolds) 
averaging.  While  the  method's  utility  for  many  engineering  appli- 
cations does  not  appear  to  be  significantly  hindered  by  these 
unresolved  Issues,  a hybrid  linear-stability/turbulence-model  (in 
contrast  to  linear-stability/e9)  simulation  of  transition  may  pro- 
vide a more  physically  sound  method  than  is  currently  provided 
either  by  the  model  or  by  linear-stability  methods. 
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Section  III 


ON  THE  APPLICATION  OF 
MODELS  TO  BOUNDARY 


SECOND-ORDER  CLOSURE 
LAYER  TRANSITION 
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III-3  ON  TIIE  APPLICATION  OF  SECOND-ORDER  CLOSURE 
MODELS  TO  BOUNDARY  LAYER  TRANSITION 


Michael  L.  Finson 

Physical  Sciences  Inc. 
Woburn,  MA  01801 


Second-order  closure  models  have  been  applied  to  a wide  variety 
of  turbulent  flows  in  recent  years.  There  are  several  advantages  in 
applying  such  methods  to  transition  problems.  One  obvious  advantage  is 
that  second-order  models  are  non-linear,  and  may  provide  useful  results 
in  the  later  stages  of  the  transition  process  as  the  flow  becomes  fully 
turbulent.  Also,  the  governing  equations  can  be  solved  with  straightforward 
numerical  techniques,  even  in  rather  complicated  situations.  However, 
it  has  yet  to  be  demonstrated  that  the  required  closure  approximations 
can  be  made  in  a manner  that  is  physically  sound  for  transition  regimes. 

For  example,  the  second-order  models  cannot  be  expected  to  reproduce, 
linear  stability  theory.  It  might  be  anticipated  that  second-order  models 
would  have  greatest  applicability  in  "by-pass"  situations,  where  transition 
is  triggered  by  relatively  high  levels  of  initial  fluctuations,  introduced  by 
free -stream  turbulence  or  surface  roughness. 

Our  studies  have  utilized  a "five  equation"  model-;  the  dependent 

fluctuating  quantities  are  the  three  components  of  the  fluctuating  kinetic 

• 7 IT  ’"To  i 1 r 

energy  u , v , w , the  Reynolds  stress  u v , and  either  the  dissipation 

rate  $ or  the  length  scale  A*  A typical  equation,  for  the  kinetic  energy 

= (u'2  + v '2  + w'2)/2,  is 


D gl- 

ut 


f t 
U V 


au 

ay 


JL 

ay 


(V  + c 


2 2 
q v 


- 2 
An__' 

dy 


(1) 


where  the  terms  represent  convection,  production,  dissipation,  and  molecular 
and  turbulent  diffusion.  The  required  closure  approximations  for  turbulent 
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flow:;  have  been  developed  from  basic  laboratory  turbulence  measure- 
ments, wherever  possible. 

A careful  examination  of  the  adequacy  of  the  various  closure 

requirements  for  the  transition  regime  has  been  performed.  The 

production  terms  are  treated  exactly  in  our  formulation  (such  is  not  the 

case  with  "2 -equation"  models),  and  so  are  the  molecular  diffusion 

terms.  Concerning  the  dissipation,  we  claim  (without  proof,  here)  that 

adequate  treatments  can  be  developed,  since  exact  solutions  are  known 

in  the  limit  Re  -*  0.  

1 5 p' 

The  difficult  terms  are  the  "redistribution"  terms  (u. ) that 

l ox. 

J 

appear  in  the  equations  for  the  individual  components  of  the-  Reynolds 

stress  tensor.  Turbulent  closure  approximations  have  been  developed 
1 2 

by  Rotta  and  Hanjalic  and  Launder  , but  these  approximations  become 
suspect  near  solid  walls,  for  length  scales  that  are  not  small  compared  to 
mean  flow  dimensions,  or  for  non -turbulent  spectra.  Such  terms  play  an 

important  role  in  the  production  of  fluctuating  energy.  With  the  direct 

'2  i 1 j u ' ~ 

production  terms,  u is  produced  by  u v — — and  u v is  produced  by 

v — — but  (under  the  boundary  layer  approximations)  there  are  no  terms 

® y qT  i 2 

to  produce  v or  w directly.  The  direct  production  terms  will  not  yi  d 

. nr  i ~ 

transition  since  v is  not  produced  and  the  production  of  u v will  hence  be 
small.  However  the  turbulent  closure  for  the  p terms  introduces  additional 
terms  of  the  form  u:  u.  . These  terms  tend  to  isotropize  the  production 

1 j ay  • 

process,  and  cause  an  extremely  rapid  amplification  in  a transition  calculation. 
In  reality,  the  pressu re -strain  terms  must  be  considerably  smaller  in  the 
transition  regime. 

Several  potential  modifications  to  the  pressure -strain  terms  have 
been  investigated.  First  there  are  wall  effects,  which  are  important  in  a 
transitional  boundary  layer  or  in  the  viscous  sublayer.  If  A/y  is  not  small, 
the  normal  fluctuations  v'  arc  reduced  by  the  presence  of  the  wall.  To 
obtain  the  proper  behavior  as  y ->»  0,  the  turbulent  closure  terms  for 


2 2 

pressure -strain  should  be  reduced  by  (1  I A /y  ) • Second,  the 

pressure -strain  terms  arc  generally  the  first  terms  in  a scries  expansion 
in  A/6.  Additional  terms  such  as 


h u.  u. 
i 1 


become  important  in  transitional  situations  where  A/6  is  not  necessarily 


small.  Note  that  this  term  introduces  a dependence  on 


The  net  result 


of  these  modifications  is  a more  reasonable  production  cycle  in  the  transi- 
tion regime.  However,  it  should  be  emphasized  that  the  present  treatment 
is  rather  tentative  and  cannot  yet  be  considered  adequate. 

Figure  1 compares  the  results  of  the  present  model  for  the  effect 
of  free -stream  turbulence  on  flat  plate  transition  with  the  data  of  various 
investigators.  In  the  calculations  the  level  of  the  free -stream  turbulence 
(assumed  isotropic)  was  specified  as  a boundary  condition,  at  y , and  the 
length  scale  in  the  free -stream  was  taken  to  be  sufficiently  large  that  the 
turbulence  there  would  not  decay  over  the  distances  considered.  The  results 
implicate  a somewhat  more  rapid  transition  than  observed,  but  the 
comparison  should  be  qualified  by  the  uncertain  nature  of  the  free -stream 

turbulence  in  most  tests  (probably  acoustic  rather  than  vortical). 

3 

Spangler  and  Wells  observed  a smaller  transition  Reynolds  number, 
for  a given  rms  fluctuation  level,  when  the  free -stream  turbulence  was 
generated  by  a grid,  consistent  with  our  results. 

The  effect  of  distributed  surface  roughness  on  transition  has  also 
been  examined.  In  so  doing,  we  modeled  the  manner  by  which  roughness 
elements  introduce  disturbances  into  the  flow  by  assuming  that  the  flow 
around  individual  elements  is  attached.  A distributed  sink  term  (for  y ik) 
is  included  in  the  mean  momentum  equation  to  account  for  the  drag  of  the 
elements,  and  a distributed  source  term  in  the  kinetic  energy  equation  describes 
the  introduction  of  fluctuations  downstream  of  each  clement.  The  details 
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inay  be  found  in  Ref.  4.  Figure  2 shows  the  resulting  computations 

5 

for  transition,  compared  with  Feindt's  data  on  flat  plates  roughened 
with  sandpaper.  At  roughness  Reynolds  numbers  below  about  250,  free- 
stream  turbulence  appears  to  be  responsible  for  transition. 

The  model  described  here  shows  promise  for  predicting  transition, 
at  least  for  some  conditions  such  as  those  where  there  are  relatively  high 
levels  of  initial  disturbances.  As  indicated  above,  the  closure  approxi- 
mations required  for  the  pressure-strain  terms  are  critical  in  the 
transition  region.  Our  current  treatment  can  only  be  considered  to  be 
an  initial  attempt,  and  further  work  is  required  to  assure  a physically 
sound  closure. 
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IV-1  AN  ANALYTICAL  STUDY  OF  THE  EFFECT 
OF  SURFACE  ROUGHNESS  ON  THE  STABILITY 
OF  A HEATED  WATER  BOUNDARY  LAYER* 


By 

Denny  R.  S.  Ko  and  M.  Kosecoff 

Physical  Dynamics,  Inc. 
Torrance,  Calif. 


This  study  analyzes  the  effects  of  surface  roughness  on  low  speed  boundary 
layer  stability  and  transition  in  a liquid  environment.  An  analytical  model 
for  distributed  surface  roughness  elements,  spaced  closely  compared  to  the 
characteristic  wavelength  of  the  boundary  layer  disturbances,  is  presented. 

The  model  postulates  an  enhanced  momentum  and  heat  exchange  near  the  wall 
induced  by  the  presence  of  these  roughness  elements.  A "turbulent  roughness 
layer"  with  a thickness  of  the  order  of  the  mean  roughness  height  is  assumed 
to  be  imbedded  within  the  ordinary  laminar  boundary  layer.  Enhanced  momentum 
and  energy  exchange  within  this  layer  is  modeled  by  an  effective  eddy  diffu- 
sity.  The  presence  of  the  turbulent  roughness  layer  modifies  the  mean  flow 
(velocity  and  temperature)  profiles  which,  in  turn,  alter  the  stability 
characteristics.  The  effect  on  boundary  layer  transition  is  then  inferred 
by  adopting  a semi-empi ri cal  en  criterion. 

To  demonstrate  the  effect  of  the  distributed  surface  roughness  on  a 
heated  water  boundary  layer,  the  model  is  incorporated  into  a set  of  similarity 
equations  for  estimating  the  mean  flow  profiles.  These  profiles  are  then 
used  in  linear  stability  calculation,  which,  when  coupled  with  the  e11  criteria, 
gives  the  transition  location.  The  linearized  stability  equations  use  a 
parallel  flow  assumption  and  incorporate  the  variation  of  viscosity  and  con- 
ductivity with  temperature.  These  equations  are  solved  using  a code  developed 
by  Lowell  and  Reshotko.  Amplification  rates  are  calculated  as  functions  of 
k/e  and  AT  for  fixed  values  of  Re.  and  dp/dx,  where  k is  the  roughness  height, 
e is  the  momentum  thickness,  and  aT  is  the  difference  between  the  wall  and 

W 

ambient  temperatures. 


* This  work  is  supported  by  the  DARPA  under  ONR  Contract  Number  N00014-76-C-0967. 
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For  a given  dp/dx,  it  is  found  that  there  exists  a range  of  values  of  k/e 
for  which  the  wall  may  be  considered  essentially  smooth.  However,  the  more 
favorable  the  pressure  gradient,  the  narrower  this  range  is.  As  anticipated, 
heating  and/or  a favorable  pressure  gradient  tend  to  stabilize  the  boundary 
layer  for  small  roughness,  and  these  effects  are  reduced  for  moderate 
roughness.  However,  when  k/e  is  sufficiently  large,  heating  tends  to  de- 
stabilize the  boundary  layer.  Furthermore,  as  the  pressure  gradient  becomes 
more  favorable,  the  criteria  for  "large  k/e"  reduces  significantly.  Finally, 
roughness  substantially  increases  the  range  of  frequencies  which  are  unstable 
at  a given  ReQ.  A semi -quantitative  summary  of  these  results  is  given  in 
Table  1. 
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Table  1:  Effect  of  Roughness  on  Amplification  Rates 


Falkner-Skan 
pa rameter.6 

k 

e 

AT  = 0PF 

AT=  30°  F 

0 

< .25 

negligible 

negligible 

.5 

small 

moderate 

.75 

drastic 

drastic 

-0.05 

< .5 

negligible 

negligible 

.75 

moderate 

drastic 

1.0 

drastic 

drastic 

0.2 

< .15 

negl igible 

negligible 

.25 

modera te 

drastic 

.4 

drastic 

drastic 

• .•A'-' 
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Section  IV 


ON  THE  EFFECT  OF  FREESTREAM  TURBULENCE  ON 
BOUNDARY-LAYER  TRANSITION 


l 1 
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IV-2  On  the  Effect  of  Freestream  Turbulence  on  Boundary-Layer  Transition 

Leslie  M.  Mack 
Jet  Propulsion  Laboratory 

The  results  of  several  experiments  are  summarized  in  figure  1 and  show 
a strong  influence  of  freestream  turbulence  on  the  transition  of  a flat-plate 
laminar  boundary  layer.  Most  of  the  data  in  figure  1 are  from  the  Schubauer- 
Skramstad  experiment,  where  tunnel  noise  dominated  transition  below  u^/Uj  of 
about  0.1%.  (u*  and  U are  the  rms  and  mean  x-velocity  components;  the  subscript 

1 refers  to  freestream  conditions.)  Only  the  data  for  Uj //U1  >0.1%  represent 
the  effect  of  freestream  turbulence.  Figure  2 interprets  these  data  in  terms 
of  linear  stability  theory  to  emphasize  just  how  strongly  turbulence  affects 
transition.  The  maximum  amplitude  ratio  of  linear  stability  theory,  in  (A/A^)^  , 
decreases  from  a value  of  8.1  at  the  transition  Reynolds  number,  Re^  = Uxx/v, 
for  ux  '/Ux  = 0.17.,  to  2.6  at  ux  '/Uj  = 17..  Consequently,  if  A at  transition  is 
a constant,  Aq  must  at  the  same  time  increase  by  a factor  of  245.  The  power-law 
relation  Aq  ~ (uj'/Uj)  ’ accounts  quite  well  for  this  behavior  of  Aq. 

The  most  straightforward  way  of  using  linear  stability  theory  as  a means 

of  transition  prediction  is  the  so-called  e9  method  of  Smith-Gamberoni  and  Van 

Ingen.  A value  of  n = in  (A/Aq)  is  chosen,  and  transition  is  considered  to  occur 

whenever  the  maximum  of  in  (A/A  ) at  some  Reynolds  number  equals  n.  Obviously 

o 

a fixed  value  of  n is  incompatible  with  figure  1.  The  easiest  way  to  incorporate 
the  effect  of  turbulence  into  the  e9  method  is  to  let  n be  a function  of  ux  ' /Ux . 
The  data  of  figure  1 are  fitted  closely  by 

n - - 8.43  - 2.4  in  (ux  //U1)  . 
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As  the  above  relation  is  based  solely  on  flat-plate  data,  another 
method  has  been  developed  which  is  intended  to  have  wider  application. 

This  method,  which  can  be  referred  to  as  an  amplitude  method  to  distinguish 
it  from  the  amplitude-density  e9  or  modified  e9  methods,  is  based  on  the 
following  ideas : 

a)  it  uses  the  actual  disturbance  amplitude  which  is  made  up  of  harmonic 
components  of  all  frequencies  and  orientations: 

b)  the  initial  amplitude  density  Ao(<i),  0) , where  to  is  the  frequency  and  0 
the  transverse  wave  number,  is  related  to  the  corresponding  component  of 
the  freestream  turbulence  by  an  interaction  relation; 

c)  isotropic  turbulence  theory  is  used  to  characterize  the  freestream 
turbulence . 

With  the  individual  harmonic  components  of  the  boundary- layer  disturbance 
assumed  to  have  random  phase  relations,  the  disturbance  amplitude  is  given 
as  a function  of  Reynolds  number  by 


where  L is  an  integral  length  scale  of  the  turbulence . If  advantage  is  taken 
of  the  sharp  tuning  of  the  boundary- layer  response,  this  expression  can  be 
simplified  to 


A/(R)  = A 


'(“W  gmax)  fQ  d(St)Z  [ao(v’  P)]  dP' 
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Because  the  function  A/Aq  always  has  roughly  the  same  shape,  a further 
simplification  can  be  made  to 


The  two  factors  A (uw/Uj2)  and  A3  are  consistently  defined  response  band- 
widths,  and  along  with  (A/A Q)max  and  (u)v/U1s)max  are  obtained  from  stability 
theory  with  little  more  computational  effort  than  is  required  for  the  e9 
method.  For  a two-dimensional  incompressible  boundary  layer,  the  maximum 
response  is  always  at  3 = 0*  Since  it  has  been  observed  for  these  same 
boundary  layers  that  A3  is  only  slightly  dependent  on  the  Reynolds  number, 
this  factor  will  be  absorbed  into  the  constant  C* . The  scale  Reynolds 
number  U1L/v  is  set  by  the  freestream  turbulence,  and  AQ  is  to  be  found 
from  an  interaction  model. 

In  the  absence  of  compelling  theoretical  or  experimental  arguments 
as  to  the  mechanism  by  which  freestream  turbulence  produces  laminar  insta- 
bility waves,  the  following  arbitrary  assumptions  are  made: 

a)  instability  waves  originate  near  the  solid  boundary; 

b)  their  amplitude  is  proportional  to  the  normal  viscous  velocity  of  the 
forced  response  in  the  viscous  sublayer  of  the  freestream  turbulence; 

c)  the  forced  response  is  given  by  Prandtl's  model  of  the  sublayer  (boundary- 
layer  equations  with  U * 0) . 
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The  result  of  the  analysis  for  the  rms  normal  viscous  velocity  vy  ' of 
component  u>,  8 is  given  by 


where  i|r  is  the  angle  between  the  wave-number  vector  in  the  plane  of  the 
surface  and  the  freestream  direction.  The  component  rms  pressure  fluc- 
tuation p'(u>,  8)  is  represented  by  a two-dimensional  power-spectrum  function 
F2  0)  • 

p'Cto,  0)  3 = / E_^\g[UiFa(m»  0) 

. P^8  \pU1s/  l L2 


The  wide-band  pressure  fluctuation  p*  is  taken  to  be  that  of  isotropic 
turbulence,  and  is  thus  related  to  the  velocity  fluctuation  by 
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With  Aq  assumed  proportional  to  v *, 
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The  final  expression  for  the  disturbance  amplitude  is,  with  8=0  and 
C - CjCs , 


V 
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The  only  remaining  task  is  to  develop  an  expression  for  the  spectrum 
function  Fa(ou).  The  two-dimensional  wave  number  spectrum  is  given  in  terms 
of  the  three-dimensional  spectrum  by 


Fa(ka,a) 


F3  (k)  dk3 


> 


where  ks  = kx2  + k32  + k32 , and  klj32  = kj2  + k32 . The  one-dimensional 
spectrum  is 


k F3(k)  dk  . 


Consequent ly. 


Mk1>3)  - - 


1 

(k2  - k1>22)* 


dk  . 


The  normalized  one-dimensional  pressure  spectrum  which  corresponds  to  the 
von  Karman  interpolation  formula  for  the  one-dimensional  velocity  spectrum 
la 

FjCkj)  - 4 fl  + 


(Ml 
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Use  of  this  expression  leads  to  the  following  interpolation  formula  for 
**a  (kx  t a)  5 


Fa(k1>a)  = 1.78  [l  + (0.82  k1>a)a] 


Both  Fx  and  Fa  are  shown  in  figure  3.  For  the  8=0  component,  which  is 
the  only  one  used  in  the  present  calculations,  k1>a  can  be  replaced  by 
uuL/Uj  . 

Calculations  have  been  made  on  the  basis  of  the  amplitude  and  modified 
e9  methods  for  Falkner-Skan  boundary  layers  with  the  pressure-gradient  parameter 
8 equal  to  0 (flat  plate)  and  0.1.  All  of  the  numerical  results  were  obtained 
from  the  newly  developed  JPL  stability  and  transition  computer  program  BLSTR. 

This  program  starts  with  a single  damped  eigenvalue,  and  first  locates  the 
unstable  region.  Then,  as  it  advances  downstream,  it  adds  and  subtracts 
frequencies  as  needed  to  cover  only  the  unstable  region  and  keep  the  number 
of  damped  frequencies  to  a minimum.  After  the  eigenvalues  have  been  computed 
at  each  Reynolds  number,  least-squares  fits  to  in  (A/Aq)  are  used  to  evaluate 
in  (A/A^^^,  corresponding  dimensionless  frequency  F‘max  = (cuv/Uj2)^^, 

and  the  bandwidth  AF  defined  as  the  difference  in  F in  the  direction  of  de- 
creasing F such  that  in  (A/A  ) - in  (A/A  ) = 1.  The  interaction  relation 

o max  o 

is  evaluated  for  Fmax  and  the  disturbance  amplitude  computed.  If  the  transition 

criterion  of  A.  = 0.04  is  not  satisfied,  the  program  advances  to  the  next 
a 

Reynolds  number.  Typical  execution  times  on  the  medium-speed  UNIVAC  1108 

computer  are  about  1%  sec  for  each  eigenvalue,  and  also  a total  of  1%  sec  to 

do  all  of  the  transition-prediction  calculations  (neutral  points,  A/A  , (A/A  ) , 

o o max 

Aq,  A^j)  for  a dozen  Reynolds  numbers. 
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The  numerical  results  from  stability  theory  for  the  two  boundary  layers 
are  given  In  figure  4.  It  is  to  be  noted  how  close  to  a straight  line  are 
the  two  envelope  curves  £n  (A/A^)mgy  for  the  larger  amplifications.  Figure  5 
shows  the  disturbance  amplitude  A^  for  the  flat-plate  boundary  layer  as  a 
function  of  Re^  for  several  turbulence  levels  and  UxL/v  = 4 X 104 . The  constant 
C in  the  relation  for  A^  was  evaluated  once  and  for  all  by  requiring  that  A^  = 
0.04  at  Re  = 2.8  X 10°  with  uj  '/Uj  = 0.001.  Finally,  figure  6 gives  Ret  as  a 
function  of  ut  '/Uj  for  8 = 0 and  0.1  according  to  the  two  prediction  methods. 

For  the  amplitude  method,  UjL/v  = 4 X 10*  as  in  figure  5.  Further  comparisons 
with  experiment  are  needed  to  evaluate  the  worth  of  the  two  methods. 

This  paper  presents  the  results  of  one  phase  of  research  carried  out  at 
the  Jet  Propulsion  Laboratory,  California  Institute  of  Technology,  under  Contract 
No.  NAS7-100,  sponsored  by  the  National  Aeronautics  and  Space  Administration. 


Jet  Propulsion  Laboratory 
Pasadena,  California  91103 

September  10,  1976 
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STABILITY  THEORY  AND  FLAT-PLATE  TRANSITION  EXPERIMENTS 
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PRESSURE  FLUCTUATION 
INTERPOLATION  SPECTRA 
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EFFECT  OF  FREESTREAM  TURBULENCE 
ON  TRANSITION-THEORY 
AND  EXPERIMENT 
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Section  IV 

THE  EFFECTS  OF  FINITE-AMPLITUDE  DISTURBANCES 
ON  THE  MEAN  FLOW  PROFILES  IN  A LAMINAR  BOUNDARY  LAYER 
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The  use  of  the  linear  stability  characteristics  of  boundary  layers  to 
indicate  the  location  of  transition  from  laminar  to  turbulent  flow  has 
proven  to  be  one  of  the  best  available  methods  for  predicting  transition 
on  bodies  of  general  engineering  interest.  The  underlying  reason  for 
this  success  is  almost  certainly  related  to  the  superior  theoretical 
foundation  of  the  stability  approach  as  compared  to  that  of  other  exist- 
ing methods.  Despite  the  relative  success  of  the  method,  the  actual 
relationship  between  stability  and  transition  is  poorly  understood,  and 
is  generally  incorporated  in  the  predictions  in  an  empirical  manner  in 
conjunction  with  an  implicit  assumption  that  the  ensuing  nonlinear  ampli- 
fication region  which  follows  the  linear  growth  region  is  very  brief 
and  leads  quickly  to  transition.  The  most  stringent  limitation  on  this 
approach  arises  because  there  is  no  mechanism  in  the  linear  theory  for 
deducing  when  the  disturbances  become  large  enough  that  nonlinear 
effects  begin  to  be  important.  The  present  approach  is  directed  towards 
including  the  appropriate  additional  terms  in  the  stability  analysis 
so  that  the  effects  of  the  disturbances  and  their  growth  can  be  deduced 
directly  from  the  stability  analysis  itself  without  recourse  to  addi- 
tional considerations.  Since  only  the  leading-order  finite-amplitude 
effects  are  included,  the  analysis. still  requires  an  empirical  calibra- 
tion to  relate  the  finite-amplitude  effects  to  the  observed  location 
of  transition.  The  present  approach  also  requires  a knowledge  of  the 
characteristics  of  the  initial  disturbance  levels  inside  the  boundary 
layer,  but  their  origin  and  character  are  not  addressed  in  the  current 
study. 

The  extension  of  stability  theory  to  include  finite-amplitude  disturbances 
starts  from  a formal  perturbation  expansion  which  includes  the  complete 
nonlinear  stability  characteristics  of  the  boundary  layer.  This  expan- 
sion procedure  demonstrates  that,  for  a single-frequency  disturbance, 
the  leading-order  effects  of  a finite-amplitude  wave  result  in  a distortion 


I 
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of  the  mean  flew  profile.  In  particular,  the  distortion  of  the  mean 
flow  profile  is  a phenomenon  which  is  of  second  order,  0 fa2),  in  the 
disturbance  amplitude,  whereas,  the  nonlinear  disturbance-disturbance 
interaction  is  of  third  order,  0(A?).  These  observations  suggest 
that  the  distortion  of  the  mean  flow  profile  can  be  accurately  computed 
while  still  retaining  a linear  equation  system  for  the  disturbances. 
Extension  of  the  analysis  to  more  realistic  multi-frequency  and  three- 
dimensional  disturbance  environments  shows  that  although  the  nonlinear 
disturbance-disturbance  interactions  now  become  of  the  same  order  as 
the  mean  flow  distortion  (namely,  order  A2),  the  second-order  accurate 
description  of  the  mean  flow  can  still  be  obtained  from  a first-order 
accurate  (linear)  representation  of  the  disturbances  (i.e.,  second- 
order  changes  in  the  disturbance  growth  rate  cause  thi rd-order  changes 
in  the  mean  flow  profiles).  Such  a simplification  to  a linear  wave 
system  is,  of  course,  fundamental  to  the  analysis  from  the  viewpoint 
of  computational  economy.  Further,  the  concept  that  the  mean  flow  profile 
distortions  can  be  accurately  computed  from  a linear  disturbance  analysis 
is  in  general  agreement  with  experimental  evidence  which  has  indicated 
that  the  generation  of  harmonic  disturbances  is  of  only  minor  signifi- 
cance in  the  transition  process. 

In  contrast  to  these  results  for  boundary  layers,  the  analogous  deduc- 
tions for  the  more  intensively-studied  Poiseui lie-flow  case  appear  to 
be  considerably  different.  Although  an  identical  expansion  scheme  has 
been  used  for  most  nonlinear  stability  studies  in  Poiseuille  flows, 
the  attention  in  this  latter  family  of  flows  has  been  almost  unanimously 
focused  on  predicting  the  changes  in  the  growth  rate,  and  has  required 
a knowledge  of  both  the  mean  flow,  distortion  and  the  wave-wave  interaction. 
Specifically,  these  analyses  have  been  carried  to  third  order,  0 (A3), 
accuracy.  The  differences  in  the  emphasis  of  the  nonlinear  approach  in 
these  two  situations  comes  not  from  the  theory  itself,  but  from  the 
comparison  of  the  theory  with  experiment.  In  boundary  layers,  the  pre- 
dicted region  in  which  disturbances  are  amplified  is  in  very  good  agree- 
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ment  with  extensive  experimental  evidence.  For  Poiseuille  flow,  transition 
to  turbulence  occurs  below  the  critical  Reynolds  number  at  which  distur- 
bances first  begin  to  be  amplified,  implying  that  amplification  can  take 
place  in  this  low  Reynolds  number  regime.  As  a result  of  this  discre- 
pancy, the  major  impetus  of  nonlinear  stability  studies  in  Poiseuille 
flows  has  been  towards  modifying  the  qualitative  results  of  linear 
stability  theory  and  using  the  nonlinear  terms  to  show  that  regions  which 
are  stable  to  infintesimal  linear  waves  are  unstable  to  finite-amplitude 
waves.  Such  a qualitative  change  in  the  stability  properties  of  boundary 
layers  is  not  needed,  and  would  disagree  with  experiment  if-found. 

Some  initial  numerical  results  which  are  based  on  single  two-dimensional 
disturbances  in  boundary  layers  are  presented  to  demonstrate  the  types 
of  mean  flow  distortions  which  can  be  generated  by  finite-amplitude 
disturbances.  The  curvatures  of  the  velocity  profiles  'have  been  used 
as  particularly  sensitive  indices  of  the  cumulative  effects  of  the  inter- 
action between  the  disturbance  and  the  mean  flow.  In  the  presence  of 
finite-amplitude  disturbances,  the  curvature  near  the  wall  decreases 
while  that  farther  from  the  wall  increases,  eventually  leading  to  a 
point  of  inflection.  Comparisons  of  zero-pressure-gradient  boundary 
layers  with  favorable  and  unfavorable  pressure  gradient  cases  indicates 
that  significant  distortion  begins  at  lower  values  of  the  disturbance 
amplitude  when  the  pressure  gradient  is  favorable  as  compared  to  when 
it  is  unfavorable;  however,  at  all  times,  substantial  interaction  occurs 
while  the  disturbance  is  very  small  (typically  u'v'  less  than  10~5  times 
the  free-stream  velocity  squared).  The  current  results  suggest  that  either 
the  velocity  profile  curvature  or  the  ratio  of  the  disturbance-generated 
Reynolds  stresses  to  the  local  shear  stress  could  serve  as  effective 
practical  indicators  of  transition.  Finally,  we  note  that  the  use  of 
a transition  indicator  which  is  based  on  the  distortion  of  the  mean 
profiles  by  the  disturbances  in  the  laminar  boundary  layer  seems  par- 
ticularly attractive  for  applications  such  as  advanced  laminar  flow 
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vehicles  where  extrapolation  of  the  current  transition  data  base  to 
higher  Reynolds  numbers,  through  the  use  of  pressure  gradient,  heating 
and  suction  is  required.  The  reduced  degree  of  reliance  on  the  empirical 
information  in  the  transition  analysis  should  increase  the  reliability 
of  this  extrapolation. 
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Section  IV 


FLUCTUATIONS  IN  A BOUNDARY  LAYER  INTRODUCED  BY 
TRAVELING-WAVE  IRROTATIONAL  FREESTREAM  DISTURBANCES 
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IV-4  FLUCTUATIONS  IN  A BOUNDARY  LAYER  INTRODUCED  BY 
TRAVELING-WAVE  IRROTATI ONAL  FREESTREAM  DISTURBANCES* 


Harold  L.  Rogler 

Department  of  Mechanical  and  Aerospace  Engineering 
Case  Western  Reserve  University;  Cleveland,  Ohio  44106 


INTRODUCTION  TO  IRR0TATIONAL  DISTURBANCES 

If  you  wave  your  hand  above  your  desk,  the  fluctuations  which  are 
generated  arise  from  the  bound  circulation  associated  with  any  lift  produced 
by  your  hand  (which  might  be  represented  as  a sheet  of  vorticity  or  a 
potential  vortex),  from  the  thickness  of  your  hand  (which  might  be  represented 
as  a distribution  of  sources  and  sinks),  and  from  the  trailing  wake  of 
vorticity  shed  behind  your  hand  as  you  unsteadily  wave  and  alter  that  lift 
(which  might  be  modeled  as  a vortex  sheet  of  varying  strength  or  a set  of 
free  vortices).  The  disturbances  introduced  by  the  thickness  and  bound 
circulation  effects  are  everywhere  irrotational,  while  those  introduced  by 
the  trailing  vorticity  are  effectively  irrotational  at  some  distance  away 
from  the  patches  or  layers  of  vorticity. 

Irrotational  fluctuations  are  also  generated  by  a rotor  blade  translating 
past  a stator  blade  in  a jet  engine,  by  wavy  walls,  or  by  the  small  amplitude 
waves  at  the  surface  of  the  ocean.  A fish  swimming  through  rapids  is 
subjected  to  irrotational  fluctuations  as  it  swims  through  water  which  is: 

(1)  flowing  around  rocks,  (2)  influenced  by  standing  vortices  or  eddies 
behind  those  rocks,  or  (3)  influenced  by  eddies  which  are  being  convected 
downstream.  Of  course  it  is  possible  that  the  fish  swim  (or  airplanes  fly) 
through  a patch  of  vorticity,  but  we  are  not  concerned  in  this  study  with  the 
special  effects  arising  from  rotational  disturbances. 

In  all  these  cases,  the  disturbances  are  induced  by  some  surface  or  body 
in  the  flow,  by  some  arrangement  of  sources,  sinks,  doublets,  monopoles,  etc., 
by  interfaces  between  stratified  fluids,  by  bound  circulations,  and  by 
patches,  sheets,  rings,  or  layers  of  vorticity.  If  a plate  is  inserted  into 
a flow  with  freestream  vorticity  fluctuations  and  if  these  vorticity 


*This  work  was  supported  by  the  Air  Force  Office  of  Scientific  Research  under 
Grants  AFOSR-74-2577A&B.  The  support  and  encouragement  by  Dr.  E.  Reshotko 
are  gratefully  acknowledged.  The  author  thanks  Dr.  J.  Paul  for  the  use  of 
his  streamline  plotting  program. 
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fluctuations  are  sufficiently  far  away  from  the  plate  such  that  the  vorticity 
does  not  convect  into  or  diffuse  into  the  boundary  layer,  then  as  far  as  the 
boundary  layer  is  concerned,  these  disturbances  are  irrotational. 

HOW  FAST  DO  IRROTATIONAL  DISTURBANCES  PROPAGATE? 

The  phase  speeds  of  irrotational  fluctuations  are  intimately  connected  to 
the  propagation  speed  of  the  source.  If  that  source  is  the  traveling,  wavy- 
wall  of  a wind  tunnel,  then  the  fluctuation  propagates  at  that  speed  also. 

In  turbomachinery,  the  relative  motion  between  blades  would  be  important  as 
well  as  the  convection  of  vortex  sheets  downstream.  If  the  source  is  a set 
of  low- intensity  vortices  being  convected  with  the  mean  flow,  then  the 
propagation  speed  is  the  local  mean  velocity  at  the  vortex  centers.  If  the 
source  is  a vortex  street  downstream  of  a bluff  body,  then  the  nonlinear 
interactions  between  the  vortices  cause  the  vortices  to  propagate  at  a speed 
either  slower  or  faster  than  the  frees tream,  depending  on  the  case  considered. 
It  is  possible  that  the  disturbances  propagate  upstream  as  well  as  downstream, 
or  they  can  be  standing  waves  and  not  propagate  at  all.  The  response  in  the 
boundary  layer  is  very  sensitive  to  the  phase  speed,  at  least  for  some  wave- 
numbers  and  Reynolds  numbers. 


INVISCID  INTERACTION  OF  AN  IRROTATIONAL  DISTURBANCE  WITH  A SEMI- INFINITE  PLATE 
If  the  freestream  disturbance  has  the  normal  velocity 


v(a) (x,y,t;k,c) 


0(k,c)ekyeik(x-Ct) 


(1) 


with  (real  and  positive)  wavenumber  k and  phase  speed  c,  and  a semi- inf inite 
plate  is  Inserted  into  the  flow  along  the  x-axis  with  leading  edge  at  the 
origin,  then  by  conformal  mapping  the  inviscid,  irrotational  solution  for  all 
amplitudes  can  be  found  which  satisfies  impermeability  along  the  plate.  The 
streamline  pattern  for  this  flowfleld  is  shown  in  Figure  1.  This  pattern 
reflects  only  the  disturbance  flow;  the  uniform  mean  flow  from  left  to  right 
Is  not  Included. 

Two  features  are  important  here.  First  note  that  the  fluctuations 

vanish  beneath  the  plate.  As  far  as  the  boundary  layer  developing  on  the 

bottom  side  of  the  plate  is  concerned,  any  Influence  of  the  irrotational 

kv+ik (x-ct) 

fluctuation  generated  above  the  plate  of  the  form  e 7 very  probably 

arises  near  the  leading  edge.  Far-downstream  of  the  leading  edge  on  the 
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Figure  1.  Disturbance  streamlines  of  art  irrotational  free- 
stream  disturbance  generated  by  some  source  far- 
above  the  plate  which  propagates  along  and  inter- 
acts with  a semi-infinite  plate.  The  source  might 
be  a vortex  sheet  or  a wavy  wall. 


<• 


-128- 


bottom  side,  there  is  no  fluctuation.  At  the  leading  edge,  a singularity 
arises  with  infinite  speeds.  Secondly,  note  on  the  top  side  of  the  plate 
that  a pattern  rapidly  develops  which  subsequently  propagates  along  without 
further  significant  change.  This  solution  is  called  the  "far-downstream" 
solution  and  develops  within  a wavelength  of  the  leading  edge.  On  the  bottom 
side  of  the  plate,  the  analogous  far- downstream  solution  is  identically  zero. 

To  distinguish  between  whether  you  have  a far-downstream  solution  as  on 
the  top-side  or  on  the  bottom-side,  you  must  know  which  side  of  the  plate  the 
source  of  the  fluctuation  lies.  It  is  not  adequate  to  merely  know  the 
amplitude  of  a sinusoidal  freestream  disturbance  along  the  x-axis.  If 
vorticity  disturbances  or  wavy  walls  exist  both  above  and  below  the  plate, 
then  as  far  as  the  far-downstream  solutions  are  concerned,  you  must  sort  out 
the  contributions  arising  from  above  and  below  the  plate. 

The  singularity  at  the  leading  edge  and  the  complications  of  disturbances 
generated  by  sources  on  the  opposite  side  introducing  fluctuations  near  the 
leading  edge  as  discussed  before  can  be  avoided  by  modifying  the  freestream 
disturbance.  If  the  irrotational  fluctuation  consists  of  a combination  of 
growing  and  decaying  waves  (in  the  y-direction) 


a.  ky  -ky. 
v(e  - e ) 


ik(x-ct) 


= v sinh(ky)  e 


ik(x-ct) 


as  shown  in  Figure  2,  then  the  "far-downstream"  solution  again  appears  except 
now  it  is  valid  everywhere  along  the  x-axis,  even  upstream. 

While  perhaps  appearing  at  first  as  a restrictive  case  of  irrotational 
fluctuations,  eqn.  (2)  is  the  far-downstream  behavior  of  many  disturbances  of 
scientific  and  engineering  significance.  It  is  the  irrotational  freestream 
disturbance  which  will  be  used  to  Introduce  fluctuations  in  the  boundary  layer. 


FORMULATION  FOR  THE  RESPONSE  IN  A BOUNDARY  LAYER  TO  IRROTATIONAL  FLUCTUATIONS 
If  the  velocities  and  pressure  are  separated  into  mean  and  disturbance 
components  in  the  Navier-Stokes  and  continuity  equations  for  a 2-D,  incom- 
pressible flow,  and  the  disturbance  velocities  and  pressure  are  further 
separated  into  the  two  parts 

i (a)  i (a)  , , (a) 

v •«  v -v,  u = u - u,  and  p'  **  p - p 

t — change  to  the  freestream  disturbance  caused  by  the  boundary 
layer  and  no- slip 

freestream  disturbance  which  exists  in  the  absence  of  the  plate 
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the  parallel-flow  representation  of  the  boundary  layer  is  introduced,  and  the 
results  are  linearized,  then  the  following  nonhoroogeneous  equation  results 

(a) 

<1?  + ”<>->!;  - v72}'2’'  - Vy)£  ‘ - 0.4  - Vy)*~ 

which  is  of  form  *^v  = £v^3^ 


The  left-hand  side  is  (the  partial-differential  form  of)  the  Orr-Sommerfeld 
operator  acting  on  the  normal  velocity.  This  normal  velocity  represents  the 
change  to  the  freestream  disturbance.  The  right-hand  side  is  a forcing 

function  which  depends  on  the  structure  of  the  freestream  disturbance 

(a)  — — 

v (x,y,t)  as  well  as  the  mean  velocity  profile,  U and  U . This  forcing  is 

generally  non- zero  inside  the  boundary  layer,  but  vanishes  outside.  The 

equation  is  subject  to  the  impermeability  and  no-slip  conditions  (v,=0  or 
(a)  (a) 

v=v  and  9v'/8y=0  or  9v/9y=9v  /9y  at  y*0)  and  the  conditions  that  the 

(a) 

freestream  disturbances  are  recovered  far-away  from  the  plate  (v’-»vv  or  v-K) 
(a) 

and  gv'/By-^v  /9y  or  9v/9y*>0  as  y-*») . Other  conditions  necessary  to  yield 
a well-posed  problem  will  be  handled  by  periodic  conditions  in  x and  t as 
specified  later. 

If  the  freestream  disturbances  are  irrotational,  V2v^=0,  then  the 
forcing  function  simplifies  to  £vv  =-Uy^3v  / 9x.  Thus  the  forcing  depends 

on  the  variation  of  mean  vorticity  in  the  boundary  layerdf (y) /9y=-U^  and  the 
wall  value  of  the  longitudinal  freestream  fluctuation. 

Seeking  solutions  of  the  form  v(x,y,t)*4>(y)exp[ik(x-ct)] where  <f>  is  complex, 
then  the  mathematical  system  which  results  for  the  freestream  disturbance  of 
eqn.  (2)  is 

{ [U(y)  - c]  (D2 -k2)  - U„„(y)  + (D4  - 2k2D2  +k‘4)}<j>(y)  « -U„„(y)sinh(ky) 


yy 


6* 


yy 


(4) 


subject  to  $*0,  D<J>=*k  at  y“0  and  <Jr>0,  D$-+0  as  y-**>. 


where  D=d/dy,  R^Jk=U<jo6*/v,  and  U(y)  is  the  mean  velocity  nondimensionalized 
with  respect  to  U^.  The  dimensionless  wavenumber  is  related  to  the  dimensional 
wavelength  and  displacement  thickness  by  k“2w6*/A.  The  phase  speed  is  also 
nondimensionalized  against  U^. 

In  general,  a numerical  solution  must  be  obtained  because  the  coefficients 
U(y)  and  U (y)  are  variable  and  are  available  only  as  data.  However,  two 
analytical  solutions  have  been  obtained  for  simple  velocity  profiles. 
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EXPLICIT  SOLUTION  WHEN  U=1  EVERYWHERE 

If  the  mean  flow  is  uniform,  as  would  arise  if  the  top  surface  of  the 
plate  were  composed  of  a belt  which  moved  along  at  the  freestream  speed,  then 
the  equation  reduces  to  { (1-c)  (D2-k2)  +~- -(D2-2k2D2+k4)  }<j>*0.  An  explicit 

"ca* 

solution  for  <f>  can  be  found  and  v'={sinh(ky)-iKy) }elk^x_ct^  formed  as 


v ' = ft eal 


sinh(ky)  +- 


rky-(.-(1+1R6*(1-c)/k)iSkn 


1-  (l  + lR^U-O/k)3*  J 


ik(x-ct) 


for  c?4! 


6*  is  a dummy  parameter  kept  only  so  that  this  result  can  be  compared  with 
solutions  where  a boundary  layer  is  present.  When  c=l,  a modified  form  of 
this  solution  results.  The  streamlines  plotted  in  Figure  3 show  that  the 
no-slip  condition  is  satisfied  through  the  creation  of  a vortex  near  the 
wall.  The  importance  of  viscosity  is  dictated  by  the  combination  of 
parameters  RgJt(l-c)/k=R^(l-c)/2rH^>.  Note  the  dependence  on  phase  speed  even 
with  this  simple  mean  velocity  profile.  Hence  disturbances  near  the  plate 
traveling  at  speeds  near  the  freestream,  c=*l,  are  highly  viscous  and  have 
relatively  thick  viscous  layers. 


PARTLY  ANALYTICAL,  PARTLY  NUMERICAL  SOLUTION  WHEN  U-CONSTANT  AND 
U = ANOTHER  CONSTANT 

yy 

If  the  mean  boundary  layer  is  modeled  by 


and 


U(y)  = Uq  and  U^(y)  **  U2  for  y < <5  (inside  the  layer) 

U = 1 and  U = 0 for  y > 6 (exterior  to  the  layer) 


where  Uq  and  U2  are  constants,  then  another  explicit  solution  of  eqn.  (4)  is 
possible  in  the  form 


(<  J A .V> 

j-i  J 


where 


= ± 


,2  + 


iaR6*(VC> 


i°R6*(Do-c) 


“ 4iU2  ~ 

-l  + «W<=>2- 


when  c?*U 


Effect  of  no-s/<t> 

?r.(/-c)/k  = /ooo 


ftno/  srreo/n/tnes 


Above  the  thin  viscous  layer, the 
flow  is  irrotational  V2iJj=0 

1 


V2V2ij(=0  (above  pattern  only) 


Disturbance  streamline  pattern, 
ip,  for  the  flow  established  by 
the  no- si ip  condition 


Final  streamline  pattern  v =tp 
consisting  of  a superposition  of  the 
freestream  disturbance  and  the  flow- 
field  illustrated  on  the  left. 


Figure  3.  Disturbance  streamline  patterns  for  the  flow  near  the  wall 
when  the  mean  velocity  is  uniform  everywhere,  U * 1,  for 
several  values  of  the  parameter  R = Rx.(l-c)/k  ■ R.  (1-c)/2tt. 
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where  the  four  constants  and  the  two  constants  Ci  and  C2  from  the  exterior 
solution  are  evaluated  numerically  from  a 6x6  matrix.  The  system  holds  great 
promise  for  economically  revealing  many  features. 


NUMERICAL  SOLUTIONS  WITH  A BLASIUS  MEAN  BOUNDARY  LAYER 

With  U(y)  and  U^(y)  represented  as  data  for  a Blasius  layer, 
has  the  complexities  and  numerical  instabilities  of  the  eigenvalue 
for  the  stability  of  a boundary  layer.  Here,  a numerical  solution 
obtained  by  an  expansion  in  Chebyshev  polynomials 


eqn.  (4) 

problem 

is 


N-l 

<Kn)  = I a T (n) 

m=0  m m 


where  T (n)  is  the  mth  Chebyshev  polynomial,  a is  the  (unknown)  mth  coefficient, 

2y  m m 

n = — 1 is  the  transformed  y-coordinate  with  the  range  -1  < n < 1 representing 
Ye 

0 < y < yg,  and  y^  is  the  "edge"  of  the  boundary  layer.  A matrix  (often 
45x46)  is  set  up  representing  the  expansions  of  the  differential  equation  and 


conditions.  Solutions  are  obtained  efficiently  for  u' , v' , -u'v',  -u'v'  dU/dy, 
dissipation,  vorticity,  kinetic  energy,  and  pressure.  Streamlines,  equi- 
vorticity  contours,  and  pressure  isobars  are  also  drawn. 


TYPICAL  NUMERICAL  RESULTS 

The  streamline  patterns  for  several  phase  speeds  are  shown  in  Figure  4 on 
the  right-hand  side.  A better  view  of  what's  going  on  follows,  however,  if 
only  the  alteration  to  the  freestream  disturbance  caused  by  the  boundary  layer 
and  no-slip  condition  are  plotted  (the  impermeability  condition  is  satisfied 
by  the  freestream  disturbance  itself).  This  alteration  is  plotted  on  the 
left-hand  side.  For  the  case  of  c=0.35,  a double  pair  of  rotating  vortices 
are  formed,  but  when  c“0.95  the  structure  has  changed  to  a single  pair  of 
vortices.  When  the  freestream  disturbances  convect  at  the  freestream  speed, 
c=l,  the  third  row  of  streamlines  result.  Note  the  change  in  patterns  when 
the  phase  speed  is  Increased  from  0.95  to  1.05. 

The  amplitude  of  the  longitudinal  disturbance  velocity  is  plotted  in 
Figure  5 for  various  phase  speeds.  Again  note  the  sensitivity  to  phase  speed 
near  c*l.  When  c»l,  the  amplitudes  inside  the  boundary  layer  are  small.  The 
Reynolds  stress  is  plotted  in  Figure  6 which  shows  that  a fairly  constant 
stress  is  produced.  This  constant  stress  continues  far  into  the  freestream. 


The  streamlines  for  the  flowfield 
representing  the  alteration  of 
the  frees tream  disturbance  by  the 
no-slip  condition  and  the  mean 
boundary  layer. 


The  superposition  of  the  freestream 
disturbance  of  Fig.  2 and  the  flow- 
field  illustrated  on  the  left  yields 
the  final  disturbance  flowfield  with 
streamlines  as  shown  above. 


Figure  A.  Disturbance  streamlines  with  a Blasius  mean  boundary  layer 
for  various  phase  speeds  and  for  k=0.2  and  R_  -1000. 


0*0001  = 


Figure  5.  Amplitude  of  the  longitudinal  velocity  for  several  phase  speeds  for  an 

irrotational  freestream  disturbance  interacting  with  a Blasius  boundary  layer. 
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where  presumably  it  vanishes  across  another  viscous  layer.  The  energy 
production  and  dissipation  are  plotted  in  Figure  7.  The  production  for  this 
set  of  parameters  (c=0.35,  R^-1000.0  and  k»0.2)  has  a maximum  at  n^l.5, 
while  the  dissipation  is  largest  in  the  viscous  region  near  the  wall. 
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Figure  7.  The  production  and  dissipation  of  disturbance  energy  of  an  irrotational 
freestream  disturbance  interacting  with  a Blasius  boundary  layer. 
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Section  V 

STABILITY  OF  A HEATED  WATER  BOUNDARY  LAYER 
FOR  NON-UNIFORM  WALL  TEMPERATURE  DISTRIBUTIONS 
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FOR  NON-UNIFORM  WALL  TEMPERATURE  DISTRIBUTIONS 
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ABSTRACT 

The  effect  of  non-uniform  wall  temperature  distributions  on  the 
stability  of  a heated  flat  plate  laminar  boundary  layer  in  water  is  studied 
experimentally.  The  experiment  is  performed  in  a closed  circuit  low 
turbulence  water  tunnel  with  free  stream  turbulence  intensities  of 
0.1-0. 2%.  Temperature-compensated  hot  film  anemometry  is  used  to  measure 
the  mean  flow  field  as  well  as  the  spatial  growth  characteristics  of 
sinusoidal  velocity  disturbances  introduced  into  the  boundary  layer  using 
a vibrating  ribbon. 

Two  types  of  non-uniform  wall  temperature  distributions  are  studied: 

step  changes  in  wall  temperature  of  the  form  T -T  =0  for  x < x , 

w °°  s 

^W-Ta>  = ^or  x > xs»  anc*  Power  law  wall  temperature  distributions  of 

the  form  T (x)  - T = Ax11.  In  order  to  Isolate  the  effects  of  the  para- 

w oo 

meters  n and  xg  on  the  boundary  layer  stability,  the  local  wall  temper- 
ature at  a reference  location,  T (x  ,)  - T , is  held  constant  as  n and 

w ret  °° 

xg  are  varied.  In  the  present  work  xre£  “ 5.5  inches,  which  is  the  center 
of  the  region  in  which  stability  characteristics  are  measured.  The  non- 
dimensional  location  of  a step  change  in  wall  temperature  is  specified  by 
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the  parameter  S = x /x 

s ref 

A similarity  solution  of  the  boundary  layer  momentum  and  energy 
equations  for  power  law  wall  temperature  distributions,  valid  for  fluid 
property  variations,  is  obtained  using  a modified  Howarth-Dorodnitsyn 
transformation.  The  resulting  ordinary  differential  equations  are  inte- 
grated numerically.  An  approximate  solution  of  the  boundary  layer 
equations  for  step  changes  in  wall  temperature,  valid  for  constant  fluid 
properties,  is  obtained  by  assuming  that  the  thermal  boundary  layer 
thickness  is  small  compared  to  the  velocity  boundary  layer  thickness.  A 

1 

closed  form  solution  of  the  resulting  equations  is  obtained  in  terms  of 
incomplete  gamma  functions. 

i 

Mean  temperature  profiles  measured  with  the  hot  film  anemometer 

operating  as  a resistance  thermometer  are  compared  to  predicted  profiles 

for  various  values  of  n and  S in  Figures  1 and  2,  where  H = (T-T^) / (T  -T^) 

and  n = yAje/vx.  The  case  S = 0 corresponds  to  uniform  wall  heating 

beginning  at  the  leading  edge,  while  S = 1 corresponds  to  a step  change 

in  temperature  at  x , = 5.5  inches.  The  case  n = -.5  corresponds  to 
ref 

input  of  a finite  amount  of  heat  at  the  leading  edge  followed  by  an 
adiabatic  wall.  For  free  stream  velocity  distributions  of  the  form 
* C^"1,  case  occurs  whenever  n = -(nri-l)/2.  For  values  of 

n < -(nrt-l)/2  the  direction  of  heat  transfer  changes  sign,  as  shown  in 
Figure  1 for  n = -0.6. 

Both  analytic  and  experimental  results  Indicate  that  the  mean 


velocity  profile  is  virtually  independent  of  variations  in  n and  S over 
the  range  of  wall  temperatures  used  in  the  present  work.  For  values 
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of  n in  the  range  -0.5  < n < 1.0,  the  wall  shear  stress  and  profile  form 

factor  vary  by  less  than  1.5%  at  a fixed  wall  temperature  level. 

Growth  rate  characteristics  of  ribbon-generated  disturbances  in 

the  boundary  layer  are  determined  in  the  following  manner.  The  Reynolds 

number  R_.  = U 5*/v  is  first  determined  using  the  displacement  thickness, 
o*  e 

6*,  of  the  mean  velocity  profile  measured  at  xref  At  a fixed  Reynolds 
number  and  ribbon  frequency  the  disturbance  amplitude  is  then  measured 
at  five  stations  between  x * 5 in.  and  x = 6 in.  The  amplitude  recorded 

at  each  station  is  the  peak  amplitude,  defined  as  A(x)  = [u'(n,x)/U  ] , < 

e max 

found  by  searching  through  the  boundary  layer  in  the  y-direction.  The 

1 dA  6* 

spatial  disturbance  growth  rate,  -a  = (—  — ) — — , is  then  calculated 

1 A ref  <$* 

using  the  slope  of  a polynomial  curve  fit  of  the  A(x)  data.  By  repeating 
the  above  process  for  several  different  frequencies  the  disturbance 
growth  rate  vs.  frequency  characteristics  are  determined  for  a fixed 
Reynolds  number  and  temperature  distribution. 

Disturbance  growth  rate  characteristics  for  uniform  wall  temper- 
ature distributions,  n = 0,  at  Rr.  = 800  are  shown  in  Figure  3.  The 

ox 

solid  lines  are  smooth  curves  drawn  through  the  measured  points,  is 

defined  as  w = 2irfv/U  ^ , and  f is  the  ribbon  frequency.  These  results 
r e 

provide  a measure  of  the  effectiveness  of  uniform  wall  heating  in  the 

reduction  of  disturbance  growth  rates  at  a fixed  Reynolds  number. 

Stability  characteristics  for  varying  values  of  n,  with  R = 800, 

o* 

Tw(xfef)  - Tm  ■ 5°F  are  shown  in  Figure  4.  Experimental  data  points  are 
not  shown  for  the  sake  of  clarity,  while  unheated  results  have  been 
included  as  a reference.  For  n < 0 the  maximum  disturbance  growth  rates 


«• 


A 


-143- 


and  the  band  of  frequencies  which  receive  amplification  are  both  reduced 
relative  to  the  uniform  wall  temperature  case.  For  n = 1 the  maximum 
disturbance  growth  rate  is  greater  than  that  for  n = 0,  and  the  band  of 
amplified  frequencies  has  increased  relative  to  the  uniform  wall  temper- 
ature case.  The  results  shown  here  are  typical  of  those  obtained  at 

T (x  ,)  - T = 3°F  and  8°F  as  well.  The  behavior  of  the  maximum  growth 
w ref  00 

rate  as  a function  of  n and  T (x  _)  - T is  shown  in  Figure  5. 

w ref  °° 

The  effect  of  varying  S at  R_.  = 800,  with  T (x  c)  - T * 3°F, 

o*  w ret  °° 

is  shown  in  Figure  S.  The  results  shown  here  are  similar  to  those  obtained 

at  T (x  ,)  - T = 5°F  as  well.  The  maximum  disturbance  growth  rate 

displays  a minimum  between  S = 0 and  S = 1 for  each  temperature  level, 

as  shown  in  Figure  7.  This  minimum  occurs  near  the  minimum  critical 

Reynolds  number  of  the  unheated  boundary  layer.  The  measured  value  of 

(R_.)  is  400,  which  corresponds  to  S = .25,  while  the  predicted 
0*  m.c. 

AT=0 

parallel-flow  value  of  (R..)  = 520  corresponds  to  S = .42. 

6*  m.c. 

AT=0 


«• 
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Fit.  A DISTURBANCE  GROWTH  CHARACTERISTICS  FOR  POWER  LAW 
VARIATIONS  IN  WALL  TEMPERATURE.  RJ#  • 800 


Fl«.  6 DISTURBANCE  GROWTH  CHARACTER  I SITCS  FOR  STEP  CHANGES 
IN  WALL  TEMPERATURE,  R,.  * 000 


-146- 


Section  V 

TRANSITION  IN  AN  AXISYMMETRIC  BOUNDARY  LAYER 
WITH  WALL  HEATING 
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TRANSI.T10N  IN  AN  AXISYMMETRIC  BOUNDARY  I.AYER 
WITH  WALL  HEATING 


Recent  numerical  calculations  of  boundary  layer  stability  have 
suggested  the  possibility  of  obtaining  transition  at  very  high  Reynolds 
numbers  in  water  by  using  wall  heating.  The  flow-tube  experiment  is 
intended  to  check  these  predictions,  and  to  determine  the  feasibility 
of  using  heat  combined  with  favorable  pressure  gradient  to  delay 
transition  on  submerged  bodies. 

The  flow  tube  is  basically  an  inside-out  analog  of  the  boundary 
layer  over  a submerged  body  of  circular  cross-section  moving  through 
the  water.  A free-stream  flow  with  very  low  disturbance  levels  is 
provided  through  a 24-inch  diameter,  12-foot  long  settling  chamber 
(see  figure  1).  This  flow  then  passes  through  a 34.5  to  1 smooth 
contraction  into  a 4-inch  diameter  circular  tube.  The  boundary  layer 
in  the  4-inch  tube  remains  thin  relative  to  the  tube  diameter  throughout 
the  length  of  the  tube  (6*  jc  0.05  D)  . 

Heat  is  applied  to  the  inside  wall  of  the  flow  tube  by  means  of 
electric  strip  heaters  wrapped  around  the  outside.  These  heaters  are 
servo-controlled  using  thermocouples  near  the  inside  wall  of  the  tube. 

Thus  the  inside  wall  temperature  is  maintained  at  pre-selected  values. 

The  length  of  the  tube  is  divided  into  20  sections  whose  wall  temperatures 
are  controlled  independently.  Thermocouples  are  also  located  on  the 
outside  wall  of  the  tube  at  about  thirty  locations  along  its  length. 

The  temperature  difference  between  the  inside  and  outside  wall  at  a 
given  location  provides  a measurement  of  the  heat  flux  into  the  boundary 
layer.  Since  this  heat  flux  will  increase  by  a factor  of  about  ten  at 
the  region  of  transition,  it  is  easy  to  locate  transition  with  temperature 
measurement.  There  is  no  instrumentation  penetrating  the  inside  wall 
of  the  tube  which  could  trip  the  boundary  layer.  The  wall  roughness  is 
less  than  four  micro-inches  and  the  wavincss  is  less  than  one  mil  per  inch. 
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Three  different  flow  tube  test  sections  are  being  used.  The  first 
is  20  feet  in  length  and  uniform  in  diameter.  This  tube  will  simulate 
a nearly  zero  pressure  gradient  boundary  layer.  (Displacement  effect 
yields  a very  weak  favorable  gradient.)  At  the  maximum  flow  velocity  of 
90  feet  per  second,  the  length  Reynolds  number  of  this  tube  is  125 
million.  The  second  and  third  flow  tube  will  be  of  non-uniform  diameter 
to  simulate  two  axisymmetric  shapes  of  practical  interest.  The  longest 
of  these  tubes  will  be  44  feet  in  length. 

The  experiment  is  located  at  the  Hydro-Mechanics  Laboratory  of 
Colorado  State  University.  This  site  was  chosen  over  several  ether 
possibilities  for  the  following  reasons.  Although  the  water  is  not 
recirculated,  the  run  time  is  effectively  unlimited.  The  lab  draws 
its  water  directly  from  a large  reservoir  through  a 24-inch  diameter 
pipe.  The  total  pressure  at  the  entrance  of  the  settling  chamber  is 
about  ninety  psig.  Thus  even  at  a test  section  velocity  of  90  feet  per 
second,  the  test  section  pressure  is  about  35  psig  and  the  possibility 
of  cavitation  can  be  ignored.  The  longer  run  time  and  higher  total 
pressure  are  the  chief  advantages  of  this  facility  relative  to  others. 

The  first  series  of  measurements  in  the  experiment  was  a careful 
documentation  of  the  flow  in  the  settling  chamber  and  test  section 
without  wall  heating.  The  test  section  in  this  case  was  a 4-foot  straight 
tube  with  several  types  of  instrumentation.  A pitot-rake  was  located  near 
the  downstream  end  to  check  the  uniformity  of  the  free-stream  mean 
velocity.  A cylindrical  hot  film  sensor  was  used  to  measure  freestream 
turbulence  intensity  and  to  determine  the  intermittency  in  the  boundary 
layer.  This  probe  was  used  to  measure  the  transition  Reynolds  number 
in  the  absence  of  heat,  which  will  be  discussed  below.  Two  flush-mount 
wall  pressure  transducers  were  also  used  to  locate  transition  in  the 
tube.  An  accelerometer  was  mounted  on  the  outside  wall  to  measure  the 
vibration  spectrum. 

A hot  film  probe  was  also  used  to  survey  the  flow  in  the  settling 


l 
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chamber  one  foot  upstream  of  the  contraction.  Mean  velocity  profiles, 
turbulence  levels,  and  turbulence  spectra  were  measured  for  various 
configurations  of  turbulence  management  devices.  The  devices  used 
include  1/8  inch  cell  aluminum  honeycomb,  10  pores  per  inch  porous  foam, 
and  wire  screen  of  34  and  60  mesh  per  inch.  For  each  turbulence 
management  configuration,  the  settling  chamber  flow  was  documented  at 
several  test  section  velocities.  In  addition,  the  unheated  transition 
Reynolds  number  was  measured  in  the  flow  tube  at  1 flow  velocity.  Since 
this  transition  Reynolds  number  has  been  measured  in  other  facilities, 
it  is  a good  indication  of  the  quality  of  the  freestreara  flow.  A plot 
of  unheated  transition  Reynolds  number  versus  settling  chamber  configu- 
ration is  shown  in  figure  2.  The  "best"  configurations  are  seen  to 
produce  transition  at  Reynolds  numbers  close  to  5 million.  This  is  very 
close  to  the  highest  value  achieved  in  the  air  flow  tube  facility  of 
C.  S.  Wells. 

Hot  film  surveys  show  that  the  boundary  layer  in  the  settling 
chamber  becomes  turbulent  at  test  section  velocities  above  30  feet 
per  second.  Since  the  boundary  layer  at  the  downstream  end  of  the 
contraction  must  be  thin  and  disturbance-free,  a boundary  layer  suction 
section  was  installed  between  the  settling  chamber  and  the  contraction. 

This  section  uses  a 4-inch  length  of  porous  wall  to  remove  the  entire 
turbulent  boundary  layer.  Measurements  have  shown  that  even  at  the 
maximum  test  section  velocity,  the  boundary  layer  at  the  entrance  of  the 
contraction  can  be  made  laminar  with  the  appropriate  amount  of  suction. 

The  turbulence  levels  in  the  contraction  entrance  are  considerably 
lower  than  the  original  goal  of  about  1 per  cent.  At  a test  section 
velocity  of  45  feet  per  second,  a turbulence  level  in  the  settling  chamber 
of  0.2  per  cent  is  typical  for  the  best  configurations.  This  configuration 
makes  use  of  1-inch  of  porous  foam  followed  by  two  sections  of  honeycomb 
and  then  three  screens.  Preliminary  results  using  a 3-inch  thickness  of 
porous  foam  show  even  lower  turbulence  levels. 


A 
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Some  tests  with  wall  heating  were  made  on  the  4-foot  flow  tube, 
The  results  of  these  were  not  strongly  positive  - only  a small  increase 
in  transition  Reynolds  number  was  seen.  However,  the  heating  and 
temperature  measuring  system  on  the  tube  were  very  crude.  Heat  was  not 
applied  far  enough  upstream  in  the  tube  and  the  suction  section  was  not 
yet  available.  At  the  date  of  this  writing,  the  20-foot  straight  flow 
tube  is  being  installed.  Early  results  with  this  tube  may  be  available 
at  the  time  of  this  meeting. 
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Section  V 

PRESSURE  FLUCTUATIONS  IN  THE  TRANSITION  REGIONS 
OF  FOREBODIES  OF  REVOLUTION 
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PRESSURE  FLUCTUATIONS  IN  THE  TRANSITION  REGIONS 
OF  FOREBODIES  OF  REVOLUTION 

by 


Thomas  T.  Huang 

David  W.  Taylor  Naval  Ship  Research  and  Development  Center 

Bethesda,  Maryland 


I 


The  transition  resulting  from  the  growth  of  Tollmein-Schlichting 
disturbances  naturally  present  in  a laminar  boundary  layer  is  defined 
as  natural.  Natural  transition  should  be  distinguished  from  transition 
triggered  by  laminar  separation.  If  laminar  separation  occurs,  either 
as  a result  of  a severe  adverse  pressure  gradient  due  to  body  shape  or 
by  the  use  of  an  artificial  turbulence  stimulator,  transition  will 
usually  occur  at  a very  short  distance  downstream  of  the  flow  reattach- 
ment , 


In  the  present  investigation,  the  amplitude  of  the  wall-pressure 
fluctuations  during  natural  transition  and  following  laminar  separation 
due  both  to  adverse  pressure  gradient  and  to  two-dimensional  trip  wires 
was  measured  by  microphones  located  in  the  transition  regions;  these 
measurements  were  extended  downstream  into  the  fully-established 
turbulent  boundary  layer.  The  fluctuating  pressure  on  the  model  sur- 
face was  sensed  through  an  0. 794-mm-diameter  pinhole  which  led  to  a 
small  cavity  enclosed  by  the  protective  cap  and  the  microphone  diaphragm. 
The  microphone  system  used  had  satisfactory  response  characteristics 
between  10  Hz  and  10  KHz. 

The  wind  tunnel  in  the  Anechoic  Flow  Facility  of  the  Center  has  a 
2.438  meters  x 2.438  meters  test  section.  The  free-stream  rms  turbu- 
lence levels  in  the  test  section  with  a model  in  place  were  found  to  be 
0.075,  0.090,  0.100,  and  0.12-0.15  percent  for  U0  = 24.4,  30.5,  30.8 
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and  45.7  tneters/sec,  respectively.  By  integrating  the  measured  noise 
spectrum  levels  in  the  test  section  from  0 to  10,000  Hz  to  background 
acoustic  noise  at  30.5  meters/sec  was  found  to  be  93  d!3:re  0.0002 
dyn/cm^,  which  corresponds  to  an  equivalent  rms  turbulence  level  of 
0.007  percent. 

The  model,  which  consisted  of  the  forebody,  together  with  a 
parallel  middle  body  and  a streamlined  afterbody,  was  supported  by  a 
single  streamlined  strut  located  roughly  at  the  middle  of  the  model. 

Three  small  guide  struts  were  used  at  the  tail  of  the  afterbody  to 
increase  the  rigidity  of  the  model. 

Two  series  of  forebodies  of  revolution  were  selected  for  this  inves- 
tigation. One  of  the  series,  designated  "Forebody  T,"  had  a natural 
transition  region,  and  the  other  series,  designated  "Forebody  S,"  had 
laminar  separation  which  triggered  transition.  Two-dimensional  rough- 
nesses (trip  wires),  located  at  various  positions  upstream  of  natural 
transition  and  laminar  separation,  were  also  investigated.  Figure  1 
shows  the  geometries  of  the  forebodies  used  and  the  computed  pressure 
coefficient.  Figure  2 shows  typical  histogram  of  microphone  output. 

The  growth  of  Tollmein-Schlichting  disturbances  in  boundary  layers 
is  approximated  by  the  spatial  stability  charts  for  the  Falkner-Skan 
boundary  layer  profile  prepared  by  Wazzan,  Okamura,  and  Smith.  Smith's 
spatial  amplification  ratio  was  chosen  as  a simple  yardstick  to  corre- 
late the  instability  characteristics  of  the  boundary  layer  with  various 
measured  stages  of  transition.  Whereas  the  amplification  ratio  is  calcu- 
lated from  linear  stability  theory,  the  latter  stage  of  transition  is  a 
nonlinear  process  and  is  three-dimensional  in  nature.  The  calculated 
spatial  amplification  ratios  A(x)  in  the  natural  transition  regions  of 
Forebodies  T-3,  T-6  and  T-8  correlated  with  the  flow  regimes  measured 
by  flush-mounted  hot-films  are  shown  in  Figures  8,  10  and  11  in  the 
paper  by  Power*.  The  measured  flow  regimes  of  T-8  do  not  agree  with 
calculated  amplification  ratios.  The  transition  characteristics  mea- 
sured in  the  present  wind  tunnel  do  agree  well  with  towing  tank  results. 

Figure  4 shows  the  correlation  between  the  measured  wall -pressure 
fluctuations  and  calculated  Smith  spatial  amplification  ratios  A(x)  in 
the  natural  transition  regions  of  Forebodies  T-3  and  T-6.  The  pressure 
fluctuations  are  small  for  A < e$  to  the  position  of  el5  where  the  flow 
has  just  become  fully  turbulent  according  to  the  outputs  of  hot-film 
probes.  The  measured  wall -pressure  fluctuations  approach  those  of  a 
ful ly-cstabl ished  turbulent  boundary  layer  at  approximately  e^®.  How- 
ever, no  large  wall-pressure  fluctuations  in  the  transition  regions  of 
Forebodies  T-8  was  measured. 

♦ 

Power,  J.L.,  "A  Comparison  Between  Measured  and  Computed  Locations 
of  Transition  on  Nine  Forebodics  of  Revolution,"  paper  presented  at 
Rand's  Workshop  II,  1976 
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The  characteristic:  lengths  K and  X,  which  is  the  dead-air  bubble 
length,  are  shown  in  Figure  5.  The  data  of  Caster  and  Arakeri  for  X 
are  also  included  in  Figure  5.  As  shown  in  Figure  6,  enormous  pressure 
fluctuations  were  measured  in  The  region?,  of  flow  reattachment  following 
laminar  separation  on  Forebodies  S-2,  S-l,  and  S-7.  If  the  distance  is 
measured  from  the  separation  point,  S-S^,  and  is  nondimcnsiorialized  by 
a characteristic  length  5,  the  distance  between  the  point  of  maximum 
wal 1-prcssure. fluctuation  and  the  separation  point,  then  the  measured 
value  of  V P,2/(pU02/2)  possess  a similar  form  with  respect  to  the 
quantity  (S-Ss)/5.  The  velocity  fluctuations  in  the  shear  layer  right 
above  the  end  of  the  dead-air  region  were  found  by  Gaster  to  be  extremely 
high  (,fu'-7uo  is  about  0.16).  Consequently,  the  pressure  fluctuations 
are  expected  to  be  extremely  high  there,  although  direct  measurements 
have  not  been  made. 

The  turbulent  stimulators  used  were  single  circular  wires  of  four 
different  diameters:  -0.05,  0.127,  0.254,  and  0.508  mm  - located  at 
different  axial  positions  on  the  Forebodies  T-3  and  S-2.  Figure  7 shows 
the  locations  of  inception  of  intermittent  disturbances  measured  by  hot- 
film  probes  and  maximum  wall  pressure  fluctuations,  respectively  c/d 
and  5/d,  for  four  wires  of  large  R^,  and  the  apparent  amplification 
factor  in  terms  of  [AlnA]j^  for  wires  having  Rj(  in  the  range  25  < Rj^  < 100. 
Here,  Rg  = Uj/K/v  is  the  wire  Reynolds  number  based  on  the  velocity  uy 
at  the  tip  of  roughness  height  K.  The  values  of  [AlnA]}(  were  obtained 
by  correlating  the  measured  wall-pressure  fluctuations  with  and  without 
trip  wires.  Preston's  minimum  Rg  of  400  to  600  is  a satisfactory 
criterion  to  determine  the  range  of  t Lp  wire  Re's  which  are  effective 
in  causing  rapid  transition.  No  trace  of  high  pressure  fluctuations 
due  to  the  trip  wire  was  measured  after  200  wire  diameter  downstream  of 
the  wire.  The  Smith  and  Clutter  criterion  of  a maximum  value  of  of  25 
so  as  not  to  prematurely  trip  transition  agrees  satisfactorily  with  the 
present  results.  For  the  Rg  range  25  < RK  < 100,  transition  will  move 
upstream. 

Cavitation  inception  of  naval  vehicles  will  be  controlled  by  the 
unsteady  pressure  fluctuations  in  regions  of  natural  transition  or  laminar 
separation,  superimposed  upon  the  potential-flow  static  pressure.  The 
location  and  magnitude  of  the  high  pressure  fluctuations  in  the  transition 
regions  have  profound  implications  relative  to  model/full-scale  cavitation- 
inception  scaling,  boundary- layer  self-noise  scaling,  and  other  naval 
applications  related  to  transition  phenomenon. 

ACKNOWLEDGEMENTS 


This  work  was  funded  by  the  David  W.  Taylor  Naval  Ship  Research  and 
Development  Center  (DTNSRDC)  under  its  Independent  Research  Program, 

Task  2R023  0101,  Work  Unit  1-1552-102. 


_0.2 


_0.4 


_0.6 


_0.8 


Figure  I- a Surface  Profiles 
Three  Forebodies 
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Typical  Histogram  of  Microphone  Output  in  Natural  Transition  Region 
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1^- 

10 

in 

ro 

0J 

o 

o 

o 

o 

O 

o 

6 

o 

6 

d 

o 

d 

a 


Amplification 


IO 

O 

w 

o 

m 

o 

CM 

CO 

8-01  x °y 

1 

• — ^ . . . — 

Model  Reynolds  Number  R, 


-165- 


Figure  3-b  Measured  Flow  Regimes  on  Forebody  T-6  as  a 
Function  of  Reynolds  Number  RD 
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Figure  3-C  Measured  Flow  Regimes  on  Forebody  T-8  as  a 
Function  of  Reynolds  Number  Rn 
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Figure  5 Measured  Separation  Bubble  Lengths  Versus  Reynolds  Number 
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A COMPARISON  BETWEEN  MEASURED  AND  COMPUTED  LOCATIONS  OF  TRANSITION 
ON  NINE  FOREBODIES  OF  REVOLUTION 
by 

John  L.  Power 

Drag  and  transition  data  are  given  for  nine  submerged  axisymmetric 
bodies.  The  bodies  have  vastly  different  forebody  geometries  ranging 
from  very  fine  to  very  blunt  profiles  and  include  two  flat-faced  forebodies. 
The  experiments  were  performed  in  a towing  basin  at  speeds  between  0.5  and 
5.1  m/s  and  at  a water  depth  of  2.7  metres.  The  experimental  arrangement 
is  shown  in  Figure  1.  The  models  were  approximately  7 metres  long  and  had  a 
maximum  diameter  of  0.62  metres.  Constant  volume  was  maintained  by  varing 
the  length  of  the  parallel  middle  body. 

The  geometries  of  the  forebodies  were  altered  by  changing  two  parameters, 
the  length/maximum  diameter  ratio  and  the  prismatic  coefficient.  The 
profiles  of  four  of  the  nine  forebodies  are  shown  in  Figure  2.  These 
forebodies  have  length -diameter  ratios  of  0.5,  1.00,  1.82  and  3.00  with  a 
constant  prismatic  coefficient  of  0.667.  The  computed  pressure  distributions 
on  these  forebodies  are  shown  in  Figure  3. 

The  contours  of  a second  set  of  five  forebodies  are  shown  in  Figure  4. 
These  shapes  are  much  blunter  than  those  shown  in  Figure  2 and  include  two 
flat -faced  bows.  The  increased  bluntness  was  obtained  by  increasing  the 
value  of  the  prismatic  coefficient.  The  computed  pressure  distributions 
on  these  bodies  are  shown  in  Figure  5. 

The  extents  of  the  laminar,  transitional  and  turbulent  flow  regimes 
were  determined  by  qualitative  analysis  of  the  signals  from  several  constant- 
temperature  hot  films  flush  mounted  on  each  forebody.  Selection  of  the  hot 
film  locations  was  assisted  by  using  transition  calculations,  computer 
pressure  distributions  and  the  Curle-Skan  laminar  separation  criterion. 

The  active  elements  of  the  hot  films  were  approximately  1.6  millimeters 
long.  They  were  mounted  on  the  end  of  small  cylinders  which  were  glued 
into  pre-drilled  holes  in  the  models. 

As  shown  in  Figure  6 analysis  of  the  hot  film  probe  outputs  revealed 
four  types  of  boundary  layer  flow:  (1)  laminar  flow,  where  the  hot-film 
signals  were  of  steady  amplitude;  (2)  a smooth  wave-like  disturbance; 

(3)  an  intermittent  turbulent  bursting  flow;  and  (4)  a full v- turbulent 
flow  of  random  nature. 
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Figure  7 gives  a summary  of  the  forebody  properties  and  locations  of 
transition  or  laminar  separation  for  the  nine  models.  The  locations  given  « 

for  laminar  separation  are  those  computed  using  the  Curle-Skan  modified 
Thwaites  criterion.  The  transition  locations  were  determined  from  the 
analysis  of  the  hot  film  records. 

i 

Laminar  separation  was  predicted  for  five  of  the  nine  models.  Hot 
film  probes  were  placed  immediately  before  and  after  these  predicted  ' 

locations.  In  all  cases,  the  probe  before  the  separation  location  predicted 
laminar  flow  (Type  1)  and  the  probe  after  the  separation  location  indicated 
turbulent  flow  (Type  4) . It  was  concluded  that  laminar-separation  initiated 
turbulence  and  the  location  of  the  laminar  separation  is  closely  predicted 
by  the  modified  Thwaites  criterion.  i 


The  flow  regimes  for  the  four  forebodies  on  which  laminar  separation 
did  not  occur  are  given  in  Figures  8,  9,  10,  and  11.  The  experimental  data 
are  compared  with  predicted  spatial  amplification  ratios  of  disturbances  as 
calculated  by  the  Smith-Gamberoni  method.  In  Figures  8 and  9,  three 
predictions  corresponding  to  spatial  amplification  ratios  of  e , e^,  and 
e!3  are  compared  with  the  measured  values  from  models  4620-3  and  4.  For 
both  models,  the  inception  of  intermittent  turbulent  bursting  of  transition 
and  the  e^  curves  agree  quite  well.  The  e9  curves  appear  to  coincide  with 
the  measured  inception  of  Type  2 wave  like  disturbances.  The  e^  curves 
correlate  very  well  with  the  measured  inception  of  fully-turbulent  flows. 

In  Figure  10  five  predictions  corresponding  to  spatial  amplification 
ratios  of  e^,  e7,  e9,  e^*,  and  e^  are  compared  with  the  measured  results 
from  model  4620-6.  For  this  model  the  e9  curve  is  located  close  to  the 
inception  of  intermittent  turbulent  bursting  of  transition  while  the  ell 
curve  is  located  well  into  the  turbulent  bursting  region.  For  the  lower 
Reynolds  numbers,  the  e*3  curve  correlates  with  locations  near  the  end  of 
the  turbulent -bursting  region  and  at  the  higher  Reynolds  numbers,  this  curve 
moves  well  into  the  fully  turbulent  (type  4)  region.  The  e5  curve  correlates 
well  with  the  inception  of  type  2 wave  like  disturbances. 


In  Figure  11,  predictions  corresponding  to  spatial  amplification  ratios 
of  e , e^,  e7,  e^,  e**,  and  e*3  are  compared  with  measured  results  from  model 
4620-8  (a  flat-faced  model).  The  correlation  between  these  predictions  and 
the  measured  results  is  radically  different  from  the  correlations  obtained 
with  the  three  models  previously  discussed.  Here  the  e^  curve  correlates 
with  the  inception  of  wave-like  disturbances.  The  e^  curve  is  located  at 
the  end  of  the  wave-like  disturbance  region  and  near  the  start  of  the 
turbulent  bursting  region.  The  start  of  fully  turbulent  flow  occurs 
between  the  e^  and  e7  curves.  The  e7  and  all  higher  exponent  curves  lie 
well  within  the  fully  turbulent  region. 


Also  shown  in  Figure  11  are  transition  data  obtained  by  T.  Huang 
in  a.  low-turbulence  wind  tunnel  using  Model  4620-8.  The  wind  tunnel  data 
are  in  substantial  agreement  with  the  towing  tank  data.  Huang  also  measured 


L 


-175- 


the  transition  flow  regimes  on  Models  4620-3  and  8.  These  data  are  not 
shown  but  similar  agreement  between  towing  tank  and  wind  tunnel  results 
were  found  for  these  models. 

Figure  12  gives  the  results  obtained  from  the  drag  measurements.  For 
each  model,  the  computed  residual  drag  coefficients  and  EHP  ratios  are  given. 
Model  4620-3  and  6 show  the  lowest  drags  and  Model  4620-7  shows  the 
highest  drag.  There  is  a difference  of  approximately  4 percent  in  EHP 
required  for  these  two  extremes. 


Ref.  1:  Huang,  T.T.,  "Pressure  Fluctuations  in  the  Transition  Regions  of 
Forebodies  of  Revolution."  Low-Speed  Boundary- Layer  Transition  Workshop:  II, 
Rand  Corp.,  Santa  Monica,  California  (Sept.  1976). 


Figure  2— Profiles  of  Four  Forebodies  of  Series  I 
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Figure  3 [Computed  Distributions  of  Pressure 
LCoef  ficient  on  the  Axisymmetric 
Forebodies. 
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Figure  7 

Forebody  Properties  and  Location  of  Transition 
or  Laminar  Separation 


Model 

VD  i' 

CP£ 

Location  of  Transition 
or  Laminar  Separation 

1 

0.500 

0.667 

L.  S., 

x£/D  = 0.472 

2 

1.000 

0.667 

L • S.  | 

xt/D  = 0.894 

3 

1.820 

0.667 

Trans . , 

1.47<  xjp/D<1.60 

4 

3.000 

0.667 

Trans . , 

1.07<  x /D<3.30 

5 

1.000 

0.850 

L • S*| 

x£/D  = 0.371 

6 

1.820 

0.850 

Trans . , 

0.57<  x /D<1 .04 

Xr 

7 

0.500 

0.850 

L.  S., 

Xj/D  = 0.184 

6* 

1.216 

0.823 

Trnas . , 

0.07<  xa/D<0.70 

9* 

1.000 

0.933 

L.  S., 

x^/D  = 0.202 

D « 0.624  m,  Volume/ D3  = 6.786 

D - model  maximum  diameter 

- location  of  transition  or  laminar  separation 
* . flat  faced  model 
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Figure  8 - Mcosured  Flow  Regimes  on  Smooth  Figure  9 - Mcosurcd  Flow  Regimes  on  Smooth 

l-E/D  = l.82  Forebody  >-E/0«3.0  Forobody 
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Figure  12 

Values  of  Residual  Resistance  Coefficient 
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VI-1  THE  COMBINED  EFFECTS  OF  PRESSURE  GRADIENT  AND  HEATING 
ON  BOUNDARY-LAYER  STABILITY  AND  TRANSITION^ 

it 

by  A.  R.  Wazzan  and  Carl  Gazley,  Jr. 

The  Rand  Corporation 
Santa  Monica,  California 

It  is  well  known  that  major  parameters  affecting  transition  from 
laminar  to  turbulent  flow  include  pressure  gradient  and  heat  transfer. 

It  is  also  commonly  accepted  that  some  insight  into  the  effect  of  these 
parameters  on  transition  can  be  obtained,  at  least  in  the  case  of  boundary 
layer  flows,  from  a study  of  the  effect  of  these  parameters  on  the  stabil- 
ity of  boundary  layers  to  small  disturbances  (Tollmien-Schlichting  waves). 
Previous  calculations^  on  the  effect  of  heat  transfer  on  the  stability 
characteristics,  and  in  particular  the  critical  Reynolds  number  R of 

laminar  water  boundary  layers  show  that  surface  heating  stabilizes  the 
flow  whereas  cooling  destabilizes  it.  In  flat  plate  flow,  e.g.,  with 
Tfi  = 60°F,  Rcr^t  increases  by  about  thirty  fold  for  a wall  temperature  of 
T^  — 135°F.  In  stagnation  flow  (Hartree  8 = 1.0)  complete  stabilization  is 
obtained  with  increasing  T^  whereas  in  flat-plate  flow  (8  = 0)  and  separating 

boundary  layers  (8  = -0.1988),  R . exhibits  a maximum  with  increasing 

crit 

surface  temperature.  These  results  show  that  heating  has  a strong  effect 
on  stability  and  that  this  stabilizing  (destabilizing)  effect  is  strongly 
dependent  on  the  interaction  between  the  pressure  gradient  and  the  surface 
temperature  distribution.  To  ascertain  the  nature  and  extent  of  this 
interaction,  further  stability  calculations,  at  other  values  of  Hartree  8, 
are  desirable.  In  addition,  such  detailed  calculations  would  he  useful 
in  assessing  the  effect  of  heating  on  the  stability  and  transition  over 
bodies  of  practical  interest  (e.g.,  bodies  of  revolution  where  the  local 
pressure  gradient  and  boundary  layer  can  be  approximated  by  an  appropriate 
wedge  flow).  For  these  reasons,  we  computed  the  spatial  stability  (Rcrit. 
spatial  amplification  rates,  etc.)  for  the  following  two-dimensional 
wedge  flows:  8 - -0.15,  -0.10,  -.05,  0,  0.10,  0.20,  0.30,  and  0.40  (at 
T - T - 0,  5,  20,  and  30°F)  for  T - 67°F. 

W G G 

- 

University  of  California,  Los  Angeles;  Consultant,  the  Rand  Corporation. 

^Papcr  presented  at  the  Low-Speed  Boundary-Layer  Transition  Workshop  II, 
The  Rand  Corporation,  Santa  Monica.  SeLptembcr  13-15,1976. 
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The  variation  of  the;  following  boundary-layer  characteristics, 

-I/O 

0*/x/Uc  , O/x/Rc  , C,/2/Ro  , and  Nu/Rc  Pr  with  (T  - T ) 
x xrx  xe  we 


are  shown  in  Fig.  1.  R. 


with  (T  - T ) is  shown  in  l'ig.  2.  The 


effect  of  surface  heating  on  boundary-layer  stability  is  due  primarily  to 
the  variation  of  viscosity  with  temperature.  The  results  are  thus  depend- 
ent no*’,  only  on  the  temperature  difference  (T  - T ) but  also  on  the 

we 

temperature  level,  Tg.  This  is  evident  in  Fig.  2 where  previous  results 
for  Tg  = 60  F are  shown  in  comparison  with  the  present  results  for  T^  = 
67°F.  Even  this  relatively  slight  change  in  ambient  temperature  results 
in  an  appreciable  difference  in  the  rate  of  change  of  viscosity  with 
temperature,  and  consequently  on  the  predicted  stability. 

Rayleigh  theorems  for  inviscid  instability  state  that  for  a boundary- 
layer  flow,  the  necessary  and  sufficient  condition  for  amplified  and 
neutral,  inviscid  instability  is  that  U"(y)  must  vanish  somewhere  in  the 
boundary  layer,  i.e.,  the  mean  velocity  profile  must  be  inflected.  In 
addition,  some  correlation  appears  to  exist  between  R and  the  location 


of  the  inflexion  point. 


In  fact,  in  adiabatic  flows,  T T , R 


decreases  as  the  location  of  the  inflexion  point  moves  away  from  the  wall. 

That  is,  in  adiabatic  flows,  R . decreases  as  8 becomes  more  negative. 

crit 

In  flows  with  favorable  pressure  gradient  (8  > 0)  where  the  profiles  are 

not  inflected,  the  boundary-layer  characteristic  of  great  importance  to 

(2) 

its  stability  characteristics  is  U"(y)  and  to  a lesser  extent  U"(0). 

In  adiabatic  flows  it  is  known  that  for  8 > 0,  where  U"  (0)  > 0,  Rcrit 

for  8 < 0 is  smaller  than  R ..  for  8 > 0.  In  fact,  when  a U" (y)  distribu- 

tion  is  the  result  of  pressure  gradient  effects  only  (adiabatic  flows) 

a strong  correlation  exists  between  R and  8 or  U"(y)  or  simply  U"(0) 

(2  3)  crit 

or  the  shape  parameter  H.  ’ Rcrit  exhikits  a simple  variation  with 
U"(C)  or  with  H;  R t increases  as  H decreases  (Fig.  3)  or  as  U"(0) 
decreases  and/or  becomes  increasingly  more  negative.  This  dependence  on 
U"  or  11,  in  fact,  is  to  be  expected.  An  inspection  of  the  Orr-Sonmerfeld 
equation  shows  the  boundary-layer  characteristics  that  directly  influence 

the  eigen  valyes,  and  hence  R , arc  U"(y)  and  U(y)  with  U"(y)  being 
the  dominant  term.  Therefore,  it  may  be  assumed  that  the  heating  of  water 
boundary  layers,  which  produces  variation  in  H similar  to  those  produced 


j 


— ISOTHERMAL  WEDGE  FLOWS 

WAZZAM,  OKAMURA,  AND  SMITH  (1968) 

O FLAT  PLATE  WITH  SUCTION  AND  BLOWING 
TSOU  AND  SPARROW  (1970) 

A FLAT  PLATE  WITH  ASYMPTOTIC  SUCTION 
HUGHES  AND  REID  (1965) 
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ig.3 — Critical  Reynolds  number  as  a function  of  the  boundary-layer  shape 
parameter,  isothermal  wedge  flows  and  flat  plate  with  suction  and  blowing 
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through  the  cfl'ect  of  pressure  gradient  alone,  leads  to  increased 

stability  and  particularly  to  increasing  critical  Reynolds  number. 

This  assumption,  however,  is  found  to  be  only  partially  true,  Figs.  4 

and  5.  These  figures  show  that  with  initial  heating  H decreases  and 

R . increases,  in  agreement  with  the  trend  observed  in  Fig.  3.  How- 
crit 

ever.  Figs.  4 and  5 show  that  although  H decreases  monotonies! ly  with 

increasing  surface  temperature,  Rx*  exhibits  a maximum  (at  least 

crit 

for  flows  with  6 < 1.0).  The  difference  in  the  variation  of  Rx* 

crit 

with  H,  when  produced  through  the  effect  of  pressure  gradient  alone 
(adiabatic  flows)  or  through  the  combined  effects  of  pressure  gradient 
and  surface  heating,  can  be  qualitatively  understood  through  an  examina- 
tion of  the  Orr-Sommerf eld  equation. 

In  the  adiabatic  case,  only  U"  and  U appear  in  the  Orr-Sommerf eld 

equation  and  hence  the  monotonic  variation  of  Rx*  with  U"  or  H.  In 

crit 

the  nonadiabatic  case,  the  small  disturbance  equation  is  a modified  Orr- 

Sommerf  eld  equation  that  contains  not  only  U"  and  U but  also  p,  p',  and 

p".  This  alters  the  nature  of  the  problem  in  two  ways.  In  the  adiabatic 

case  the  pressure  gradient  affects  Rx*  mainly  through  the  mean  veloc- 

crit 

ity  or  the  term  U"  (which  can  be  represented  by  some  function  of  H) , 

whereas  in  the  nonadiabatic  case  heating  affects  Rx*  not  only  through 

crit 

the  mean  velocity  or  the  U term  (which  can  still  be  represented  by  some 

function  of  H)  but  also  through  the  terms  p,  p',  and  p"  that  appear  in 

the  new  Orr-Soramerfeld  equation.  Second,  the  nature  of  the  eigen  function 

$ (and  consequently  all  eigen  values  and  properties  depending  on  the  eigen 

values  such  as  R . .)  is  different  in  the  two  cases.  In  the  adiabatic 
crit 

case,  the  differential  equation  for  <J>  (the  Orr-Sommerfeld  equation)  in- 
cludes only  the  function  $ and  even  the  derivatives  <J>"  and  whereas 

in  the  nonadiabatic  case,  the  differential  equation  for  <J>  (the  modified 
Orr-Sommerfeld  equation)  includes  not  only  <j>  and  the  even  derivatives 
and  but  also  the  odd  derivatives  <J>'  and  Therefore,  in  the  case 

of  heating,  although  U"  or  H still  characterize  the  boundary  layer,  they 
do  not  present  a complete  relationship  between  the  mean  flow  and  the 

eigen  function  <J>,  and  hence  Rx*  . Physically,  this  may  be  interpreted 

crit 

as  follows:  with  heating,  initially  II  decreases  rapidly  (indicating  a 

decrease  in  momentum  loss)  (Fig.  6)  and  Rx*  increases.  However,  with 

° crit 
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FJg.4 — Critical  Reynolds  number  for  heated  wedge  flows  in  water 
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Fig.  5 — Previous  results  for  heated  wedge  flows  in  water 
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TW.°F 


Fig.  6 — Variation  of  the  boundary-layer  shape  parameter  H with 
surface  temperature  for  three  wedge  flows,  Te  =60°  F 
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8till  increased  heating,  H continues,  to  decrease  but  at  a much  slower 
rate  (Fig.  6).  In  the  meantime,  p,  p',  and  p"  continue  to  vary  appre- 
ciably with  increasing  T . In  fact  p,  which  has  a destabilizing  effect, 

w 

continues  to  decrease  raonotonically  with  heating  whereas  p'and  p", 

(4) 


(4) 


and  in  particular  p' (0)  and  p"(0),  which  have  a stabilizing  effect, 

reverse  their  trend  in  variation  with  temperature  (change  from  increasing 

with  T to  decreasing  with  T ) near  che  temperature  where  R . _ exhibits 
w w ent 

a maximum  for  wedge  flows  with  B < 1.0.  Therefore,  for  high  heating 

rates  the  variation  of  H with  T becomes  negligible  and  the  variation 

w 


of  p,  p',  p"  with  T dominate  the  effect  of  heating  on  R 

w crit 


at  high  heating  rates  it  is  expected  that  Rg* 


crit 


Therefore, 
will  exhibit  a maximum 


with  T 


In  the  case  of  B = 1.0,  the  maximum  in  R*  with  T is  not  ob- 

°crit  w 

served  because  the  unstable  zone  (region  contained  within  the  neutral 

curve)  is  rather  limited,  and  in  the  initial  stages  of  heating  when  H 

is  fast  decreasing  with  T the  unstable  zone  is  fast  approaching  a point. 

w 

In  fact,  just  at  about  the  temperature  when  the  variation  of  H with  T^ 
begins  to  slow  down  considerably  and  p,  p',  and  p"  begins  to  dominate, 
the  unstable  region  shrinks  to  zero  and  the  flow  becomes  totally  stable. 

In  spite  of  this  entire  discussion  of  the  relative  importance  of  H 
and/or  the  p,  p',  and  p"  terms  on  the  stability  characteristics  of  a 
given  boundary  layer,  Figs.  4 and  5 show  that  over  the  large  range  of 
H values,  a decrease  in  H results  in  increased  stability,  e.g. , increas- 
ing R<$*  j • 

crit 

Of  more  practical  importance  than  the  variation  of  Rcrit  with  sur- 
face temperature  is  the  effect  of  heat  transfer  on  transition  in  boundary- 
layer  flows.  Although  the  transition  process  is  a fairly  complex  one, 
arid  the  relation  between  transition  and  linear  stability  theory  is  at 
best  only  qualitatively  understood,  it  appears  that  in  certain  instances 
linear  instability  theory  is  capable  of  providing  limits  on  the  transition 
Reynolds  number. ^ These  limits  must  be  the  critical  Reynolds  number 
RX  ^ as  a lower  bound,  and  the  Reynolds  number  where  the  most  dangerous 
disturbance  (or  frequency)  attains  a given  amplitude,  say  A,  or  when  the 
disturbance  with  initial  amplitude  Aq,  is  amplified  by  a certain  a = A/A q . 
The  relation  between  transition  and  this  amplification  ratio  may  be 
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* 


established  as  follows.  A disturbance  growing  according  to  linear 
theory,  sooner  or  later  reaches  a state  where  (1)  the  linear  theory 
ceases  to  be  valid,  and  nonlinear  processes  commence;  (2)  the  boundary 
layer  becomes  locally  turbulent — turbulent  spots  are  formed  and  grow  and 
increase  in  number;  and  (3)  these  spots  spread  into  the  neighboring 
laminar  flow  until  the  mean  flow  becomes  fully  turbulent.  Therefore, 
satisfactory  assessment  of  the  beginning  of  transition  for  approximately 
two-dimensional  boundary  layers  requires  at  least  three  elements 
(a)  adequate  knowledge  of  the  input  disturbance  and  the  corresponding 
boundary-layer  receptivity;  (b)  knowledge  of  the  development  of  the  mean 
profiles  and  access  to  their  stability  characteristics;  (c)  information 
on  the  length  of  the  nonlinear  processes  and  secondary  instability  as 
dependent  on  pressure  gradient,  heat  transfer,  etc.  Since  the  information 
required  in  element  (a)  often  is  not  available  in  the  literature,  one 
usually  characterizes  a disturbance  in  terms  of  of  the  ratio  of  its 

amplitudes  (A  /a  ) at  two  locations,  x,  and  x„. 

2 1 (61  1 2 
According  to  Klebanoff  et  al.,  stage  (1)  is  reached,  for  a flat 

plate,  when  the  rms  velocity  fluctuation  u'  in  the  boundary  layer 

reaches  (u'/U)  *.  .015,  but  that  the  first  appearance  of  turbulence 

max 

spots  is  expected  at  (u'/U)  = 0.2.  That  is,  beyond  the  onset  of  non- 

raaX  25  (51 

linearity  an  amplification  factor  of  10  to  15  times  (=e  * ) is  required.  ' 

Liepmann^2^  hypothesized  that  at  the  breakdown  to  turbulence  the  Reynolds 

stress  t = -puv,  due  to  the  amplified  fluctuations  u',  becomes  comparable 

in  magnitude  to  the  maximum  mean  laminar  shear  stress,  x.  = p 5u/3y  in 

the  boundary  layer. 


t/tl 


{kb(u'/U  )2[a(x)]2} 

e max 


(1) 


where 


b - v'/u' 


and 


(2) 


k ■ u'v'/uv 


(3) 
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where  u and  v are  the  velocities  in  the  x and  y direction,  respectively. 
For  a given  frequency  w*  the  amplification  a is  given  by 


a(x,w*)  = exp[ 


x 


(o  /R) (U  /U  )dx] 
1 e <*> 


(4) 


n 


where  R «-  freestream  length  Reynolds  number,  = (IML/v^) , 

= spatial  amplification  rate, 
x = x*/L, 

2 

<o*  = nondimensional  frequency  = ui  U /v,  and 

e 

L - characteristic  length. 

/ ON 

Smith'  reduced  Liepmann's  criterion  for  transition  to  an  explicit 
dependence  on  the  local  laminar  skin  friction  coefficient,  th°  disturb- 
ance input  at  the  neutral  point  xr,  and  the  total  amplification  ratio 

(A„  /Av  ) where  n refers  to  the  neutral  point  and  t to  the  transition 
*t  xn 

point.  Smith  studied  available  information  on  transition  of  attached 

boundary  layers  where  the  freestream  turbulence  level  was  low.  Assuming 

linear  theory  valid  up  to  the  transition  point.  Smith  showed  that  the 

ratio  of  the  disturbance  amplitude  at  transition  Ax  to  that  at  the 

6 £ g 

neutral  point  A^  is  given  by  (AX{./AXn)  = a(xt,wt>  - e . Later,  more 

accurate  calculations^  showed  /Ax  ) = e^®.  In  any  event,  since 

t n 2 5 

in  the  nonlinear  zone  the  amplification  to  transition  is  =e  , we  find 

that  for  low  freestream  disturbance  levels  the  linear  TS  amplification 

of  about  e ***  does  control  to  a large  extent  the  major  part  of  the 

development  of  the  disturbance  to  the  beginning  of  transition  and  that 

element  (b) , of  the  transition  process,  paragraph  1,  page  11,  appears, 

at  least  in  this  case,  to  dominate  elements  (a)  and  (c). 

In  spite  of  the  dominant  role  of  element  (b),  however,  the  role  of 

element  (a)  remains  extremely  important.  For  example,  when  Spangler 

and  Wells  minimized  sound  disturbances  in  their  measurements  of 

transition  in  a low-speed  boundary-layer  channel,  their  Ry  exceeded 

*t 

five  millions!  These  results,  where  a mixture  of  vorticity  and  sound 
disturbances  arc  present,  cannot  be  predicted  using  any  of  the  present 
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aval  lab]  c forecasting  techniques.  Hence  a knowledge  of  the  receptivity 

of  the  boundary  layer  to  vorticity  and  sound  disturbances  is  needed  for 
further  progress,  and  much  attention  need  be  given  element  (a).  In  the 
absence  of  such  information  and  in  view  of  the  fact  that  the  TS  mechanism 
may  describe,  in  the  absence  of  the  effects  of  surface  roughness,  vibration, 
and  sound,  the  substantial  growlh  of  disturbances  up  to  the  emergence  of 
the  final  three-dimensional  turbulent  spots  and  wedges  and  the  beginning 
of  transition,  it  is  not  unreasonable  to  employ,  for  the  present,  linear 
theory  as  outlined  by  Smith  in  bracketing  the  Reynolds  number  at  the 
beginning  of  transition  R^  for  two-dimensional  and  axisymmetric 
boundary-layer  flows  (in  axisymmetric  flows  x is  replaced  by  s,  the  dis- 
tance measured  along  the  body  surface) . 

Therefore,  in  certain  boundary-layer  flows  where  the  linear  mechanism 
dominates  the  growth  of  disturbances  to  transition,  the  transition  Reynolds 

number  can  be  bound  by  R^  on  the  lower  side  and  by  (R*)  Q on  the  upper 

crit  ev 

side.  Hence,  a plot  of  R and  (R  ) for  two-dimensional  wedge  flows 

xcrit  x e“ 

with  and/or  without  heating  (as  shown  in  Fig.  7)  can  perhaps  be  used  as  a 


guideline  ip  bracketing  R on  a body  of  revolution,  for  example. 

xtrans 

This  can  be  accomplished  by  computing  for  the  body  of  revolution  H vs.  Re^. 

If  this  locus  of  R vs.  H falls  between  the  two  curves  labeled  Rv  . and 
x xcrit 

(Rx)fi9  In  Fig.  7,  then  the  flow  over  the  body  may  be  considered  to  be 
completely  laminar.  If  the  locus  of  R^  vs.  H crosses  the  (R^g  curve, 
then  it  is  assumed  to  undergo  transition  at  the  R^  of  the  intersection 
point.  If  the  locus  lies  very  close  to  the  Rjj  ^ curve,  then  the  body 

is  overheated,  whereas  if  the  locus  lies  very  close  to  the  (R^)  curve, 
then  more  heating  would  be  preferable  (to  maintain  laminar  flow).  These 
remarks,  of  course,  may  not  hold  completely  since,  as  stated  earlier,  in 
heated  flows  H alone  does  not  totally  control  the  stability,  and  hence 
the  transition  behavior  of  the  body.  Further  confidence  in  this  analysis 


can  be  gained  as  some  measurements  of  R,  vs.  H become  available 

xtrans 

and  are  used  to  check  the  validity  of  the  trends  indicated  in  Fig.  7. 


/ 
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Section  VI 


LAMINARI ZAT ION  OF  WATER  BOUNDARY  LAYERS 
OVER  BODIES  OF  REVOLUTION 
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VI-2  LAMINARIZATION  OF  WATER  BOUNDARY  LAYERS  OVER  BODIES 

of  REVOLUTION 

A.M.O.  Smith  and  A.R.  Wazzan 
ENERGY  AND  KINETICS  DEPARTMENT 
University  of  California 
Los  Angeles 

I . INTRODUCTION 

This  paper  is  in  the  nature  of  a progress  report  and  description  of 
the  objectivesof  a contract  with  the  Navy  (N000024-76-C-7071) . The  three 
main  task  statements  of  the  contract  are: 

A.  Establish  and  evaluate  design  configurations 

B.  Determine  the  effect  of  pressure  gradient,  body  shape,  surface  heating 
rates,  wall  temperature  distribution  and  suction  on  boundary  layer 
stability  and  transition  characteristics. 

C.  Calculate  and  analyze  the  effects  of  free  stream  turbulence,  wall  rough- 
ness, vibration  and  acoustics  disturbances. 

Statement  B mentions  the  effect  of  pressure  gradient.  The  effect  is 

being  examined  in  two  ways  (1)  by  the  systematic  family  of  bodies  to  be 

described  and  (2)  by  study  of  the  two  dimensional  wedge  flow  family.  Some 

of  the  results  of  this  study  will  be  reported  separately  in  Wazzan 's  and 

Gazley's  paper.  Concerning  Item  C,  while  we  are  interested  in  the  subject 

and  it  is  important,  the  main  responsibility  for  studying  the  problem  has 

been  given  to  Dr.  Leslie  Mack  under  subcontract. 

Therefore  our  attention  here  will  be  with  a family  of  bodies  that  we 

have  developed  and  are  now  studying.  All  the  transition  studies  have  used 

g 

the  e method.  Hence  the  value  of  this  study  stands  or  falls  with  the 

9 

accuracy  of  the  e method,  which  is  reviewed  in  the  Wazzan-Gazley  paper. 
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Since  no  shortcut  transition  method  has  been  found  for  bodies  of  revolution, 
the  systematic  family  of  shapes  at  least  offers  the  chance  for  interpolation. 
The  study  is  engineering  oriented,  trying  to  show  general  answers  based  on 
the  state  of  the  art  of  the  theory  rather  than  trying  to  advance  a special 
piece  of  theory  - Les  Mack's  part  being  an  exception.  The  project  started 
about  Feb.  1 and  lasts  for  a year,  so  it  is  a bit  more  than  halfway  through. 
Also  it  might  be  said  that  the  work  is  not  necessarily  new.  Many  similar 
studies  have  been  made  but  results  are  not  available  in  the  public  domain. 
Furthermore,  it  is  hoped  that  the  studies  will  be  more  systematic  and 
comprehensive  than  those  of  the  past.  Therefore  comprehensiveness  and 
general  availability  are  the  justification  for  this  investigation. 

II.  THE  FAMILY  OF  BODIES 

We  will  first  describe  the  general  characteristics  of  the  whole  family 
and  then  in  III  describe  certain  specific  attributes  found  so  far.  The 
family  of  six  bodies  is  shown  in  Figures  1 through  6.  The  family  is 
capable  of  giving  some  idea  of  the  effect  of  fineness  ratio,  pressure 
gradient  and  prismatic  coefficient  upon  performance.  Except  for  PR-9  all 
are  derivaties  of  the  Reichardt  body,  hence  the  names.  Since  there  was  much 
interest  in  a 9:1  fineness  ratio  body,  REIK-9,  Fig.  1,  was  chosen  as  one  of 
the  bodies,  the  9 being  for  fineness  ratio.  Figures  1-6  are  complete  in  a 
sense.  They  show  the  shape,  the  velocity  distribution,  the  prismatic 
coefficient,  volume,  and  the  basic  defining  equations.  The  basic  defining 
equations  are  not  written  in  terms  of  x/l  but  instead  use  the  coordinate 
system  shown  on  the  body.  A tail  fairing  is  added  to  the  basic  body.  The 
junction  point  and  the  length  are  noted  in  the  drawing.  After  working  out 
the  shapes  in  these  simple  coordinates  the  bodies  were  then  scaled  into  the 
x/l  system.  Since  Reichardt 's  body  was  defined  by  one  simple  equation. 
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except  for  the  tail,  this  course  was  followed  for  all  the  others  except 
PR-9.  For  the  bodies  with  favorable  pressure  gradient  like  REIK-4  a 
considerable  number  of  modification  methods  were  examined  but  since  we  do 
not  know  exactly  what  we  want  as  a pressure  distribution,  the  simplest  factor, 
(1  + ax)  that  gave  reasonable  shape  variations  was  chosen.  REIK-9  and 
REIk-2  show  the  effect  of  thickness  ratio.  REIK-2,  3,  4,  and  5 show  the 
effect  of  pressure  gradient,  REIK-4  and  5 show  the  effect  of  pressure  gradient 
and  thickness. 

PR-9  does  not  fall  into  the  REIK-family.  There  is  so  much  advantage  to 
high  prismatic  coefficients  that  it  was  decided  to  see  what  can  be  done  by 
heat  in  overcoming  adverse  pressure  gradients.  A cylindrical  body  always 
has  adverse  gradients  at  the  nose,  which  normally  are  ruinous  to  laminar 
flow.  The  goal  therefore  was  to  take  a desirable  practical  shape  and  see 
whether  it  is  possible  by  judicious  use  of  heat  to  overcome  the  injurious 
effects  of  the  adverse  nose  gradients.  The  nose  is  a modified  leminiscate 
The  junction  point  is  noted.  The  leminiscate  matched  values  and  first  and 
second  derivatives  at  the  junction  point.  This  body  is  of  9:1  fineness  ratio. 

Velocity  distribution  for  all  the  bodies  were  calculated  by  Hess's 
Higher  Order  Neumann  Method. 

III.  SPECIFIC  BODIES 

Because  of  Navy  interest  and  because  of  uncertainties  about  scaling,  a 
basic  condition  for  calculation  was  U^/v  * 6.86  x 10^  per  ft.  This  number 
corresponds  to  75  fps  in  67°F  fresh  water.  67°F  fresh  water  was  used  because 
that  is  representative  of  a towing  tank.  Also  a specific  length  was  used, 
9.5625  ft.,  giving  a length  Reynolds  number  of  65.6  x 10**. 
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A.  REIK-9 

Figures  7,  8,  9 and  10  are  representative  results  for  this  body.  The 

object  was  to  stabilize  the  boundary  layer  to  get  laminar  flow  clear  back 

to  separation.  Because  the  velocity  distribution  is  almost  perfectly  flat, 

cold  transition  is  about  that  of  a flat  plate,  i.e.,  RgTr  = 5 x 10^.  This 

number  is  for  £n  A/Aq  = 9.  Now  in  order  to  have  a margin  of  safety  and  to 

stay  in  the  linear  range  of  Tollmien-Schlichting  disturbances,  when  we  are 

looking  for  heat  requirements  we  have  chosen  to  try  to  hold  the  amplification 
7 

ratio  below  e . Figures  7 and  8 are  two  results.  Figure  7 is  the  first  try 
and  Fig.  8 the  fourth  and  successful  try.  What  kind  of  heat  distribution  to 
use  is  an  unknown  science.  Such  considerations  as  the  growth  of  T-S  waves 
cold,  where  they  begin  growing,  whether  they  stop  growing,  shape  factor  and 
the  like  all  enter  in.  The  first  heat  distribution  saw  the  growth  of  T-S 
waves  begin  at  about  S = 0.35  ft.  and  the  growth  exceeded  the  desired  7.0. 
Then  after  two  intermediate  tries  VAR  IV,  Fig.  8 was  tried  and  the  goal  of 
e^  was  successfully  met. 

It  has  generally  been  more  convenient  and  realistic  to  specify  heat 
input  rather  than  wall  temperature.  Figure  9 shows  results  for  the  two 
cases  just  studied.  Note  that  the  beginning  of  heating  in  VAR  IV  was 
slightly  further  ahead  than  in  VAR  I.  Examination  of  the  amplification 
curves  indicated  more  heat  was  needed  forward  and  less  after.  Figure  10 
shows  the  wall  temperature  distribution  for  the  two  cases  as  well  as  the 
total  heat. 

Out  of  curiosity  for  what  it  might  show  we  devised  a dimensionless 
heating  coefficient,  C^,  (H  for  heat)  defined  as  follows: 

. S = (energy  in  to  wat_erj.  (ft  . lbs/sec) 

H l/2pV3  (Vol)Z ' 3 
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Values  are  shown  on  Fig.  9. 

A different  kind  of  result  is  shown  in  Fig.  10.  Assuming  = U^k/V  = 25 
we  calculated  allowable  heights  of  discrete  roughness  for  REIK-9  moving  at 
Hx/V  = 6*86  x 106/ft  = 75  fps.  The  values  are  very  small.  Heat,  as  you 
see,  hurts  slightly.  However  at  1/4  speed,  i.e.,  about  19  fps,  roughness 
can  be  considerable.  The  most  sensitive  part  of  these  bodies  is  just  behind 
the  stagnation  point,  where  there  is  a minimum.  Now  it  is  found  that  this 
minimum  always  seems  to  be  ahead  of  the  point  where  heat  is  first  applied. 
Therefore  the  heat  does  not  affect  the  minimum  and  the  minimum  can  be 
determined  by  cold  calculations  alone,  see  Fig.  10. 

B.  RE IK- 2 

While  calculations  have  been  done,  it  is  a rather  uninteresting  body 
again  having  8=0  over  most  of  its  length.  Therefore  it  will  not  be 
discussed  further,  except  to  say  that  if  the  flow  is  laminar  to  x/l  = 0.7 
at  Rg=  20  x 106  separation  calculates  to  be  at  x/fe  = 0.96. 

C.  REIK-3,  4 

Some  calculations  have  been  done,  but  work  on  these  bodies  is  mostly 
in  the  future.  REIK-2,3,4,5  form  a family  having  increasingly  favorable 
pressure  gradients.  For  the  four  respectively,  they  develop  maximum 
Hartree  8's  near  x/&  = .7  of  0,  0.)  '2  and  0.3  approximately.  However, 

none  of  these  8's  are  sufficiently  high  to  ge‘  i./  real  extent  of  laminar 
flow  at  high  Reynolds  numbers,  without  the  assistance  of  heat. 

D.  REIK-5 

Because  it  is  doubtful  that  the  full  study  will  be  completed  because  of 
time  and  money  limitations  we  have  chosen  to  look  at  extremes  so  as  to  »ee 
what  the  boundaries  are.  Then  later  as  time  and  money  allow  we  can  fill  in 
with  intermediate  cases.  One  might  say  we  are  aiming  for  two-point  interpola- 
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tion  and  then  later  as  conditions  allow,  proceed  to  three-point.  Therefore 
even  though  not  identified  with  any  practical  concept  we  looked  next  at 
REIK-5 . 

Only  two  heat  distribution  tries  were  made.  While  the  fa  orable  pres- 
sure gradient  alone  will  not  give  an  extended  run  of  laminar  flow,  it  does 
indeed  help.  As  a result  the  first  try  led  to  excessive  stability  and  a 

second  try  was  made.  The  result  is  shown  in  Fig.  11.  It  is  still  too  stable 

7 

but  the  amplification  is  near  enough  e that  this  try  has  been  taken  as  final, 
although  we  know  that  better  can  be  done.  The  heating  distribution,  temper- 
ature, and  shape  factor  are  shown  in  Fig.  12.  The  heat  intensity  is  much 
less  than  that  found  needed  for  REIK-9.  Also  the  total  heat  requirement 
is  low  although  the  body  has  a much  larger  volume.  Notice  that  the  shape 
factor  becomes  very  low,  2.32.  Also  that  Cu  = .0C78,  or  less  than  1/4  of 
that  on  REIK-9.  This  body  was  deliberately  made  thick,  knowing  that 
separation  would  occur  at  the  tail.  The  tail  would  either  have  to  be 
lengthened  or  suction  would  have  to  be  used  to  prevent  separation.  At 
Rp  = 20  x 10^  and  laminar  flow  to  x/&  = 0.70  separation  occurs  at  x/£  = .825. 
Some  of  the  boundary  layer  parameters  for  the  cold  case  are  shown  in  Fig.  13. 

A roughness  calculation  was  done  for  the  cold  case  to  compare  with  the 
REIK-9.  No  surprises  were  expected  but  one  was  found.  On  REIK-9  the 

minimum  is  very  far  forward  and  after  that  the  allowable  value  slowly  but 

1/4 

definitely  grows.  (On  a flat  plate  it  would  grow  as  x ).  On  REIK-5 
there  is  indeed  a minimum  very  near  the  nose  but  a large  portion  of  the 
body  has  a higher  but  nearly  constant  level.  The  reason  is  that  an 
expanding  body  thins  the  boundary  layer.  Hence  bodies  like  this  should 
be  slightly  more  sensitive  to  roughness  than  a straight  Reichardt  type. 
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E.  PR-9 

This  shape  has  adverse  gradients  near  the  nose.  Now  the  name  of  the 

game  is  to  see  what  can  be  done  with  heat.  Notice  on  Fig.  15  that  3 

reaches  a value  of  -.135  near  the  nose.  This  figure  shows  the  third  try. 

2 

Notice  that  = 20  Kw/ft  near  the  nose  and  that  wall  temperatures  reach 

nearly  600°R.  Figure  16  shows  amplification  growth,  which  is  a little  too 

much.  To  try  to  reduce  it  slightly  the  last  part  of  the  heat  distribution 

2 

was  held  constant  at  2.5  Kw/ft  , instead  of  reducing  linearly  to  zero  at 

6.5  ft.  Compare  Figs  15  and  17.  This  small  change  led  to  large  changes 

in  the  wave  growth.  The  maximum  was  reduced  more  than  desired  from 

Jin  A/Ao  equal  to  7.9  in  Fig.  16  to  3.9  in  Fig.  18.  Hence  according  to  our 

calculation  heat  is  making  a poor  shape  work.  But  a larger  quantity  is 

necessary  than  for  a good  shape.  For  REIK-9,  = 0.0334,  for  PR-9, 

C„  = 0.0671.  That  is,  for  the  same  fineness  ratio  body,  over  twice  as 
H 

much  heat  was  necessary  to  fix  up  the  poorly  shaped  body.  For  REIK-5 
which  has  the  most  favorable  shape  CH  = 0.0078.  It  is  easily  shown  that 
the  following  formula  applies. 

C 

Q AVAILABLE  1 " "th  D Vol 

Q NEEDED  = np  nth  * CH 

where  Q is  the  heat  rate,  n^  the  thermal  efficiency  of  the  engine  that 
runs  the  device  and  np  the  propeller  efficiency.  The  ratio  of  course 
must  be  greater  than  unity. 

IV.  SURFACE  SMOOTHNESS 

The  last  item  we  present  is  our  work  on  the  subject  of  allowable  surface 
roughness.  First  we  will  discuss  Westinghouse's  recommended  line,  see 
Fig.  19.  Westinghouse  has  assumed  that  surface  smoothness  acts  the  same 
way  as  discrete  roughness,  that  is  a small  roughness  develops  an  unsteady 
wake  behind  it.  They  basically  use  the  number  R^  « 25.  Then  to  relate 
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surface  smoothness  to  this  number  they  estimate  that  there  will  be  a reasonable 

number  of  roughnesses  which  are  5 times  the  statistical  RMS  value.  Then 

they  assume  that  the  velocity  in  the  boundary  layer  at  the  top  of  the 

roughness  is  1/5  of  the  velocity  at  the  edge.  Putting  all  together  they 

getCU^/v  x RMS  < 25  but  actually  they  use  20.  This  is  their  line  on  Fig.  19. 

Now  it  can  be  shown  that  if  we  follow  their  basic  R^  assumption  more 

-3/4 

carefully  we  get  the(Uoo/v)  line.  Let  us  assume  that  R^  is  to  be  a 
constant,  say  25,  and  for  simplicity  let  us  look  at  a flat  plate  flow. 

The  flow  will  be  as  sketched  below.  A boundary 


layer  approaches  the  roughness  as  shown.  The  roughness  is  very  low  in 
the  boundary  layer  so  the  profile  will  be  assumed  to  be  straight  in  this 
region,  an  excellent  assumption.  Then  we  can  write,  for  a flat  plate 

V k ( 3~y  ^ w = kUJ  2vx  _ F"(0)  = . 332kv[ ^ 


\ k 

Now  using  Rj^  = — — 
to  obtain 

»k 


25  we  can  introduce  the  above  expression  for 


Uk 


.332  k2U 
v 


25 


Now  if  we  solve  for  k we  find 


3/4 


- - -1/4  (*)' 
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The  difference  from  Westinghouse  analysis  is  that  we  solved  for 

instead  of  assuming  it.  The(U/v)-3^line  in  Fig.  19  was  simply  put 

through  the  same  point  as  Westinghouse  at  U/v  = 106,  Hence  the  authors 

feel  that  at  higher  Reynolds  numbers  per  foot,  the  Westinghouse  criteria 

i is  unduly  conservative. 

However,  we  are  not  sure  the  above  is  the  right  approach.  The  surface 

may  be  so  smooth  that  no  unsteady  wakes  are  formed.  Then  a slightly  rough 

surface  "tickles"  the  boundary  layer  somewhat  like  small  external  turbulence 

does.  If  that  is  the  case,  in  view  of  stability  theory  the  characteristic 

Reynolds  number  would  be  a thickness  Reynolds  number  like  R^  instead  of 

a length  Reynolds  number.  Now  . According  to  this  theory  then, 

-1/2 

a line  should  slope  asQJ^/v)  . This  is  the  third  line  plotted  in 

Fig.  19.  There  is  hardly  any  data  to  verify  these  lines.  The  NACA  points 

are  those  cited  by  Westinghouse.  At  Douglas  about  1956  we  roughened  a surface 

with  several  grades  of  sand  paper,  measured  the  roughness  with  a profilometer 

and  found  the  effect  on  transition.  The  tests  were  far  better  designed  than 

those  of  the  NACA.  The  points  are  shown.  The  point  at  U^/v  <=  1.1  x 106 

* —1/2 

did  not  affect  transition.  The  other  two  seem  to  confirm  the  (U  /v) 

00  w 

slope  but  there  is  too  little  data  to  be  sure.  In  any  case  thefU^/v)"1 
curve  is  not  right. 
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Prismatic  coeff.  = 0.4965  Volume  = 0.015207 


cos2  (0+45°)  X^pcosd,  Y'=?sin  0 -3O°^0£O° 
o = 4.4669000 ( 0. 80592 74-X')  , r/r0=(  1-2.1491400Y' 
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VI -3  LOW- SPEED  BOUNDARY  LAYER  TRANSITION  WORKSHOP 
(Sept.  13-14,  1976:  Rand  Corp.) 


A STUDY  OF  HEAT  ADDITION  EFFECTS  ON  LAMINAR, 
WATER  BOUNDARY  LAYERS  OF  AXISYMMETRIC  BODIES 


by 

Joseph  J.  Eisenhuth 
Applied  Research  Laboratory 
The  Pennsylvania  State  University 
P.  0.  Box  30 
State  College,  PA  16801 


Summary  of  Presentation 

The  presently  popular  subject  of  stabilizing  boundary  layers  by  the  addition 
of  heat  at  the  wall  has  received  attention  at  the  Applied  Research  Laboratory  (ARL) , 
Penn  State.  The  primary  interest  here  has  been  in  the  prediction  of  the  heat 
necessary,  from  a designer's  point  of  view,  to  forestall  transition  on  axisymmetric 
bodies  and  in  the  estimation  of  the  accompanying  mean-flow  boundary  layer  character- 
istics that  result.  The  focal  point  in  this  investigation  has  been  the  implementation 
of  a code  to  provide  a solution,  by  numerical  integration  means,  of  the  boundary 
layer  equations  and  in  getting  the  kind  of  correlation  of  critical  Reynolds  number 
that  would  permit  bypassing  the  very  involved  stability  calculations  in  making 
simple  estimates  of  transition. 

Basic  Equations  and  Computer  Code 

The  basic  boundary  layer  equations  for  axisymmetric  flow  are  shown  in  Slide  1 
and  are  derived  for  an  incompressible  flow  with  heat  addition  at  the  boundary.  In 
the  development  of  the  equations,  effects  such  as  buoyancy  and  dissipation  of  energy 
by  friction  are  ignored.  The  boundary  conditions  are  shown  at  the  bottom  of  the 


slide. 
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To  put  the  equations  In  a form  suitable  for  solution,  a stream  function,  ip, 
and  a nondimens ional  normal  distance,  H,  are  defined  in  Slide  2.  In  addition, 
the  M,  R,  and  9 parameters,  which  will  appear  in  the  resulting  equations,  are  also 
defined.  M is  the  standard  Falkner-Skan  parameter,  R is  a radius  gradient  para- 
meter which  results  from  the  axisymmetric  case,  and  0 is  just  the  normalized 
temperature  difference. 

The  resulting  equations  with  boundary  conditions  are  shown  in  Slide  3.  These 
ordinary  differential  equations  are  explicitly  independent  of  x and  the  velocity 
and  temperature  profiles  that  result  are  given  by  similar  type  solutions.  The 
primes  indicate  derivatives  with  respect  to  n;  the  bars  over  the  physical 
quantities  show  normalization  with  respect  to  the  free-stream  quantities,  P is 
the  free-stream  Prandtl  number;  Cp  is  the  specific  heat  at  constant  pressure;  k 
is  the  thermal  conductivity. 

The  code  for  the  solution  of  the  two-dimensional  version  of  the  mean-flow 
equations  appearing  in  Slide  3 is  presented  in  Reference  1 as  part  of  a more 
general  computer  code.  In  this  solution  the  method  of  Nachtsheim  and  Swigert 
(Reference  2)  was  adapted  for  the  numerical  integration  process  and  was  found 
to  be  satisfactory.  In  the  present  study,  the  code  for  the  two-dimensional 
equations  was  extended  to  the  axisymmetric  case  and  was  further  altered  to 
provide  certain  boundary  layer  information  that  wasn't  explicitly  calculated 
in  the  original  version. 

Critical  Reynolds  Number  Correlation 

Stabilitv  information  in  terms  of  the  critical  Reynolds  number,  Rt*  , 

° crit 

is  available  for  two-dimensional  flow  both  for  the  case  of  the  unheated  boundary 
with  pressure  gradient  and  for  the  flow  over  a flat  plate  with  heating  at  the 
wall.  The  results  of  the  unheated  case  are  found  in  Reference  3 as  the  critical 
Reynolds  number  versus  the  Pohlhausen  parameter  A - , 


and  the  heated  case 
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results  were  obtained  from  Reference  1 as  critical  Reynolds  number  versus  temperature 

difference.  A and  AT  were  translated  into  the  shape  parameter  H = 6*/9  so  that  in 

each  case  Rt*  could  be  plotted  against  H.  The  curves  in  Slide  4 are  the  result. 

° crit 

It  can  be  seen  that  the  correlation  is  quite  good,  indicating  that  the  stability 
of  the  laminar  boundary  layer  is  strongly  dependent  upon  H,  regardless  of  whether 
it  is  obtained  by  favorable  pressure  gradient  or  by  heat.  A similar  correlation 
is  reported  in  Reference  4. 

This  correlation  provides  the  basis  for  the  development  of  the  type  of  design 
information  which  will  be  presented  here.  Although  based  on  two-dimensional 
stability  information,  the  correlation  is  considered  valid  for  the  axisymmetric 
case  because  the  stability  equations  are  the  same  for  both  cases. 

Computed  Results 

There  are  a variety  of  ways  to  present  results  depending  upon  the  user's 

needs.  Slide  5 presents  curves  showing  the  effect  of  pressure  gradient  and 

temperature  on  the  shape  parameter  for  the  two-dimensional  case  (R=0) . If  one 

wishes  to  obtain  a certain  shape  factor,  one  need  only  interpolate  between  the 

curves  of  constant  M to  determine  the  temperature  difference  necessary  to  do  so. 

If  one  would  rather  work  in  terms  of  an  effective  M,  the  additional  curve  of  H 

versus  M in  Slide  5 could  be  used  along  with  the  idea  of  the  shape  factor-stability 

correlation  to  generate  the  curves  shown  in  Slide  6.  With  the  curves  of  Slide  6, 

one  need  only  interpolate  between  curves  to  arrive  at  a temperature  difference 

to  produce  a chosen  effective  pressure  gradient  parameter.  In  lieu  of  trying  to 

work  with  temperature  as  an  additional  parameter,  it  was  felt  that  the  effective  M 

could  be  used  in  some  simplified  transition  prediction  method.  As  an  example  of 

local  shear  stress  and  local  heat  flux  values  that  result  from  a choice  of 

stability  level.  Slides  7 and  8 are  presented  for  the  specific  values  of  Ug  = 50 

ft/sec,  x ■ 4.0  ft,  and  again  the  two-dimensional  case.  The  values  of  shear  stress 

and  heat  flux  can  be  scaled  from  these  curves  for  other  values  of  x and  U . 

e 

tk  ’ 

. : - - — * 
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Slides  9 and  10  are  presented  to  show  the  effect  of  having  different 
ambient  temperatures.  The  interest  in  this  kind  of  effect  stemmed  from  consider- 
ations of  a water  tunnel  test  program  in  which  the  Reynolds  number  range  could  be 
extended  by  heating  the  tunnel  water.  Slide  11  is  developed  to  provide  values 
of  temperature  difference  to  maintain  a certain  level  of  effective  M.  In  this 
case  the  =0.3  level  represents  close  to  the  maximum  heat  that  can  be  added 

before  the  critical  Reynolds  number  will  start  decreasing. 

Slides  12  and  13  are  presented  to  show  the  difference  between  the  two- 
dimensional  and  axisymmetric  cases.  Similar  plots  can  be  shown  for  other 
values  of  R.  Slide  14  is  similar  to  Slide  11  in  that  it  shows  the  temperature 
difference  necessary  to  maintain  a certain  level  of  stability,  in  this  case  for 
various  values  of  R and  M.  The  level  chosen  here  was  again  close  to  the  maximum 
heat  that  can  be  added  before  the  critical  Reynolds  number  will  again  decrease. 
Laminar  Separation  Limits 

An  attempt  was  made  to  obtain  some  limits,  with  regard  to  the  various 
parameters,  beyond  which  laminar  separation  would  occur.  Slide  15  represents  one 
way  of  pr  ■ ing  this  information.  Values  of  R and  M below  and  to  the  left  of  the 
curves  represent  separated  flows.  Another  presentation  of  this  same  information 
is  found  in  Slide  16.  It  can  be  seen  that  temperature  difference  has  relatively 
little  effect  on  laminar  separation.  This  is  consistent  with  the  calculations 
reported  in  Reference  5 for  the  two-dimensional  case  and  are  shown  as  the  curve 
marked  "Rand".  Although  showing  slightly  more  effect  of  temperature,  the  Rand 
calculations  also  indicate  the  minor  role  of  heating  in  the  prevention  of  laminar 
separation. 

Transition  Criteria 

A purpose  in  employing  the  R,*  - H correlation  was  in  being  able  to  use 

0 crit 

the  effective  M to  enter  experimental  transition  curves  based  on  unheated  operation. 
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Presumably,  theoretical  stability  information  for  the  unheated  case  might  also  be 
reduced  to  a form  that  would  make  use  of  the  effective  Falkner-Skan  parameter. 
Aspects  of  this  work  are  still  in  progress,  so  only  a first  attempt  at  developing 
some  simplified  transition  curves  is  presented. 

It  would  seem  logical  that  the  parameters  that  affect  H and  thus  R^*  , 

0 crit 

would  also  be  controlling  as  far  as  transition  is  concerned.  The  first  attempt 
therefore  made  use  of  the  parameters  M and  R.  Reference  6 provides  transition  • 
information  for  several  bodies  over  a range  of  Reynolds  numbers.  The  geometric 
characteristics  and  the  potential  flow  calculations  for  one  of  these  bodies  were 
used  to  generate  R and  M distributions  along  the  body  length.  Combining  this 
information  with  the  experimental  transition  curve  given  for  that  body,  a plot 
like  the  one  shown  in  Slide  17  is  possible.  Very  little  variation  of  M along  the 
body  occurs  and  the  major  parameter  appears  to  be  R.  This  is  similar  to  the  kind 
of  criterion  used  in  Reference  7 except  that  instead  of  R a parameter,  (D/r) (dr/dJl) 
is  used.  It  is  anticipated  that  M can  be  an  equally  effective  parameter  in 
influencing  transition  and  that  bodies  with  greater  variations  in  M would  show 
different  trends.  In  other  words,  there  should  be  some  cross  curves  of  constant 
M appearing  in  a general  plot  for  all  types  of  bodies. 
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Vi-4  DESIGN  OF  A LAMINAR,  LOW-DRAG  RESEARCH  VEHICLE 

R.  S.  Scotti , B.  H.  Carmichael,  and  D.  A.  King, 
Autonetics  Group,  Rockwell  International,  Anaheim,  CA 


An  effort,  sponsored  by  DARPA,  is  underway  to  perform  the  preliminary 
design  of  a low-drag  research  vehicle  with  a laminar  forebody  at  a transition 
length  Reynolds  number  to  minimum  pressure  of  250  x 102 * * * 6.  Boundary  layer 
stability  is  to  be  maintained  by  exploiting  the  effects  of  favorable  pressure 
gradient  (shape)  and  distributed  surface  heating.  The  present  level  of  under- 
standing of  these  two  effects,  particularly  in  combination,  does  not  allow 
for  their  simultaneous  optimization  without  subjective  inputs. 

The  accepted  procedure  for  determining  body  shape  to  maintain  (or 
stabilize)  laminar  boundary  layer  flow  is  one  of  iterative  numerical  compu- 
tational trials.  The  flow  field  and  corresponding  stability  characteristics 
must  be  computed  for  each  trial  shape.  The  addition  of  surface  heating  adds 
new  dimensions  to  the  task  of  optimizing  vehicle  contours  by  introducing 
several  new  variables  such  that  the  best  unheated  contour  is  probably  not 
the  best  heated  contour. 

A numerical  study  is  being  performed  at  Autonetics  following  the  above 
iterative  approach.  Guidance  and  overall  direction  is  being  provided  by 
extensive  experience  with  laminar  flows  possessed  by  the  second  author,  by 
up-to-date  results  of  laminar  flow  research,  and  by  intuition.  The  following 
specific  criteria  in  order  of  decreasing  importance  are  applied  in  the 
process  of  evaluating  candidate  designs: 

In  regard  to  shape, 

1.  Primary  objective  is  to  achieve  a transition  Reynolds  number  (based 
on  length  and  local  free-stream  velocity)  of  250  x 10&.  Design  is 
to  be  conservative  as  judged  by  all  available  transition  predictive 
methods.  “Conservative"  is  quantified  by  the  next  item. 

2.  Allow  for  flexibility  of  both  vehicle  operating  conditions  and  test 
plan.  The  design  must  accommodate  surface  temperature  distributions 
such  that, 

0 < AT(x)  5 ATmax(x) 

where  ATmax(x)  is  in  excess  of  that  predicted  to  completely  damp 

growth  of  all  small  disturbances.  Moreover,  the  design  must  allow 
for  a progressive  approach  to  the  Re-p  objective  by  means  of  a 

systematic  variation  of  both  vehicle  speed  and  surface  temperature 
field  from  lower  operational  transition  Reynolds  number  to  250  x 10&. 
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3.  Maintain  a prismatic  coefficient  which  provides  sufficient  volume, 
particularly  in  the  vehicle  nose  region,  to  practically  accommodate 
the  surface  heat  exchanger.  Also,  low  sensitivity  to  angle  of 
attack  is  to  be  sought  in  a rounded  nose  design. 

In  regard  to  size, 

4.  Maximize  test  time  for  a bouyancy-driven  vehicle  operating  from  a 
fixed  release  depth,  while  minimizing  the  effects  of  residual 
unsteadiness  (acceleration)  during  "test"  and  the  problem  of 
broaching  at  the  surface. 

5.  Minimize  the  effects  of  surface  roughness  and  waviness  within 
acceptable  and  reasonable  (from  the  point  of  view  of  fabrication) 

1 imits. 

6.  Maintain  contact  with  feasible/practical  engineering  practices  in 
terms  of 

-transportability 

-maneuverability 

-producibility 

The  stabilizing  influences  of  a local  favorable  pressure  gradient  are 
controlled  by-  vehicle  contour  and  fineness  ratio.  Those  of  surface  overheat 
are  apparently  controlled  by  both  the  magnitude  and  the  gradient  of  the 
temperature  field.  In  the  present  study  all  four  variables,  i.e.,  contour, 
fineness  ratio,  magnitude  of  and  gradient  of  temperature  are  being  investi- 
gated. 

Results  in  the  form  of  boundary  layer  flow  and  corresponding  stability 
characteristics  are  presented  for  several  promising  shapes.  These  are 
compared  to  those  of  other,  earlier  conceptual  designs  in  Figures  land  2— 
the  first  being  a plot  of  some  vehicle  shapes  and  corresponding  pressure 
coefficients,  while  the  second  displays  the  boundary  layer  shape  factor, 

H = <5*/6,  over  the  surface  of  several  vehicles.  Also  shown  in  Figure  1 (in 
brackets)  are  the  disturbance  amplitude  ratios  as  computed  by  Prof.  F.  Wazzan 
for  the  most  amplified  Tol lmien-Schl ichting  waves  in  the  boundary  layer  of 
each  body.  Indications  are  that  a fair  measure  of  boundary  layer  stability 
is  provided  by  the  magnitude  and  history  along  the  vehicle  surface  of  the 
shape  factor  H. 

The  vehicle  contour  designated  as  AH-12,  which  has  a rounded  nose  like 
the  AH-9  (as  compared  to  the  pointed  AH-8)  and  a fineness  ratio  of  4.0  (as 
compared  to  5.0  for  the  vehicles  shown  on  Figure  1)  apparently  offers 
significant  improvements  in  terms  of  the  criteria  listed  above  over  others 
previously  studied.  Final  refinement  of  this  design  is  in  process. 


• _____  . 
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VI-5  AXISYMMETRIC  VEHICLE  SYSTEM  DESIGN  FOR  TOTAL  PROPULSION  POWER  OPTIMIZATION 

F.  R.  Goldschmied 

Westinghouse  Research  Laboratories 
Pittsburgh,  Pennsylvania  15235 

The  design  objective  is  the  power  minimization  for  a given  useful  hull 
volume  and  operating  speed,  since  it  is  the  total  propulsion-power  of  the  under- 
seas vehicle  that  is  considered,  rather  than  the  drag  of  the  bare  hull.  The 
design  system  optimization  includes  also  the  fins  (to  provide  the  desired  degree 
of  stability  and  maneuverability)  with  their  drag  and  the  propulsor  with  its 
efficiency. 

Two  design  optimizations  have  been  approached: 

a)  Utmost  exploitation  of  laminar  boundary-layers,  whenever  mission 
and  economic  constraints  appear  to  allow  it. 

b)  Assumed  all-turbulent  boundary-layers,  when  no  useful  laminar 
extent  can  be  expected. 

The  design  method  is  based  on  the  optimum  integration  of  hull  design, 
boundary-layer  control  (single  suction  slot)  and  stern  jet  propulsion,  as  pre- 
sented by  Goldschmied.^  A 60"  test  model  is  shown  in  the  wind-tunnel  in  Fig. 

1.  Laminar  transition  is  shown  at  70%  chord  by  the  China  Clay  technique.  The 
suction  slot  is  visible  at  83%  chord;  the  BLC  flow  is  ducted  out  through  the 
hollow  support  strut.  The  experimental  pressure  distribution  at  +6°  angle  of 
attack  is  shown  in  Fig.  2. 

The  laminar /turbulent  axlsymmetric  transition  prediction  method  is  that 

(2) 

given  by  Goldschmied . 

At  the  length  Reynolds  Number  of  10^,  assuming  the  same  tail-fins 
drag  ACq  * 0.003  and  the  same  propulsor  efficiency  of  85%  for  all  cases,  the 
following  evaluation  can  be  made: 

(3) 

a)  Optimized  theoretical  all-turbulent  Parsons/Goodson/Goldschmied 


-263- 

1-36  body  with  tail  fins  and  stern  wake  propeller:  power  coefficient 
CHp  = 0.0295  (Hull  CD  = 0.022). 

b)  Wind-tunnel  test  of  all  turbulent  vehicle,  with  tail-fins  and 
ducted  impeller:  power  coefficient  C = 0.0145. 

tlx 

(3) 

c)  Optimized  theoretical  laminar  Parsons/Goodson/Goldschmied''  ' 

X-35  body  with  tail-fins  and  stern  wake-propeller:  power  coefficient  CUTJ  = 

HP 

0. 0120. (Hull  CD  = 0.0071). 

d)  Wind-tunnel  test  of  laminar  (72%  length)  vehicle,  with  tail-fins 

and  ducted  impeller:  power  coefficient  Cuo  = 0.0106. 

HP 

It  can  be  readily  seen  that  the  system  design  approach  has  yielded 
a 51%  power  reduction  for  the  all-turbulent  case  and  a 12%  power  reduction 
for  the  laminar  case.  Furthermore,  it  is  quite  significant  that  the  all- 
turbulent  vehicle  requires  only  21%  more  power  than  the  laminar  X-35  body. 

The  power  coefficients  are  plotted  against  length  Reynolds  number  in  Fig.  3. 

Considerably  better  results  may  be  obtained  from  a numerical  optimi- 
zation study;  mathematical  models  are  available  for  the  hull  geometry,  the  single 
suction-slot,  the  axial  pump  and  the  stern  propulsive  jet.  For  instance,  five 
parameters  may  be  allocated  to  the  hull  profile,  two  parameters  to  the  suction 
slot,  two  parameters  to  the  pump  and  one  parameter  to  the  jet;  the  ten  parameters 
may  then  be  simultaneously  optimized  by  the  methods  of  Ref.  3,  subject  to  a set 
of  suitable  constraints. 
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Fig.  3 - Power  coefficients  summary  plot 
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SIMPLE  RELATIONS  FOR  THE  STABILITY  OF  HEATED  LAMINAR  BOUNDARY 
LAYERS  IN  WATER:  MODIFIED  DUNN-LINN  METHOD 


I 


J.  Aroesty,  W.  S.  King,  G.  M.  Harpole, 

W.  Matyskiela,  A.  R.  Wazzan,  C.  Gazley,  Jr. 
The  Rand  Corporation 
Santa  Monica,  California 

ABSTRACT 


We  have  developed  a modification  of  Lin's  method  which  can  be  used 
to  estimate  the  minimum  critical  Reynolds  number  of  a heated  laminar 
boundary  layer  in  water.  This  method  is  both  simple  and  surprisingly 
accurate.  As  in  Lin's  original  work,  the  method  is  most  applicable 
when  c •*  0,  a ■*  0 and  aRe  <*>.  The  internal  details  of  the  method  do 
not  agree  entirely  with  exact  analyses  in  other  ranges  of  these  param- 
eters but  the  few  required  constants  have  been  chosen  so  that  pre- 
dictions of  minimum  critical  Reynolds  number  are  consistent  with  re- 
sults of  more  elaborate  numerical  integration  of  the  Orr-Summerf eld 
equations  (Fig.  1). 

This  method,  which  we  have  labeled  the  "modified  Dunn-Lin  method," 
can  be  used  with  confidence  in  engineering  studies  of  laminar  boundary 
layer  control  in  water.  The  details  of  the  method  are  reported  in 
Ref.  1 and  a summary  of  the  equations  used  to  estimate  minimum  critical 
Reynolds  number  is  given  below: 

v(c) (1  - 2X (c) ) = .58 


v(c) 


X = 


<EV> 


critical 


^ 28u' (0)u(0) 

4 

c 


The  parameters  v,  X,  and  Re^  depend  on  the  characteristics  of  the 
velocity  profiles  u(n)  and  u"(n)  and  the  viscosity  u evaluated  at  the 
wall  and  the  critical  layer. 


Local  max  at  AT~40° 


Isothermal  wedge  flows,  ® 


Ref.  2 

Wedge  flows  with  heat  transfer 
in  water,  Te  = 60°F 

P = -0.1988  Ref.  3 

(3-0  Ref.  _4_ 

(3  = 1.0  Ref.  5 
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Primary  and  Secondary  Instabilities  in 
Unstably  Stratified  Shear  Flows 

by 

R.E.  Kelly  and  R.M.  Clever 
Mechanics  and  Structures  Department 
School  of  Engineering  and  Applied  Science 
University  of  California,  Los  Angeles,  Cal.  90024 

It  has  been  known  for  many  years  that  instability  in  an  un- 
stably stratified  shear  flow  (due,  say,  to  heating  from  below)  can 
occur  in  a manner  quite  dissimilar  to  that  characteristic  of  the 
homogeneous  case.  For  the  latter  case,  Tollmien-Schlichting  waves, 
which  are  two-dimensional  in  the  plane  of  the  flow,  constitute  the 
initial  instability.  For  the  unstably  stratified  case,  however, 
disturbances  can  first  appear  which  are  periodic  in  a plane  normal 
to  the  flow  and  which  exhibit  no  streamwise  waviness.  These  distur- 
bances induce  a distortion  of  the  mean  flow,  however,  and  are  there- 
fore three-dimensional  in  this  sense. 

We  call  these  buoyancy-driven  instabilities  longitudinal  "rolls" 
or  "vortices",  in  order  to  distinguish  them  from  the  hydrodynamically- 
driven  waves.  In  appearance,  they  are  similar  to  Gortler  vortices, 
and  so  it  was  appropriate  that  Gortler  wrote  the  first  review  paper 
on  their  occurrence  (Ref.  1).  This  paper  contains  most  of  the  ref- 
erences to  the  (mainly  qualitative)  experimental  work  on  this  insta- 
bility done  in  the  period  1920-1940.  More  recent  experimental  work 
concerns  the  case  of  a very  large  Prandtl  number  fluid  (Ref.  2;  the 
results  are  pertinent  to  convection  in  the  earth's  mantle),  the  ef- 
fects of  shear  upon  doubly-dif fusive  instability  (Ref.  3),  and  con- 
vection in  the  thermal-entry  flow  (Ref.  4,5).  However,  analogous 
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work  has  been  done  in  connection  with  free  convection  flows  occur- 
ring on  heated  inclined  surfaces,  where  similar  rolls  are  observed 
for  sufficiently  large  angles  of  inclination  away  from  the  vertical 
(Refs.  6-9).  The  experiments  for  fully  developed  flows  indicate 
that  the  main  effect  of  the  shear,  at  least  for  low  values  of  the 
Reynolds  number  (Re) , is  to  organize  the  convection  in  the  manner 
described  above  (i.e.,  as  spanwise  periodic  longitudinal  rolls). 

This  conclusion  agrees  with  the  theoretical  results  for  fully 

developed  flows,  e.g.,  Couette  flow  (Ref.  10-11)  and  Poiseuille 

flow  (Ref.  12).  These  analyses  indicate  that,  for  small  amounts 

of  shear,  all  disturbances  are  stabilized  by  the  shear  except  for 

the  longitudinal  roll  disturbance  which,  because  it  is  independent 

of  the  streamwise  coordinate,  cannot  be  influenced  by  the  shear. 

Thus,  it  will  become  unstable  as  soon  as  the  Rayleigh  number  (Ra) 

exceeds  the  critical  value  (Ra  = 1708  for  a flow  contained  between 

c 

two  rigid  horizontal  surfaces).  For  such  flows,  therefore,  only 
the  Rayleigh  number  is  important  for  finding  the  condition  under 
which  the  flow  becomes  unstable  due  to  buoyancy.  The  Reynolds, 
Froude,  and  Richardson  numbers  are  not  pertinent  to  predicting  the 
onset  of  the  longitudinal  convective  mode  of  instability  for  fully 
developed  flows.  Naturally,  they  are  pertinent  to  predicting  the 
onset  of  the  two-dimensional  shear  instability  which  can  occur  at 
higher  values  of  the  Reynolds  number  with  Ra  < Rac. 

Whether  or  not  a similar  conclusion  holds  for  the  case  of  an 
unstably  stratified  boundary  layer  still  is  an  unanswered  question 
(the  experiments  of  Strazisar  et  al.  (Ref.  13)  concern  a stably 
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stratified  boundary  layer) . The  theoretical  analysis  of  the  stabil- 
ity of  such  a boundary  layer  is  complicated  by  the  fact  that  non- 
parallel flow  effects  are  likely  to  play  an  important  role  in  de- 
termining the  neutral  stability  curve  of  the  longitudinal  roll 
mode.  Again,  an  analogy  with  Gortler  vortices  can  be  made,  in  which 
case  the  most  critical  disturbance  on  the  basis  of  a parallel  flow 
analysis  turns  out  to  have  an  infinite  wavelength  (see  Fig.  IX. 4 of 
Ref.  14).  This  is  in  contrast  to  shear  instabilities,  partly  be- 
cause the  longitudinal  mode  has  no  frequency  and  is  therefore  inde- 
pendent of  time  along  the  neutral  curve.  Thus,  the  parallel  flow 
stability  equation  becomes  purely  diffusive  at  the  edge  of  the 
boundary  layer.  By  including  one  nonparallel  flow  term  (correspond- 
ing to  the  vertical  velocity)  at  lowest  order  for  the  heated  in- 
clined plate  stability  problem,  Haaland  and  Sparrow  (Ref.  15)  were 
able  to  obtain  a finite  wavelength  for  the  critical  disturbance. 

Such  a simple  solution  is  unlikely  for  the  case  of  a Blasius-like 
boundary  layer  because  the  vertical  velocity  at  the  edge  of  the 
boundary  layer  is  in  a direction  opposite  (outwards)  to  that  of 
the  free  convection  boundary  layer  (inwards) . The  importance  of 
effects  associated  with  streamwise  variations  might  explain  why 
the  experimental  critical  Rayleigh  numbers  for  the  thermal  entry 
flow  as  reported  by  Kamotani  and  Sparrow  (Ref.  5)  are  so  far  above 
the  theoretical  predictions  of  Hwang  and  Cheng  (Ref.  16) , who  did 
a "local"  parallel  flow  stability  analysis.  If  such  effects  are 
important  at  lowest  order,  then  the  Reynolds  and  Prandtl  numbers 
will  become  important  in  determining  the  neutral  stability  of  the 
rolls. 
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Other  theoretical  work  on  the  problem  of  thermal  convection 
in  a shear  flow  has  concerned  the  unstably  stratified  Ekman  bound- 
ary layer  flow  (Ref.  17)  in  order  to  understand  the  "cloud  band" 
structure  often  observed  in  the  earth's  planetary  boundary  layer 
(Ref.  18) . Numerical  work  has  been  done  which  indicates  how  the 
rolls  distort  the  mean  velocity  profile  (Refs.  19-21)  and  how  the 
disturbance  field  develops  in  time  (Ref.  22) . 

At  UCLA,  our  research  has  concentrated  upon  determining  what 
happens  to  the  rolls  as  the  Rayleigh  number  increases  beyond  the 
critical  value  for  the  case  of  Couette  flow  (with  a Prandtl  number 
of  0.7).  It  is  known  for  the  case  without  shear  (Ref.  23)  that 
convective  rolls  can  themselves  become  unstable  and  lead  to  three- 
dimensional  and  possibly  unsteady  flow  as  the  Rayleigh  number  so 
increases.  For  the  case  when  shear  is  present,  the  development  of 
three-dimensionality  (in  the  sense  now  that  the  roll  becomes  wavey 
in  the  streamwise  direction)  means  that  the  perturbed  roll  can  ex- 
tract energy  from  the  mean  flow.  Although  the  rolls  occur  origi- 
nally due  to  buoyant  effects,  their  secondary  instability  can  be 
controlled  by  hydrodynamic  effects  and  might  lead  therefore  to  a 
more  severe  disruption  of  the  basic  flow  than  provided  by  the  rolls 
themselves.  In  this  sense,  the  rolls  can  be  viewed  as  introducing 
a "thermal  roughness"  into  the  flow,  with  possibly  all  the  conse- 
quences that  more  ordinary  roughness  elements  provide. 

The  analysis  is  a numerical  one,  based  on  use  of  the  Galerkin 
method,  as  discussed  in  Ref.  23.  It  consists  of  the  following  steps 
(1)  Calculate  the  disturbance  field  associated  with  the  rolls (in- 
cluding distortion  of  the  mean  velocity)  for  arbitrary  Rayleigh 
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numbers; 

(2)  perturb  the  resulting  (nonplanar)  velocity  and  temperature  dis- 
tributions by  a three-dimensional; time-dependent  disturbance  and 
find  a neutral  curve  for  the  secondary  instability  as  a function 

of  the  Rayleigh  and  Reynolds  numbers; 

(3)  calculate  the  secondary  disturbance  energy  equation  so  as  to 
understand  the  mechanism  causing  the  secondary  instability. 

The  results  will  soon  be  submitted  for  publication  and  are 

% 

not  presented  here  in  detail.  However,  they  can  be  summarized 
as  follows. 

A secondary  instability  is  possible  which  is  very  much  depen- 
dent upon  the  Reynolds  number.  As  Re  increases  to  about  100,  the 
value  of  Ra  - Ra^  at  which  the  secondary  instability  occurs  de- 
creases to  about  60,  i.e. , it  occurs  when  Ra  is  about  3.5%  above 
the  critical  value  when  Re  > 100.  The  most  unstable  secondary 
disturbance  travels  with  a wavespeed  equal  to  the  mean  of  the  basic 
flow  (for  the  antisymmetric  Couette  flow  investigated,  it  was  sta- 
tionary) . A second  mode  is  possible  which  propagates  relative  to 
the  mean  of  the  basic  flow,  but  it  occurs  only  for  larger  values 

of  Ra  — Ra  . 

c 

The  eigenvalue  calculations  indicate  therefore  that  shear  can 
destabilize  the  rolls.  The  energy  calculations  indicate  that  dis- 
tortions of  the  basic  flow,  both  in  the  plane  of  the  flow  and  nor- 
mal to  it,  act  as  a major  source  for  disturbance  energy.  Thus, 
the  idea  of  the  rolls  acting  as  a source  of  "thermal  roughness" 
seems  to  make  sense.  However,  this  does  not  imply  that  an  imme- 
diate breakdown  into  turbulence  would  occur.  An  early  observation 
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of  wavey  rolls  in  Poiseuille  flow  was  made  by  Avsec  (Ref.  24) ; 
the  rolls  are  well-defined  and  are  certainly  not  turbulent.  For 
the  case  of  free  convection  flow  in  an  inclined,  differentially 
heated  box.  Hart  (Ref.  8)  shows  that  the  rolls  must  become  quite 
distorted  before  "turbulence"  occurs  (see  his  Fig.  3;  still,  the 
turbulence  occurs  for  values  of  Ra  much  less  than  those  necessary 
for  ordinary  "turbulent"  convection  in  a horizontal  container) . 

The  situation  there  is  probably  similar  to  that  observed  in  connec- 
tion with  Taylor  vortices  (Ref.  25) , in  which  a sequence  of  states 
must  be  completed  before  the  flow  becomes  turbulent.  A similar 
sequence  probably  occurs  for  the  kind  of  unstably  stratified  shear 
flows  discussed  here.  However,  the  range  of  Rayleigh  number  neces- 
sary for  the  completion  of  the  sequence  might  be  very  much  depen- 
dent upon  the  Reynolds  number.  Thus,  the  neutral  curves  which  de- 
fine the  sequence  of  states  might  bunch  together  as  the  Reynolds 
number  increases,  allowing  transition  to  occur  for  very  modest 
Rayleigh  numbers  at  high  enough  Reynolds  numbers.  The  best  way 
to  explore  this  possibility  is  by  means  of  a carefully  controlled 
experiment. 
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The  adverse  effect  of  surface  roughness  on  a stabilized  boundary  layer  is 
well  known,  but  quantitative  data  is  lacking.  The  research  discussed  herein  is 
centered  on  obtaining  transition  data  on  a constant-pressure  body  with  wall  heating 
and  a variety  of  surface  roughnesses.  The  transition  location  will  be  determined 
by  laser  anemometer  measurements  made  in  the  NUC  low-turbulence  water  tunnel , 
capable  of  a model  length  Reynolds  numbers  of  7X10^.  Surface  roughness  of  8,  32, 

64  and  120  micro-inch  will  be  tested  with  surface  overheats  of  5 and  10°F  being 
applied.  The  purpose  of  this  paper  is  to  discuss  the  experimental  facility  being 
used  and  to  present  preliminary  results  of  transition  measurements  performed  on 
smooth  and  roughened,  9:1  ellipsoids  of  revolution  at  ambient  wall  temperature. 

It  is  important  that  the  experiments  be  carried  out  in  a well-characterized 
experimental  facility.  The  NUC  water  tunnel  has  a turbulence  level  of  less  than 
0.2%  over  the  entire  speed  range  of  0-16  m/sec.  Transition  Reynolds  numbers  of 
approximately  2 - 3 X 10^  are  observed  for  this  freestream  turbulence  level.  The 
freestream  velocity,  acoustic  and  vibration  variations  have  been  well-characterized 
in  terms  of  level  and  spectra.  The  test  section  configuration  is  an  open-jet 
type,  with  a jet  diameter  of  approximately  300  mm.  The  tunnel  is  a recirculating 
type  with  an  essentially  infinite  run  time. 

Transition  location  is  determined  by  analysis  of  laser  anemometer  measurements 
of  velocity  fluctuations  made  near  the  body  surface  at  sequential  arc  length 
locations.  Determination  of  the  freestream  velocity  for  each  location  (X/l^) 
at  which  the  turbulence  level  exceeds  some  prescribed  value  (~1Z)  defines  the 

transition  Reynolds  number,  R,  . 

X2 

A nine-to-one  ellipsoid  of  revolution  (50  mm  maximum  diameter)  was  chosen 
as  the  basic  test  body  due  to  its  previously  well-characterized  behavior  and 
its  hydrodynamic  simplicity.  The  body  has  an  essentially  zero  pressure  gradient 
history  for  a large  portion  of  its  length,  with  a minimal  region  of  adverse 
pressure  gradient  over  the  afterbody.  The  confounding  of  transition  delay  due  to 
pressure  gradient  with  the  effects  of  wall  heating  and  roughness  is  minimized. 

Two  models  have  been  constructed  thus  far.  The  first  is  the  baseline  model.  It 
is  a smooth  (surface  roughness  of  8 microinches  rms)  anodized  aluminum  body 
instrumented  with  internally-mounted  thermistors  for  the  measurement  of  the  axial 
temperature  distribution.  The  body  is  heated  via  an  axisymnetrlc  parallel-flow 
heat  exchanger  using  hot  water  pumped  in  from  outside  the  tunnel  as  the  working  fluid. 
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The  temperature  distribution  can  be  controlled  by  varying  the  heating  fluid 
mass  flow  rate,  inlet  temperature  and  flow  direction.  Figure  1 shows  the 
model  configuration.  The  second  body  is  identical  to  the  first  except  for 
material  (stainless  steel)  and  surface  roughness.  The  first  roughened  model 
has  a lath-machined  120  rms  finish  (two-  dimensional  distributed  roughness). 

The  succeeding  models  will  have  the  same  overall  configuration  but  with  64,  32 
and  16  rms  microinch  finishes. 

Test  results  obtained  thus  far  are  preliminary  in  nature  and  do  not  include 
the  effects  of  wall  heating.  Partial  results  of  the  flow  field  characterization 
measurements  are  shown  in  Figure  2,  which  also  describes  the  body  geometry.  The 
calculated  Ue/Uco  and  C as  a function  of  x are  done  at  a Re,  » 6.5X10^,  the 

" 6 2 

measured  data  points  were  taken  at  4.3X10  . Considering  the  Re  difference,  all 
measurements  are  within  experimental  error  of  the  predicted  values. 

Transition  measurements  made  thus  far  are  summarized  in  Figure  3.  The  solid 
and  dashed  curves  represent  the  transition  measurements  obtained  in  the  NASA-Ames 
and  Northrop  Aircraft  wind  tunnels  on  a similar  9:1  ellipsoidal  model.  The 
turbulence  level  for  those  measurements  was  reported  to  be  0.02  and  0.12% 
respectively.  The  data  points  for  a smooth  body  of  revolution  in  the  NUC  water  tunnel 
(turbulence  level  « 0.2%)  are  in  reasonable  agreement  with  these  data.  A transition 
value  predicted  by  an  McDonnell-Douglas  analysis  using  the  Transition  Analysis 
Program  System  (TAPS)  code  is  also  shown.  This  code  uses  a spatial  linear 
stability  analysis  to  predict  transition. 

The  limited  data  taken  on  the  roughened  model  under  the  same  freestream 
conditions  is  described  by  the  closed  data  points  in  Figure  3.  The  k/0  values 
at  the  measured  transition  point  (X/L2)  range  from  .021  to  .041  values  smaller 
than  one  might  expect  to  trigger  transition  according  to  Ko  and  Kosecoff  (1). 

If  one  selects  k/0  of  0.2  (see  ref.  (1))  as  a value  which  will  "tickle"  a 
laminar  boundary  layer  and  ultimately  result  in  transition,  the  corresponding 
axial  locations  are  3-7  mm  from  the  nose.  It  should  also  be  noted  that  the  data 
reported  here  qualitatively  agree  with  the  surface  smoothness  data  presented  by 
Smith  (2). 
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Figure  1.  Schematic  of  thermally  stabilized  boundary  layer  transition  study. 


Flowfield  characteristics  of  9: 1 ellipsoid. 


Figure  3.  Transition  data  for  smooth  and  roughened  ellipsoids. 


